Study of the anomalous acceleration of Pioneer 10 and 11 
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Our previous analyses ol radio Doppler and ranging data Irom distant spacecraft in the solar 
system indicated that an apparent anomalous acceleration is acting on Pioneer 10 and 11, with a 
magnitude ap ~ 8 x 10~* cm/s^, directed towards the Sun. Much effort has been expended looking 
for possible systematic origins of the residuals, but none has been found. A detailed investigation of 
effects both external to and internal to the spacecraft, as well as those due to modeling and compu- 
tational techniques, is provided. We also discuss the methods, theoretical models, and experimental 
techniques used to detect and study small forces acting on interplanetary spacecraft. These include 
the methods of radio Doppler data collection, data editing, and data reduction. 

There is now further data for the Pioneer 10 orbit determination. The extended Pioneer 10 data 
set spans 3 January 1987 to 22 July 1998. [For Pioneer 11 the shorter span goes from 5 January 
1987 to the time of loss of coherent data on 1 October 1990.] With these data sets and more 
detailed studies of all the systematics, we now give a result, of ap = (8.74 ± 1.33) x 10~* cm/s^. 
(Annual/diurnal variations on top of ap, that leave ap unchanged, are also reported and discussed.) 

PACS numbers: 04.80.-y, 95.10.Eg, 95.55.Pe 



I. INTRODUCTION 

Some thirty years ago, on 2 March 1972, Pioneer 10 
was launched on an Atlas/Centaur rocket from Cape 
Canaveral. Pioneer 10 was Earth's first space probe to 
an outer planet. Surviving intense radiation, it success- 
fully encountered Jupiter on 4 December 1973 In 
trail-blazing the exploration of the outer solar system. 
Pioneer 10 paved the way for, among others, Pioneer 
11 (launched on 5 April 1973), the Voyagers, Galileo, 
Ulysses, and the upcoming Cassini encounter with Sat- 
urn. After Jupiter and (for Pioneer 11) Saturn encoun- 
ters, the two spacecraft followed hyperbolic orbits near 
the plane of the ecliptic to opposite sides of the solar sys- 
tem. Pioneer 10 was also the first mission to enter the 
edge of interstellar space. That major event occurred in 
June 1983, when Pioneer 10 became the first spacecraft 
to "leave the solar system" as it passed beyond the orbit 
of the farthest known planet. 

The scientific data collected by Pioneer 10/11 has 
yielded unique information about the outer region of the 
solar system. This is due in part to the spin-stabilization 
of the Pioneer spacecraft. At launch they were spinning 
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at approximately 4.28 and 7.8 revolutions per minute 
(rpm), respectively, with the spin axes running through 
the centers of the dish antennae. Their spin-stabilizations 
and great distances from the Earth imply a minimum 
number of Earth-attitude reorientation maneuvers are re- 
quired. This permits precise acceleration estimations, to 
the level of 10~* cm/s^ (single measurement accuracy av- 
eraged over 5 days). Contrariwise, a Voyager- type three- 
axis stabilized spacecraft is not well suited for a precise 
celestial mechanics experiment as its numerous attitude- 
control maneuvers can overwhelm the signal of a small 
external acceleration. 

In summary, Pioneer spacecraft represent an ideal sys- 
tem to perform precision celestial mechanics experiments. 
It is relatively easy to model the spacecraft's behavior 
and, therefore, to study small forces affecting its motion 
in the dynamical environment of the solar system. In- 
deed, one of the main objectives of the Pioneer extended 
missions (post Jupiter/Saturn encounters) ^ was to per- 
form accurate celestial mechanics experiments. For in- 
stance, an attempt was made to detect the presence of 
small bodies in the solar system, primarily in the Kuiper 
belt. It was hoped that a small perturbation of the 
spacecraft's trajectory would reveal the presence of these 
objects 0-0 ■ Furthermore, due to extremely precise 
navigation and a high quality tracking data, the Pioneer 
10 scientific program also included a search for low fre- 
quency gravitational waves [lol [Tl| . 

Beginning in 1980, when at a distance of 20 astronom- 
ical units (AU) from the Sun the solar-radiation-pressure 
acceleration on Pioneer 10 away from the Sun had de- 
creased to < 5 X 10~^ cm/s^, we found that the largest 
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systematic error in the acceleration residuals was a con- 
stant bias, ap, directed toward the Sun. Such anoma- 
lous data have been continuously received ever since. Jet 
Propulsion Laboratory (JPL) and The Aerospace Corpo- 
ration produced independent orbit determination analy- 
ses of the Pioneer data extending up to July 1998. We 
ultimately concluded that there is an unmodeled 

acceleration, ap, towards the Sun of ~ 8 x 10~^ cm/s^ 
for both Pioneer 10 and Pioneer 11. 

The purpose of this paper is to present a detailed expla- 
nation of the analysis of the apparent anomalous, weak, 
long-range acceleration of the Pioneer spacecraft that we 
detected in the outer regions of the solar system. We at- 
tempt to survey all sensible forces and to estimate their 
contributions to the anomalous acceleration. We will dis- 
cuss the effects of these small non-gravitational forces 
(both generated on-board and external to the vehicle) on 
the motion of the distant spacecraft together with the 
methods used to collect and process the radio Doppler 
navigational data. 

We begin with descriptions of the spacecraft and other 
systems and the strategies for obtaining and analyzing 
information from them. In Section ^ we describe the 
Pioneer (and other) spacecraft. We provide the reader 
with important technical information on the spacecraft, 
much of which is not easily accessible. In Section ITTll we 
describe how raw data is obtained and analyzed and in 
Section |W1 we discuss the basic elements of a theoretical 
foundation for spacecraft navigation in the solar system. 

The next major part of this manuscript is a description 
and analysis of the results of this investigation. We first 
describe how the anomalous acceleration was originally 
identified from the data of all the spacecraft in Section 
IVl 12, 13]. We then give our recent results in Section 
IVII In the following three sections we discuss possible 
experimental systematic origins for the signal. These in- 
clude systematics generated by physical phenomena from 
sources external to (Section IVIIII and internal to (Sec- 
tion lVIli|) the spacecraft. This is followed by Section llXl 
where the accuracy of the solution for ap is discussed. In 
the process we go over possible numerical/calculational 
errors/systematics. Sections lVIIIIXI are then summarized 
in the total error budget of Section Ixl 

We end our presentation by first considering possible 
unexpected physical origins for the anomaly fSection lXI|l . 
In our conclusion. Section [XIII we summarize our results 
and suggest venues for further study of the discovered 
anomaly. 



II. THE PIONEER AND OTHER SPACECRAFT 

In this section we describe in some detail the Pioneer 
10 and 11 spacecraft and their missions. We concentrate 
on those spacecraft systems that play important roles in 
maintaining the continued function of the vehicles and 
in determining their dynamical behavior in the solar sys- 
tem. Specifically we present an overview of propulsion 



and attitude control systems, as well as thermal and com- 
munication systems. 

Since our analysis addresses certain results from the 
Galileo and Ulysses missions, we also give short descrip- 
tions of these missions in the final subsection. 



A. General description of the Pioneer spacecraft 

Although some of the more precise details are often 
difficult to uncover, the general parameters of the Pi- 
oneer spacecraft are known and well documented 
The two spacecraft are identical in design . At launch 
each had a "weight" (mass) of 259 kg. The "dry weight" 
of the total module was 223 kg as there were 36 kg of 
hydrazine propellant 0, 0| . The spacecraft were de- 
signed to fit within the three meter diameter shroud of 
an added third stage to the Atlas/Centaur launch vehi- 
cle. Each spacecraft is 2.9 m long from its base to its 
cone-shaped medium-gain antenna. The high gain an- 
tenna (HGA) is made of aluminum honeycomb sandwich 
material. It is 2.74 m in diameter and 46 cm deep in the 
shape of a parabolic dish. (See Figures Q and |21) 




FIG. 1: NASA photo #72HC94, with caption "The Pioneer 
F spacecraft during a checkout with the launch vehicle third 
stage at Cape Kennedy." Pioneer F became Pioneer 10. 

The main equipment compartment is 36 cm deep. 
The hexagonal flat top and bottom have 71 cm long 
sides. The equipment compartment provides a thermally 
controlled environment for scientific instruments. Two 
three-rod trusses, 120 degrees apart, project from two 
sides of the equipment compartment. At their ends, each 
holds two SNAP-19 (Space Nuclear Auxiliary Power, 
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FIG. 2: A drawing of the Pioneer spacecraft. 



model 19) RTGs (Radioisotope Thermoelectric Gener- 
ators) built by Teledyne Isotopes for the Atomic En- 
ergy Commission. These RTGs are situated about 3 m 
from the center of the spacecraft and generate its electric 
power. [We will go into more detail on the RTGs in Sec- 
tion lVIlil ] A third single-rod boom, 120 degrees from the 
other two, positions a magnetometer about 6.6 m from 
the spacecraft's center. All three booms were extended 
after launch. With the mass of the magnetometer being 
5 kg and the mass of each of the four RTGs being 13.6 
kg, this configuration defines the main moment of inertia 
along the z-spin-axis. It is about Iz « 588.3 kg m^ [l^ . 
[Observe that this all left only about 164 kg for the main 
bus and superstructure, including the antenna.] 

Figures ^ and |21 show the arrangement within the 
spacecraft equipment compartment. The majority of the 
spacecraft electrical assemblies are located in the cen- 
tral hexagonal portion of the compartment, surrounding 
a 16.5-inch-diameter spherical hydrazine tank. Most of 
the scientific instruments' electronic units and internally- 
mounted sensors are in an instrument bay ("squashed" 
hexagon) mounted on one side of the central hexagon. 
The equipment compartment is in an aluminum honey- 
comb structure. This provides support and meteoroid 
protection. It is covered with insulation which, together 
with louvers under the platform, provides passive thermal 
control. [An exception is from off-on control by thermal 
power dissipation of some subsystems. (See Sec. IVIIIII .] 



B. Propulsion and attitude control systems 

Three pairs of these rocket thrusters near the rim of the 
HGA provide a threefold function of spin-axis precession, 
mid-course trajectory correction, and spin control. Each 
of the three thruster pairs develops its repulsive jet force 
from a catalytic decomposition of liquid hydrazine in a 
small rocket thrust chamber attached to the oppositely- 
directed nozzle. The resulted hot gas is then expended 
through six individually controlled thruster nozzles to ef- 
fect spacecraft maneuvers. 

The spacecraft is attitude-stabilized by spinning about 
an axis which is parallel to the axis of the HGA. The 
nominal spin rate for Pioneer 10 is 4.8 rpm. Pioneer 11 
spins at approximately 7.8 rpm because a spin-controlling 
thruster malfunctioned during the spin-down shortly af- 
ter launch. [Because of the danger that the thruster's 
valve would not be able to close again, this particular 
thruster has not been used since.] During the mission 
an Earth-pointing attitude is required to illuminate the 
Earth with the narrow-beam HGA. Periodic attitude ad- 
justments are required throughout the mission to com- 
pensate for the variation in the heliocentric longitude 
of the Earth-spacecraft line. [In addition, correction of 
launch vehicle injection errors were required to provide 
the desired Jupiter encounter trajectory and Saturn (for 
Pioneer 11) encounter trajectory.] These velocity vector 
adjustments involved reorienting the spacecraft to direct 
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the thrust in the desired direction. 

There were no anomahes in the engineering telemetry 
from the propulsion system, for either spacecraft, dur- 
ing any mission phase from launch to termination of the 
Pioneer mission in March 1997. From the viewpoint of 
mission operations at the NASA/ Ames control center, 
the propulsion system performed as expected, with no 
catastrophic or long-term pressure drops in the propul- 
sion tank. Except for the above-mentioned Pioneer 11 
spin-thruster incident, there was no malfunction of the 
propulsion nozzles, which were only opened every few 
months by ground command. The fact that pressure was 
maintained in the tank has been used to infer that no 
impacts by Kuiper belt objects occurred, and a limit has 
been placed on the size and density distribution of such 
objects 0, another useful scientific result. 

For attitude control, a star sensor (referenced to Cano- 
pus) and two sunlight sensors provided reference for ori- 
entation and roll maneuvers. The star sensor on Pioneer 
10 became inoperative at Jupiter encounter, so the sun 
sensors were used after that. For Pioneer 10, spin calibra- 
tion was done by the DSN until 17 July 1990. From 1990 
to 1993 determinations were made by analysts using data 
from the Imaging Photo Polarimeter (IPP). After the 6 
July 1993 maneuver, there was not enough power left to 
support the IPP. But approximately every six months 
analysts still could get a rough determination using in- 
formation obtained from conscan maneuvers jlSj on an 
uplink signal. When using conscan, the high gain feed 
is off-set. Thrustcr firings are used to spiral in to the 
correct pointing of the spacecraft antenna to give the 
maximum signal strength. To run this procedure (con- 
scan and attitude) it is now necessary to turn off the 
traveling-wave-tube (TWT) amplifier. So far, the power 
and tube life-cycle have worked and the Jet Propulsion 
Laboratory's (JPL) Deep Space Network (DSN) has been 
able to reacquire the signal. It takes about 15 minutes 
or so to do a maneuver. [The magnetometer boom in- 
corporates a hinged, viscous, damping mechanism at its 
attachment point, for passive nutation control.] 

In the extended mission phase, after Jupiter and Sat- 
urn encounters, the thrusters have been used for preces- 
sion maneuvers only. Two pairs of thrusters at opposite 
sides of the spacecraft have nozzles directed along the 
spin axis, fore and aft (See Figure|21) In precession mode, 
the thrusters are fired by opening one nozzle in each pair. 
One fires to the front and the other fires to the rear of the 
spacecraft in brief thrust pulses. Each thrust pulse 
precesses the spin axis a few tenths of a degree until the 
desired attitude is reached. 

The two nozzles of the third thruster pair, no longer 
in use, are aligned tangentially to the antenna rim. One 
points in the direction opposite to its (rotating) velocity 
vector and the other with it. These were used for spin 
control. 



C. Thermal system and on-board power 

Early on the spacecraft instrument compartment is 
thermally controlled between « F and 90 F. This is 
done with the aid of thermo-responsive louvers located at 
the bottom of the equipment compartment. These lou- 
vers are adjusted by bi-metallic springs. They are com- 
pletely closed below ~ 40 F and completely open above ~ 
85 F. This allows controlled heat to escape in the equip- 
ment compartment. Equipment is kept within an opera- 
tional range of temperatures by multi-layered blankets of 
insulating aluminum plastic. Heat is provided by electric 
heaters, the heat from the instruments themselves, and 
by twelve one- watt radioisotope heaters powered directly 
by non- fissionable plutonium (^||Pu^^g|U-|-|He). 

^•^^Pu, with a half life time of 87.74 years, also provides 
the thermal source for the thermoelectric devices in the 
RTGs. Before launch, each spacecraft's four RTGs deliv- 
ered a total of approximately 160 W of electrical power 
mini. Each of the four space-proven SNAP-19 RTGs 
converts 5 to 6 percent of the heat released from pluto- 
nium dioxide fuel to electric power. RTG power is great- 
est at 4.2 Volts; an inverter boosts this to 28 Volts for 
distribution. RTG life is degraded at low currents; there- 
fore, voltage is regulated by shunt dissipation of excess 
power. 

The power subsystem controls and regulates the RTG 
power output with shunts, supports the spacecraft load, 
and performs battery load-sharing. The silver cadmium 
battery consists of eight cells of 5 ampere-hours capacity 
each. It supplies pulse loads in excess of RTG capability 
and may be used for sharing peak loads. The battery 
voltage is often discharged and charged. This can be 
seen by telemetry of the battery discharge current and 
charge current 

At launch each RTG supplied about 40 W to the input 
of the ^ 4.2 V Inverter Assemblies. (The output for other 
uses includes the DC bus at 28 V and the AC bus at 61 
V) Even though electrical power degrades with time (see 
Section lVlIIDI) . at —41 F the essential platform temper- 
ature as of the year 2000 is still between the acceptable 
limits of -63 F to 180 F. The RF power output from the 
traveling-wave-tube amplifier is still operating normally. 

The equipment compartment is insulated from extreme 
heat influx with aluminized mylar and kapton blankets. 
Adequate warmth is provided by dissipation of 70 to 120 
watts of electrical power by electronic units within the 
compartment; louvers regulating the release of this heat 
below the mounting platform maintain temperatures in 
the vicinity of the spacecraft equipment and scientific in- 
struments within operating limits. External component 
temperatures are controlled, where necessary, by appro- 
priate coating and, in some cases, by radioisotope or elec- 
trical heaters. 

The energy production from the radioactive decay 
obeys an exponential law. Hence, 29 years after launch, 
the radiation from Pioneer lO's RTGs was about 80 per- 
cent of its original intensity. However the electrical power 
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delivered to the equipment compartment has decayed at 
a faster rate than the ^'^*Pu decays radioactively. Specif- 
ically, the electrical power first decayed very quickly and 
then slowed to a still fast linear decay By 1987 the 
degradation rate was about —2.6 W/yr for Pioneer 10 
and even greater for the sister spacecraft. 

This fast depletion rate of electrical power from the 
RTGs is caused by normal deterioration of the thermo- 
couple junctions in the thermoelectric devices. 

The spacecraft needs 100 W to power all systems, in- 
cluding 26 W for the science instruments. Previously, 
when the available electrical power was greater than 100 
W, the excess power was either thermally radiated into 
space by a shunt-resistor radiator or it was used to charge 
a battery in the equipment compartment. 

At present only about 65 W of power is available to 
Pioneer 10 Therefore, all the instruments are no 

longer able to operate simultaneously. But the power 
subsystem continues to provide sufhcient power to sup- 
port the current spacecraft load: transmitter, receiver, 
command and data handling, and the Geiger Tube Tele- 
scope (GTT) science instrument. As pointed out in Sec. 
Ill El the science package and transmitter are turned off 
in extended cruise mode to provide enough power to fire 
the attitude control thrusters. 



D. Communication system 

The Pioneer 10/11 communication systems use S-band 
(A ~ 13 cm) Doppler frequencies |24| . The communica- 
tion uplink from Earth is at approximately 2.11 GHz. 
The two spacecraft transmit continuously at a power of 
eight watts. They beam their signals, of approximate fre- 
quency 2.29 GHz, to Earth by means of the parabolic 2.74 
m high-gain antenna. Phase coherency with the ground 
transmitters, referenced to H-maser frequency standards, 
is maintained by means of an S-band transponder with 
the 240/221 frequency turnaround ratio (as indicated by 
the values of the above mentioned frequencies). 

The communications subsystem provides for: i) up- 
link and down-link communications; ii) Doppler coher- 
ence of the down-link carrier signal; and iii) generation 
of the conscan signal for closed loop precession of 
the spacecraft spin axis towards Earth. S-band carrier 
frequencies, compatible with DSN, are used in conjunc- 
tion with a telemetry modulation of the down-link signal. 
The high-gain antenna is used to maximize the teleme- 
try data rate at extreme ranges. The coupled medium- 
gain/omni-directional antenna with fore and aft elements 
respectively, provided broad-angle communications at in- 
termediate and short ranges. For DSN acquisition, these 
three antennae radiate a non-coherent RF signal, and for 
Doppler tracking, there is a phase coherent mode with a 
frequency translation ratio of 240/221. 

Two frequency-addressable phase-lock receivers are 
connected to the two antenna systems through a ground- 
commanded transfer switch and two diplexers, provid- 



ing access to the spacecraft via either signal path. The 
receivers and antennae are interchangeable through the 
transfer switch by ground command or automatically, if 
needed. 

There is a redundancy in the communication systems, 
with two receivers and two transmitters coupled to two 
traveling-wave-tube amplifiers. Only one of the two re- 
dundant systems has been used for the extended mis- 
sions, however. 

At launch, communication with the spacecraft was at 
a data rate 256 bps for Pioneer 10 (1024 bps for Pioneer 
11). Data rate degradation has been —1.27 mbps/day for 
Pioneer 10 (-8.78 mbps/day for Pioneer 11). The DSN 
still continues to provide good data with the received sig- 
nal strength of about —178 dBm (only a few dB from the 
receiver threshold). The data signal to noise ratio is still 
mainly under 0.5 dB. The data deletion rate is often be- 
tween and 50 percent, at times more. However, during 
the test of 11 March 2000, the average deletion rate was 
about 8 percent. So, quality data are still available. 



E. Status of tlie extended mission 

The Pioneer 10 mission officially ended on 31 March 
1997 when it was at a distance of 67 AU from the Sun. 
(See FigureOl) At a now nearly constant velocity relative 
to the Sun of ~12.2 km/s. Pioneer 10 will continue its 
motion into interstellar space, heading generally for the 
red star Aldebaran, which forms the eye of Taurus (The 
Bull) Constellation. Aldebaran is about 68 light years 
away and it would be expected to take Pioneer 10 over 2 
million years to reach its neighborhood. 




FIG. 3: Ecliptic pole view of Pioneer 10, Pioneer 11, and 
Voyager trajectories. Pioneer 11 is traveling approximately 
in the direction of the Sun's orbital motion about the galactic 
center. The galactic center is approximately in the direction 
of the top of the figure. [Digital artwork by T. Esposito. 
NASA ARC Image # AC97-0036-3.] 
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A switch failure in the Pioneer 11 radio system on 1 
October 1990 disabled the generation of coherent Doppler 
signals. So, after that date, when the spacecraft was ~ 30 
AU away from the Sun, no useful data have been gen- 
erated for our scientific investigation. Furthermore, by 
September 1995, its power source was nearly exhausted. 
Pioneer 11 could no longer make any scientific observa- 
tions, and routine mission operations were terminated. 
The last communication from Pioneer 11 was received 
in November 1995, when the spacecraft was at distance 
of ^ 40 AU from the Sun. (The relative Earth motion 
carried it out of view of the spacecraft antenna.) The 
spacecraft is headed toward the constellation of Aquila 
(The Eagle), northwest of the constellation of Sagittarius, 
with a velocity relative to the Sun of ^11.6 km/s Pioneer 
11 should pass close to the nearest star in the constella- 
tion Aquila in about 4 million years ^ . (Pioneer 10 and 
11 orbital parameters are given in the Appendix.) 

However, after mission termination the Pioneer 10 ra- 
dio system was still operating in the coherent mode when 
commanded to do so from the Pioneer Mission Opera- 
tions center at the NASA Ames Research Center (ARC). 
As a result, after 31 March 1997, JPL's DSN was stiU able 
to deliver high-quality coherent data to us on a regular 
schedule from distances beyond 67 AU. 

Recently, support of the Pioneer spacecraft has been on 
a non-interference basis to other NASA projects. It was 
used for the purpose of training Lunar Prospector con- 
trollers in DSN coordination of tracking activities. Un- 
der this training program, ARC has been able to main- 
tain contact with Pioneer 10. This has required careful 
attention to the DSN's ground system, including the in- 
stallation of advanced instrumentation, such as low-noisc 
digital receivers. This extended the lifetime of Pioneer 
10 to the present. [Note that the DSN's early estimates, 
based on instrumentation in place in 1976, predicted that 
radio contact would be lost about 1980.] 

At the present time it is mainly the drift of the space- 
craft relative to the solar velocity that necessitates ma- 
neuvers to continue keeping Pioneer 10 pointed towards 
the Earth. The latest successful precession maneuver to 
point the spacecraft to Earth was accomplished on 11 
February 2000, when Pioneer 10 was at a distance from 
the Sun of 75 AU. [The distance from the Earth was ~ 76 
AU with a corresponding round-trip light time of about 
21 hour.] The signal level increased 0.5-0.75 dBm 0] as 
a result of the maneuver. 

This was the seventh successful maneuver that has 
been done in the blind since 26 January 1997. At that 
time it had been determined that the electrical power to 
the spacecraft had degraded to the point where the space- 
craft transmitter had to be turned off to have enough 
power to perform the maneuver. After 90 minutes in 
the blind the transmitter was turned back on again. So, 
despite the continued weakening of Pioneer lO's signal, 
radio Doppler measurements were still available. The 
next attempt at a maneuver, on 8 July 2000, turned out 
in the end to be successful. Signal was tracked on 9 July 



2001. Contact was reestablished on the 30th anniversary 
of launch, 2 March 2002. 



F. The Galileo and Ulysses missions and spacecraft 

1. The Galileo mission 

The Galileo mission to explore the Jovian system 
was launched 18 October 1989 aboard the Space Shuttle 
Discovery. Due to insufficient launch power to reach its 
final destination at 5.2 AU, a trajectory was chosen with 
planetary flybys to gain gravity assists. The spacecraft 
flew by Venus on 10 February 1990 and twice by the 
Earth, on 8 December 1990 and on 8 December 1992. The 
current Galileo Millennium Mission continues to study 
Jupiter and its moons, and coordinated observations with 
the Cassini flyby in December 2000. 

The dynamical properties of the Galileo spacecraft are 
very well known. At launch the orbiter had a mass of 
2,223 kg. This included 925 kg of usable propellant, 
meaning over 40% of the orbiter's mass at launch was 
for propellant! The science payload was 118 kg and the 
probe's total mass was 339 kg. Of this latter, the probe 
descent module was 121 kg, including a 30 kg science 
payload. The tensor of inertia of the spacecraft had the 
following components at launch: Jxx = 4454.7, Jyj = 
4061.2, J^^ = 5967.6, Jxy = -52.9, J^z = 3.21, Jy^ = 
— 15.94 in units of kg m^. Based on the area of the 
sun-shade plus the booms and the RTGs we obtained 
a maximal cross-sectional area of 19.5 m^. Each of the 
two of the Gahleo's RTGs at launch delivered of 285 W 
of electric power to the subsystems. 

Unlike previous planetary spacecraft, Galileo featured 
an innovative "dual spin" design: part of the orbiter 
would rotate constantly at about three rpm and part of 
the spacecraft would remain fixed in (solar system) in- 
ertial space. This means that the orbiter could easily 
accommodate magnetospheric experiments (which need 
to made while the spacecraft is sweeping) while also pro- 
viding stability and a fixed orientation for cameras and 
other sensors. The spin rate could be increased to 10 rev- 
olutions per minute for additional stability during major 
propulsive maneuvers. 

Apparently there was a mechanical problem between 
the spinning and non-spinning sections. Because of this, 
the project decided to often use an all-spinning mode, 
of about 3.15 rpm. This was especially true close to the 
Jupiter Orbit Insertion ( JOI) , when the entire spacecraft 
was spinning (with a slower rate, of course). 

Galileo's original design called for a deployable high- 
gain antenna (HGA) to unfurl. It would provide approx- 
imately 34 dB of gain at X-band (10 GHz) for a 134 kbps 
downlink of science and priority engineering data. How- 
ever, the X-band HGA failed to unfurl on 11 April 1991. 
When it again did not deploy following the Earth fly- 
by in 1992, the spacecraft was reconfigured to utilize the 
S-band, 8 dB, omni-directional low-gain antenna (LGA) 
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for downlink. 

The S-band frequencies are 2.113 GHz - up and 
2.295 GHz - down, a conversion factor of 240/221 at 
the Doppler frequency transponder. This configuration 
yielded much lower data rates than originally scheduled, 
8-16 bps through JOI 27). Enhancements at the DSN 
and reprogramming the flight computers on Galileo in- 
creased telemetry bit rate to 8-160 bps, starting in the 
spring of 1996. 

Currently, two types of Galileo navigation data are 
available, namely Doppler and range measurements. As 
mentioned before, an instantaneous comparison between 
the ranging signal that goes up with the ranging signal 
that comes down would yield an "instantaneous" two- 
way range delay. Unfortunately, an instantaneous com- 
parison was not possible in this case. The reason is that 
the signal-to-noise ratio on the incoming ranging signal 
is small and a long integration time (typically minutes) 
must be used (for correlation purposes). During such 
long integration times, the range to the spacecraft is con- 
stantly changing. It is therefore necessary to "electron- 
ically freeze" the range delay long enough to permit an 
integration to be performed. The result represents the 
range at the moment of freezing |2^ |2^ . 



2. The Ulysses mission 

Ulysses was launched on 6 October 1990, also from 
the Space Shuttle Discovery, as a cooperative project of 
NASA and the European Space Agency (ESA) . JPL man- 
ages the US portion of the mission for NASA's Office of 
Space Science. Ulysses' objective was to characterize the 
heliosphere as a function of solar latitude [s^l ■ To reach 
high solar latitudes, its voyage took it to Jupiter on 8 
February 1992. As a result, its orbit plane was rotated 
about 80 degrees out of the ecliptic plane. 

Ulysses explored the heliosphere over the Sun's south 
pole between June and November, 1994, reaching maxi- 
mum Southern latitude of 80.2 degrees on 13 September 
1994. It continued in its orbit out of the plane of the 
ecliptic, passing perihelion in March 1995 and over the 
north solar pole between June and September 1995. It 
returned again to the Sun's south polar region in late 
2000. 

The total mass at launch was the sum of two parts: a 
dry mass of 333.5 kg plus a propellant mass of 33.5 kg. 
The tensor of inertia is given by its principal components 
J^^ = 371.62, Jyy = 205.51, J^z = 534.98 in units kg m^. 
The maximal cross section is estimated to be 10.056 m^. 
This estimation is based on the radius of the antenna 1.65 
m (8.556 m^) plus the areas of the RTGs and part of the 
science compartment (yielding an additional ~ 1.5 m^). 
The spacecraft was spin-stabilized at 4.996 rpm. The 
electrical power is generated by modern RTGs, which 
are located much closer to the main bus than are those 
of the Pioneers. The power generated at launch was 285 
W. 



Communications with the spacecraft are performed at 
X-band (for downlink at 20 W with a conversion factor 
of 880/221) and S-band (both for uplink 2111.607 MHz 
and downlink 2293.148 MHz, at 5 W with a conversion 
factor of 240/221). Currently both Doppler and range 
data are available for both frequency bands. While the 
main communication link is S-up/X-down, the S-down 
link was used only for radio-science purposes. 

Because of Ulysses' closeness to the Sun and also be- 
cause of its construction, any hope to model Ulysses for 
small forces might appear to be doomed by solar radia- 
tion pressure and internal heat radiation from the RTGs. 
However, because the Doppler signal direction is towards 
the Earth while the radiation pressure varies with dis- 
tance and has a direction parallel the Sun-Ulysses line, 
in principle these effects could be separated. And again, 
there was range data. This all would make it easier to 
model non-gravitational acceleration components normal 
to the line of sight, which usually are poorly and not sig- 
nificantly determined. 

The Ulysses spacecraft spins at ^ 5 rpm around its 
antenna axis (4.996 rpm initially). The angle of the spin 
axis with respect to the spacecraft- Sun line varies from 
near zero at Jupiter to near 50 degrees at perihelion. 
Any on-board forces that could perturb the spacecraft 
trajectory are restricted to a direction along the spin axis. 
[The other two components are canceled out by the spin.] 

As the spacecraft and the Earth travel around the Sun, 
the direction from the spacecraft to the Earth changes 
continuously. Regular changes of the attitude of the 
spacecraft are performed throughout the mission to keep 
the Earth within the narrow beam of about one degree 
full width of the spacecraft-fixed parabolic antenna. 



III. DATA ACQUISITION AND PREPARATION 

Discussions of radio-science experiments with space- 
craft in the solar system requires at least a general knowl- 
edge of the sophisticated experimental techniques used at 
the DSN complex. Since its beginning in 1958 the DSN 
complex has undergone a number of major upgrades and 
additions. This was necessitated by the needs of particu- 
lar space missions. [The last such upgrade was conducted 
for the Cassini mission when the DSN capabilities were 
extended to cover the Ka radio frequency bandwidth. 
For more information on DSN methods, techniques, and 
present capabihtics, see 0.] For the purposes of the 
present analysis one will need a general knowledge of the 
methods and techniques implemented in the radio-science 
subsystem of the DSN complex. 

This section reviews the techniques that are used to 
obtain the radio tracking data from which, after analysis, 
results are generated. Here we will briefly discuss the 
DSN hardware that plays a pivotal role for our study of 
the anomalous acceleration. 
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A. Data acquisition 

The Deep Space Network (DSN) is the network of 
ground stations that are employed to track interplanetary 
spacecraft |^ [S^l- There are three ground DSN com- 
plexes, at Goldstone, California, at Robledo de Chavela, 
outside Madrid, Spain, and at Tidbinbilla, outside Can- 
berra, Australia. 

There are many antennae, both existing and decom- 
missioned, that have been used by the DSN for space- 
craft navigation. For our four spacecraft (Pioneer 10, 
11, Gahleo, and Ulysses), depending on the time period 
involved, the following Deep Space Station (DSS) anten- 
nae were among those used: (DSS 12, 14, 24) at the 
California antenna complex; (DSS 42, 43, 45, 46) at the 
Australia complex; and (DSS 54, 61, 62, 63) at the Spain 
complex. Specifically, the Pioneers used (DSS 12, 14, 42, 
43, 62, 63), Galileo used (DSS 12, 14, 42, 43, 63), and 
Ulysses used (DSS 12, 14, 24, 42, 43, 46, 54, 61, 63). 

The DSN tracking system is a phase coherent system. 
By this we mean that an "exact" ratio exists between the 
transmission and reception frequencies; i.e., 240/221 for 
S-band or 880/221 for X-band (This is in distinction 
to the usual concept of coherent radiation used in atomic 
and astrophysics.) 

Frequency is an average frequency, defined as the num- 
ber of cycles per unit time. Thus, accumulated phase is 
the integral of frequency. High measurement precision 
is attained by maintaining the frequency accuracy to 1 
part per 10^^ or better (This is in agreement with the 
expected Allan deviation for the S-band signals.) 

The DSN Frequency and Timing System (FTS): 
The dsn's FTS is the source for the high accuracy just 
mentioned (see Figure 0J. At its center is an hydrogen 
maser that produces a precise and stable reference fre- 
quency [33I [3^ . These devices have Allan deviations IsBI 
of approximately 3 x 10~^^ to 1 x 10~^^ for integration 
times of 10^ to 10'^ seconds, respectively. 
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FIG. 4: Block-diagram of the DSN complex as used for radio 
Doppler tracking of an interplanetary spacecraft. For more 
detailed drawings and technical specifications see Ref. [3l|. 



These masers are good enough so that the quality 
of Doppler-measurement data is limited by thermal or 
plasma noise, and not by the inherent instability of the 
frequency references. Due to the extreme accuracy of 
the hydrogen masers, one can very precisely characterize 
the spacecraft's dynamical variables using Doppler and 
range techniques. The FTS generates a 5 MHz and 10 
MHz reference frequency which is sent through the lo- 
cal area network to the Digitally Controlled Oscillator 
(DCO). 

The Digitally Controlled Oscillator (DCO) and 
Exciter: Using the highly stable output from the FTS, 
the DCO, through digitally controlled frequency multi- 
pliers, generates the Track Synthesizer Frequency (TSF) 
of ~ 22 MHz. This is then sent to the Exciter Assembly. 
The Exciter Assembly multiplies the TSF by 96 to pro- 
duce the S-band carrier signal at ~ 2.2 GHz. The signal 
power is amplified by Traveling Wave Tubes (TWT) for 
transmission. If ranging data are required, the Exciter 
Assembly adds the ranging modulation to the carrier. 
[The DSN tracking system has undergone many upgrades 
during the 29 years of tracking Pioneer 10. During this 
period internal frequencies have changed.] 

This S-band frequency is sent to the antenna where 
it is amplified and transmitted to the spacecraft. The 
onboard receiver tracks the up-link carrier using a phase 
lock loop. To ensure that the reception signal does not 
interfere with the transmission, the spacecraft (e.g.. Pi- 
oneer) has a turnaround transponder with a ratio of 
240/221. The spacecraft transmitter's local oscillator is 
phase locked to the up-link carrier. It multiplies the re- 
ceived frequency by the above ratio and then re-transmits 
the signal to Earth. 

Receiver and Doppler Extractor: When the two- 
way 36] signal reaches the ground, the receiver locks on 
to the signal and tunes the Voltage Control Oscillator 
(VCO) to null out the phase error. The signal is sent 
to the Doppler Extractor. At the Doppler Extractor the 
current transmitter signal from the Exciter is multiplied 
by 240/221 (or 880/241 for X-band)) and a bias, of 1 
MHz for S-band or 5 MHz for X-band H, is added to 
the Doppler. The Doppler data is no longer modulated 
at S-band but has been reduced as a consequence of the 
bias to an intermediate frequency of 1 or 5 MHz 

Since the light travel time to and from Pioneer 10 is 
long (more than 20 hours), the transmitted frequency 
and the current transmitted frequency can be different. 
The difference in frequencies are recorded separately and 
are accounted for in the orbit determination programs we 
discuss in Section IVI 

Metric Data Assembly (MDA): The MDA con- 
sists of computers and Doppler counters where continu- 
ous count Doppler data are generated. The intermediate 
frequency (IF) of 1 or 5 MHz with a Doppler modulation 
is sent to the Metric Data Assembly (MDA). From the 
FTS a 10 pulse per second signal is also sent to the MDA 
for timing. At the MDA, the IF and the resulting Doppler 
pulses are counted at a rate of 10 pulses per second. At 
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each tenth of a second, the number of Doppler pulses 
are counted. A second counter begins at the instant the 
first counter stops. The result is continuously-counted 
Doppler data. (The Doppler data is a biased Doppler 
of 1 MHz, the bias later being removed by the analyst 
to obtain the true Doppler counts.) The Range data (if 
present) together with the Doppler data is sent separately 
to the Ranging Demodulation Assembly. The accompa- 
nying Doppler data is used to rate aid (i.e., to "freeze" 
the range signal) for demodulation and cross correlation. 

Data Communication: The total set of tracking 
data is sent by local area network to the communica- 
tion center. From there it is transmitted to the Goddard 
Communication Facility via commercial phone lines or by 
government leased lines. It then goes to JPL's Ground 
Communication Facility where it is received and recorded 
by the Data Records Subsystem. 



B. Radio Doppler and range techniques 

Various radio tracking strategies are available for de- 
termining the trajectory parameters of interplanetary 
spacecraft. However, radio tracking Doppler and range 
techniques are the most commonly used methods for nav- 
igational purposes. The position and velocities of the 
DSN tracking stations must be known to high accuracy. 
The transformation from a Earth fixed coordinate sys- 
tem to the International Earth Rotation Service (lERS) 
Celestial System is a complex series of rotations that in- 
cludes precession, nutation, variations in the Earth's ro- 
tation (UTl-UTC) and polar motion. 

Calculations of the motion of a spacecraft are made 
on the basis of the range time-delay and/or the Doppler 
shift in the signals. This type of data was used to deter- 
mine the positions, the velocities, and the magnitudes of 
the orientation maneuvers for the Pioneer, Galileo, and 
Ulysses spacecraft considered in this study. 

Theoretical modeling of the group delays and phase de- 
lay rates are done with the orbit determination software 
we describe in the next section. 

Data types: Our data describes the observations that 
are the basis of the results of this paper. We receive 
our data from DSN in closed- loop mode, i.e., data that 
has been tracked with phase lock loop hardware. (Open 
loop data is tape recorded but not tracked by phase lock 
loop hardware.) The closed-loop data constitutes our 
Archival Tracking Data File (ATDF), which we copy H^l 
to the National Space Science Data Center (NSSDC) on 
magnetic tape. The ATDF files are stored on hard disk 
in the RMDC (Radio Metric Data Conditioning group) 
of JPL's Navigation and Mission Design Section. We 
access these files and run standard software to produce 
an Orbit Data File for input into the orbit determination 
programs which we use. (See Section Ivl) 

The data types are two-way and three-way [s^ Doppler 
and two-way range. (Doppler and range are defined in 
the following two subsections.) Due to unknown clock 



offsets between the stations, three-way range is generally 
not taken or used. 

The Pioneer spacecraft only have two- and three-way 
S-band Doppler. Galileo also has S-band range data 
near the Earth. Ulysses has two- and three-way S-band 
up- link and X-band down-link Doppler and range as 
well as S-band up-link and S-band down-link, although 
we have only processed the Ulysses S-band up-link and 
X-band down-link Doppler and range. 



1. Doppler experimental techniques and strategy 

In Doppler experiments a radio signal transmit- 
ted from the Earth to the spacecraft is coherently 
transponded and sent back to the Earth. Its frequency 
change is measured with great precision, using the hy- 
drogen masers at the DSN stations. The observable is 
the DSN frequency shift Q 

Aiy{t)=iyo--, (1) 
c at 

where £ is the overall optical distance (including diffrac- 
tion effects) traversed by a photon in both directions. 
[In the Pioneer Doppler experiments, the stability of 
the fractional drift at the S-band is on the order of 
Av/i'Q ~ 10~^^, for integration times on the order of 
10'^ s.] Doppler measurements provide the "range rate" 
of the spacecraft and therefore are affected by all the 
dynamical phenomena in the volume between the Earth 
and the spacecraft. 

Expanding upon what was discussed in Section Fill Al 
the received signal and the transmitter frequency (both 
are at S-band) as well as a 10 pulse per second timing 
reference from the FTS are fed to the Metric Data Assem- 
bly (MDA) . There the Doppler phase (difference between 
transmitted and received phases plus an added bias) is 
counted. That is, digital counters at the MDA record 
the zero crossings of the difference (i.e., Doppler, or al- 
ternatively the beat frequency of the received frequency 
and the exciter frequency). After counting, the bias is 
removed so that the true phase is produced. 

The system produces "continuous count Doppler" and 
it uses two counters. Every tenth of a second, a Doppler 
phase count is recorded from one of the counters. The 
other counter continues the counts. The recording alter- 
nates between the two counters to maintain a continuous 
unbroken count. The Doppler counts are at 1 MHz for 
S-band or 5 MHz for X-band. The wavelength of each S- 
band cycle is about 13 cm. Dividers or "time resolvers" 
further subdivide the cycle into 256 parts, so that frac- 
tional cycles are measured with a resolution of 0.5 mm. 
This accuracy can only be maintained if the Doppler is 
continuously counted (no breaks in the count) and coher- 
ent frequency standards are kept throughout the pass. It 
should be noted that no error is accumulated in the phase 
count as long as lock is not lost. The only errors are the 
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stability of the hydrogen maser and the resolution of the 
"resolver." 

Consequently, the JPL Doppler records are not fre- 
quency measurements. Rather, they are digitally counted 
measurements of the Doppler phase difference between 
the transmitted and received S-band frequencies, divided 
by the count time. 

Therefore, the Doppler observables, we will refer to, 
have units of cycles per second or Hz. Since total count 
phase observables are Doppler observables multiplied by 
the count interval Tc, they have units of cycles. The 
Doppler integration time refers to the total counting of 
the elapsed periods of the wave with the reference fre- 
quency of the hydrogen maser. The usual Doppler in- 
tegrating times for the Pioneer Doppler signals refers to 
the data sampled over intervals of 10 s, 60 s, 600 s, or 
1980 s. 



2. Range measurements 

A range measurement is made by phase modulating a 
signal onto the up-link carrier and having it echoed by the 
transponder. The transponder demodulates this ranging 
signal, filters it, and then re-modulates it back onto the 
down-link carrier. At the ground station, this returned 
ranging signal is demodulated and filtered. An instanta- 
neous comparison between the outbound ranging signal 
and the returning ranging signal that comes down would 
yield the two-way delay. Cross correlating the returned 
phase modulated signal with a ground duplicate yields 
the time delay. (See and references therein.) As the 
range code is repeated over and over, an ambiguity can 
exist. The orbit determination programs are then used 
to infer (some times with great difficulty) the number of 
range codes that exist between a particular transmitted 
code and its own corresponding received code. 

Thus, the ranging data are independent of the Doppler 
data, which represents a frequency shift of the radio car- 
rier wave without modulation. For example, solar plasma 
introduces a group delay in the ranging data but a phase 
advance in the Doppler data. 

Ranging data can also be used to distinguish an ac- 
tual range change from a fictitious range change seen in 
Doppler data that is caused by a frequency error [s^ . 
The Doppler frequency integrated over time (the accu- 
mulated phase) should equal the range change except for 
the difference introduced by charged particles 



3. Inferring position information from Doppler 

It is also possible to infer the position in the sky of 
a spacecraft from the Doppler data. This is accom- 
plished by examining the diurnal variation imparted to 
the Doppler shift by the Earth's rotation. As the ground 
station rotates underneath a spacecraft, the Doppler shift 



is modulated by a sinusoid. The sinusoid's amplitude de- 
pends on the declination angle of the spacecraft and its 
phase depends upon the right ascension. These angles 
can therefore be estimated from a record of the Doppler 
shift that is (at least) of several days duration. This al- 
lows for a determination of the distance to the spacecraft 
through the dynamics of spacecraft motion using stan- 
dard orbit theory contained in the orbit determination 
programs. 



C. Data preparation 

In an ideal system, all scheduled observations would 
be used in determining parameters of physical interest. 
However, there are inevitable problems that occur in data 
collection and processing that corrupt the data. So, at 
various stages of the signal processing one must remove 
or "edit" corrupted data. Thus, the need arises for ob- 
jective editing criteria. Procedures have been developed 
which attempt to excise corrupted data on the basis of 
objective criteria. There is always a temptation to elim- 
inate data that is not well explained by existing models, 
to thereby "improve" the agreement between theory and 
experiment. Such an approach may, of course, eliminate 
the very data that would indicate deficiencies in the a pri- 
ori model. This would preclude the discovery of improved 
models. 

In the processing stage that fits the Doppler samples, 
checks are made to ensure that there are no integer cycle 
slips in the data stream that would corrupt the phase. 
This is done by considering the difference of the phase 
observations taken at a high rate (10 times a second) 
to produce Doppler. Cycle slips often are dependent on 
tracking loop bandwidths, the signal to noise ratios, and 
predictions of frequencies. Blunders due to out-of-lock 
can be determined by looking at the original tracking 
data. In particular, cycle slips due to loss-of-lock stand 
out as a 1 Hz blunder point for each cycle slipped. 

If a blunder point is observed, the count is stopped and 
a Doppler point is generated by summing the preceding 
points. Otherwise the count is continued until a spec- 
ified maximum duration is reached. Cases where this 
procedure detected the need for cycle corrections were 
fiagged in the database and often individually examined 
by an analyst. Sometimes the data was corrected, but 
nominally the blunder point was just eliminated. This 
ensures that the data is consistent over a pass. However, 
it does not guarantee that the pass is good, because other 
errors can affect the whole pass and remain undetected 
until the orbit determination is done. 

To produce an input data file for an orbit determina- 
tion program, JPL has a software package known as the 
Radio Metric Data Selection, Translation, Revision, In- 
tercalation, Processing and Performance Evaluation Re- 
porting ( RMD-ST RIPPER) Program. As we discussed 
in Section lill B II this input file has data that can be in- 
tegrated over intervals with different durations: 10 s, 60 
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s, 600 s and 1980 s. This input Orbit Determination File 
(ODFILE) obtained from the RMDC group is the initial 
data set with which both the JPL and The Aerospace 
Corporation groups started their analyses. Therefore, 
the initial data file already contained some common data 
editing that the RMDC group had implemented through 
program flags, etc. The data set we started with had 
already been compressed to 60 s. So, perhaps there were 
some blunders that had already been removed using the 
initial STRIPPER program. 

The orbit analyst manually edits the remaining cor- 
rupted data points. Editing is done either by plotting 
the data residuals and deleting them from the fit or plot- 
ting weighted data residuals. That is, the residuals are 
divided by the standard deviation assigned to each data 
point and plotted. This gives the analyst a realistic view 
of the data noise during those times when the data was 
obtained while looking through the solar plasma. Apply- 
ing an 'W-(t" (cr is the standard deviation) test, where 
N is the choice of the analyst (usually 4-10) the analyst 
can delete those points that lie outside the N-a rejec- 
tion criterion without being biased in his selection. The 
N-a test, implemented in CHASMP, is very useful for 
data taken near solar conjunction since the solar plasma 
adds considerable noise to the data. This criterion later 
was changed to a similar criteria that rejects all data 
with residuals in the fit extending for more than ±0.025 
Hz from the mean. Contrariwise, the JPL analysis edits 
only very corrupted data; e.g., a blunder due to a phase 
lock loss, data with bad spin calibration, etc. Essentially 
the Aerospace procedure eliminates data in the tails of 
the Gaussian probability frequency distribution whereas 
the JPL procedure accepts this data. 

If needed or desired, the orbit analyst can choose to 
perform an additional data compression of the origi- 
nal navigation data. The JPL analysis does not apply 
any additional data compression and uses all the orig- 
inal data from the ODFILE as opposed to Aerospace's 
approach. Aerospace makes an additional compression 
of data within CHASMP. It uses the longest available 
data integration times which can be composed from ei- 
ther summing up adjacent data intervals or by using data 
spans with duration > 600 s. (Effectively Aerospace 
prefers 600 and 1980 second data intervals and applies 
a low-pass filter.) 

The total count of corrupted data points is about 10% 
of the total raw data points. The analysts' judgments 
play an important role here and is one of the main rea- 
sons that JPL and Aerospace have slightly different re- 
sults. (See Sections Ivbnd IVTI ) In Section Ivtve will show 
a typical plot (Figure |S1 below) with outliers present in 
the data. Many more outliers are off the plot. One would 
expect that the two different strategies of data compres- 
sion used by the two teams would result in significantly 
different numbers of total data points used in the two 
independent analyses. The infiuence of this fact on the 
solution estimation accuracy will be addressed in Section 
rvJ below. 



D. Data weighting 

Considerable effort has gone into accurately estimat- 
ing measurement errors in the observations. These errors 
provide the data weights necessary to accurately estimate 
the parameter adjustments and their associated uncer- 
tainties. To the extent that measurement errors are accu- 
rately modeled, the parameters extracted from the data 
will be unbiased and will have accurate sigmas assigned 
to them. Both JPL and Aerospace assign a standard un- 
certainty of 1 mm/s over a 60 second count time for the 
S-band Pioneer Doppler data. (Originally the JPL team 
was weighting the data by 2 mm/s uncertainty.) 

A change in the DSN antenna elevation angle also di- 
rectly affects the Doppler observables due to tropospheric 
refraction. Therefore, to correct for the infiuence of the 
Earth's troposphere the data can also be deweighted for 
low elevation angles. The phcnomcnological range cor- 
rection is given as 

a = ^no.inal(l + ^Y^^), (2) 

where dnominai is the basic standard deviation (in Hz) 
and 9e is the elevation angle in degrees Each leg 

is computed separately and summed. For Doppler the 
same procedure is used. First, Eq. (|2Jl is multiplied by 
60 s/Tc, where Tc is the count time. Then a numerical 
time differentiation of Eq. Q is performed. That is, 
Eq. (|2Jl is differenced and divided by the count time, Tc- 
(For more details on this standard technique see Refs. 

Ei-iii.) 

There is also the problem of data weighting for data 
influenced by the solar corona. This will be discussed in 
Section HVDl 

E. Spin calibration of the data 

The radio signals used by DSN to communicate with 
spacecraft are circularly polarized. When these signals 
are reflected from spinning spacecraft antennae a Doppler 
bias is introduced that is a function of the spacecraft spin 
rate. Each revolution of the spacecraft adds one cycle 
of phase to the up-link and the down-link. The up-link 
cycle is multiplied by the turn around ratio 240/221 so 
that the bias equals (1-1-240/221) cycles per revolution of 
the spacecraft. 

High-rate spin data is available for Pioneer 10 only 
up to July 17, 1990, when the DSN ceased doing spin 
calibrations. (See Section FlI B\ ) After this date, in or- 
der to reconstruct the spin behavior for the entire data 
span and thereby account for the spin bias in the Doppler 
signal, both analyses modeled the spin by performing in- 
terpolations between the data points. The JPL interpo- 
lation was non-linear with a high-order polynomial fit of 
the data. (The polynomial was from second up to sixth 
order, depending on the data quality.) The CHASMP 
interpolation was linear between the spin data points. 
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After a maneuver in mid-1993, there was not enough 
power left to support the IPP. But analysts still could get 
a rough determination approximately every six months 
using information obtained from the conscan maneuvers. 
No spin determinations were made after 1995. However, 
the archived conscan data could still yield spin data at 
every maneuver time if such work was approved. Further, 
as the phase center of the main antenna is slightly offset 
from the spin axis, a very small (but detectable) sine- 
wave signal appears in the high-rate Doppler data. In 
principle, this could be used to determine the spin rate for 
passes taken after 1993, but it has not been attempted. 
Also, the failure of one of the spin-down thrusters pre- 
vented precise spin calibration of the Pioneer 11 data. 

Because the spin rate of the Pioneers was changing over 
the data span, the calibrations also provide an indication 
of gas leaks that affect the acceleration of the spacecraft. 
A careful look at the records shows how this can be a 
problem. This will be discussed in Sections IVI Al and 
IVlllFI 



IV. BASIC THEORY OF SPACECRAFT 
NAVIGATION 

Accuracy of modern radio tracking techniques has pro- 
vided the means necessary to explore the gravitational 
environment in the solar system up to a limit never be- 
fore possible The major role in this quest belongs to 
relativistic celestial mechanics experiments with planets 
(e.g., passive radar ranging) and interplanetary space- 
craft (both Doppler and range experiments). Celestial 
mechanics experiments with spacecraft have been carried 
out by JPL since the early 1960's [Mli^l. The motiva- 
tion was to improve both the ephemerides of solar system 
bodies and also the knowledge of the solar system's dy- 
namical environment. This has become possible due to 
major improvements in the accuracy of spacecraft navi- 
gation, which is still a critical element for a number of 
space missions. The main objective of spacecraft navi- 
gation is to determine the present position and velocity 
of a spacecraft and to predict its future trajectory. This 
is usually done by measuring changes in the spacecraft's 
position and then, using those measurements, correcting 
(fitting and adjusting) the predicted spacecraft trajec- 
tory. 

In this section we will discuss the theoretical founda- 
tion that is used for the analysis of tracking data from 



interplanetary spacecraft. We describe the basic physical 
models used to determine a trajectory, given the data. 



A. Relativistic equations of motion 

The spacecraft ephemeris, generated by a numerical 
integration program, is a file of spacecraft positions and 
velocities as functions of ephemeris (or coordinate) time 
(ET) . The integrator requires the input of various param- 
eters. These include adopted constants (c, G, planetary 
mass ratios, etc.) and parameters that are estimated 
from fits to observational data (e.g., corrections to plan- 
etary orbital elements). 

The ephemeris programs use equations for point-mass 
relativistic gravitational accelerations. They are derived 
from the variation of a time-dependent, Lagrangian- 
action integral that is referenced to a non-rotating, solar- 
system, barycentric, coordinate frame. In addition to 
modeling point-mass interactions, the ephemeris pro- 
grams contain equations of motion that model terrestrial 
and lunar figure effects. Earth tides, and lunar phys- 
ical librations j^-ISOl. The programs treat the Sun, 
the Moon, and the nine planets as point masses in the 
isotropic, parameterized post-Newtonian, N-body metric 
with Newtonian gravitational perturbations from large, 
main-belt asteroids. 

Responding to the increasing demand of the naviga- 
tional accuracy, the gravitational field in the solar sys- 
tem is modeled to include a number of relativistic ef- 
fects that are predicted by the different metric theories 
of gravity. Thus, within the accuracy of modern exper- 
imental techniques, the parameterized post-Newtonian 
(PPN) approximation of modern theories of gravity pro- 
vides a useful starting point not only for testing these 
predictions, but also for describing the motion of self- 
gravitating bodies and test particles. As discussed in 
detail in |51|, the accuracy of the PPN limit (which is 
slow motion and weak field) is adequate for all foresee- 
able solar system tests of general relativity and a number 
of other metric theories of gravity. (For the most general 
formulation of the PPN formalism, see the works of Will 
and Nordtvedt [Hil.) 

For each body i (a planet or spacecraft anywhere in 
the solar syst em) , the point-mass acceleration is written 
as EHIilElillll 
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where /i^ is the "gravitational constant" of body i. It 
actually is its mass times the Newtonian constant: fii — 
Grrii. Also, ri(t) is the barycentric position of body i, 
Tij = |rj— r^l and Vi = |ri|. For planetary motion, each of 
these equations depends on the others. So they must be 
iterated in each step of the integration of the equations 
of motion. 

The barycentric acceleration of each body j due to 
Newtonian effects of the remaining bodies and the aster- 
oids is denoted by r^. In Eq. Q, (3 and 7 are the PPN 
parameters |U lE^- General relativity corresponds to 
/3 = 7 = 1, which we choose for our study. The Brans- 
Dicke theory is the most famous among the alternative 
theories of gravity. It contains, besides the metric tensor, 
a scalar field (p and an arbitrary coupling constant w, re- 
lated to the two PPN parameters as 7 = (3 — \. 
Equation Q allows the consideration of any problem in 
celestial mechanics within the PPN framework. 



B. Light time solution and time scales 

In addition to planetary equations of motion Eq. Ipj|l. 
one needs to solve the relativistic light propagation equa- 
tion in order to get the solution for the total light time 
travel. In the solar system, barycentric, space-time frame 
of reference this equation is given by: 



h-h = 
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where fj,Q is the gravitational constant of the Sun and /i^ 
is the gravitational constant of a planet, an outer plan- 
etary system, or the Moon, rf ^r^ andrf2 are the he- 
liocentric distances to the point of RF signal emission 
on Earth, to the point of signal reflection at the space- 
craft, and the relative distance between these two points. 
Correspondingly, r\,r\, and r\2 are similar distances rel- 
ative to a particular z-th body in the solar system. In 
the spacecraft light time solution, ti refers to the trans- 
mission time at a tracking station on Earth, and ^2 refers 
to the reflection time at the spacecraft or, for one-way 
|3fi| | data, the transmission time at the spacecraft. The 
reception time at the tracking station on Earth or at an 
Earth satellite is denoted by ^3. Hence, Eq. is the 
up-lcg light time equation. The corresponding down-leg 
light time equation is obtained by replacing subscripts as 
follows: 1 ^ 2 and 2 ^ 3. (See the details in |43|.) 

The spacecraft equations of motion relative to the so- 
lar system barycenter are essentially the same as given by 
Eq. The gravitational constants of the Sun, planets 
and the planetary systems are the values associated with 
the solar system barycentric frame of reference, which are 
obtained from the planetary ephemeris ^j. We treat a 
distant spacecraft as a point-mass particle. The space- 
craft acceleration is integrated numerically to produce 



the spacecraft ephemeris. The ephemeris is interpolated 
at the ephemeris time (ET) value of the interpolation 
epoch. This is the time coordinate t in Eqs. Q and 
i.e., t = ET. As such, ephemeris time means coor- 
dinate time in the chosen celestial reference frame. It is 
an independent variable for the motion of celestial bod- 
ies, spacecraft, and light rays. The scale of ET depends 
upon which reference frame is selected and one may use a 
number of time scales depending on the practical applica- 
tions. It is convenient to express ET in terms of Interna- 
tional Atomic Time (TAI). TAI is based upon the second 
in the International System of Units (Si). This second 
is defined to be the duration of 9,192,631,770 periods 
of the radiation corresponding to the transition between 
two hyperfine levels of the ground state of the cesium- 133 
atom [53. 

The differential equation relating ephemeris time (ET) 
in the solar system barycentric reference frame to TAI at 
a tracking station on Earth or on Earth satellite can be 
obtained directly from the Newtonian approximation to 
the N-body metric . This expression has the form 



dTAI 

dET ^ c 



-2{u-{U) + \v'-\{v^))+0{^\ (5) 

where U is the solar system gravitational potential eval- 
uated at the tracking station and v is the solar system 
barycentric velocity of the tracking station. The brack- 
ets ( ) on the right side of Eq. ||SJ) denote long-time 
average of the quantity contained within them. This av- 
eraging amounts to integrating out periodic variations in 
the gravitational potential, U , and the barycentric veloc- 
ity, , at the location of a tracking station. The desired 
time scale transformation is then obtained by using the 
planetary ephemeris to calculate the terms in Eq. (jSJ. 

The vector expression for the ephemeris/coordinate 
time (ET) in the solar system barycentric frame of ref- 
erence minus the TAI obtained from an atomic clock at 
a tracking station on Earth has the form |53| 

ET-TAI = 32.184 s + ^(i-Q •rO)-)-l(rr.r|)-f 



1 
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(rf.r 



c'^ifiQ + fiSa) 



ir1a-r^Sa) + i(j'i'-rf), (6) 



where and position and velocity vectors of point 
i relative to point j (they are functions of ET); super- 
script or subscript SSB denotes solar system barycenter; 
stands for the Sun; B for the Earth-Moon barycen- 
ter; J, Sa denote the Earth, Jupiter, and Saturn corre- 
spondingly, and A is for the location of the atomic clock 
on Earth which reads TAI. This approximated analytic 
result contains the clock synchronization term which de- 
pends upon the location of the atomic clock and five 
location-independent periodic terms. There are several 
alternate expressions that have up to several hundred 
additional periodic terms which provide greater accura- 
cies than Eq. (0 . The use of these extended expressions 



14 



provide transformations of ET - TAI to accuracies of 1 ns 

For the purposes of our study the Station Time (ST) is 
especially significant. This time is the atomic time TAI 
at a DSN tracking station on Earth, ST=TAIstation- This 
atomic time scale departs by a small amount from the 
"reference time scale." The reference time scale for a 
DSN tracking station on Earth is the Coordinated Uni- 
versal Time (UTC). This last is standard time for 0° lon- 
gitude. (For more details see [i^-lssj.! 

All the vectors in Eq. © except the geocentric position 
vector of the tracking station on Earth can be interpo- 
lated from the planetary ephemeris or computed from 
these quantities. Universal Time (UT) is the measure of 
time which is the basis for all civil time keeping. It is an 
observed time scale. The specific version used in JPL's 
Orbit Determination Program (ODP) is UTl. This is used 
to calculate mean sidereal time, which is the Greenwich 
hour angle of the mean equinox of date measured in the 
true equator of date. Observed UTl contains 41 short- 
term terms with periods between 5 and 35 days. They 
are caused by long-period solid Earth tides. When the 
sum of these terms, AUTl, is subtracted from UTl the 
result is called UTIR, where R means regularized. 

Time in any scale is represented as seconds past 1 
January 2000, 12^, in that time scale. This epoch is 
J2000.0, which is the start of the Juhan year 2000. The 
Julian Date for this epoch is JD 245,1545.0. Our analyses 
used the standard space-fixed J2000 coordinate system, 
which is provided by the International Celestial Reference 
Frame (ICRF). This is a quasi-inertial reference frame de- 
fined from the radio positions of 212 extragalactic sources 
distributed over the entire sky [s^. 

The variability of the earth-rotation vector relative to 
the body of the planet or in inertial space is caused by 
the gravitational torque exerted by the Moon, Sun and 
planets, displacements of matter in different parts of the 
planet and other excitation mechanisms. The observed 
oscillations can be interpreted in terms of mantle elas- 
ticity, earth flattening, structure and properties of the 
core-mantle boundary, rheology of the core, underground 
water, oceanic variability, and atmospheric variability on 
time scales of weather or climate. 

Several space geodesy techniques contribute to the con- 
tinuous monitoring of the Earth's rotation by the Inter- 
national Earth Rotation Service (lERS). Measurements 
of the Earth's rotation presented in the form of time de- 
velopments of the so-called Earth Orientation Param- 
eters (EOP). Universal time (UTl), polar motion, and 
the celestial motion of the pole (precession/nutation) 
are determined by Very Long-Baseline Interferometry 
(VLBI). Satellite geodesy techniques, such as satellite 
laser ranging (SLR) and using the Global Positioning 
System (GPS), determine polar motion and rapid varia- 
tions of universal time. The satellite geodesy programs 
used in the lERS allow determination of the time varia- 
tion of the Earth's gravity field. This variation reflects 
the evolutions of the Earth's shape and of the distribution 



of mass in the planet. The programs have also detected 
changes in the location of the center of mass of the Earth 
relative to the crust. It is possible to investigate other 
global phenomena such as the mass redistributions of the 
atmosphere, oceans, and solid Earth. 

Using the above experimental techniques, Universal 
time and polar motion are available daily with an accu- 
racy of about 50 picoseconds (ps). They are determined 
from VLBI astrometric observations with an accuracy of 
0.5 milliarcseconds (mas). Celestial pole motion is avail- 
able every flve to seven days at the same level of accuracy. 
These estimations of accuracy include both short term 
and long term noise. Sub-daily variations in Universal 
time and polar motion are also measured on a campaign 
basis. 

In summary, this dynamical model accounts for a num- 
ber of post-Newtonian perturbations in the motions of 
the planets, the Moon, and spacecraft. Light propaga- 
tion is correct to order c~^. The equations of motion of 
extended celestial bodies are valid to order c^"*. Indeed, 
this dynamical model has be en g ood enough to perform 
tests of general relativity |23, |5^ |52J ■ 



C. Standard modeling of small, non-gravitational 

forces 

In addition to the mutual gravitational interactions of 
the various bodies in the solar system and the gravita- 
tional forces acting on a spacecraft as a result of presence 
of those bodies, it is also important to consider a num- 
ber of non-gravitational forces which are important for 
the motion of a spacecraft. (Books and lengthy reports 
have been written about practically all of them. Consult 
Ref. [57I IHsjl for a general introduction.) 

The Jet Propulsion Laboratory's ODP accounts for 
many sources of non-gravitational accelerations. Among 
them, the most relevant to this study, are: i) solar radia- 
tion pressure, ii) solar wind pressure, iii) attitude-control 
maneuvers together with a model for unintentional space- 
craft mass expulsion due to gas leakage of the propulsion 
system. We can also account for possible influence of 
the interplanetary media and DSN antennae contribu- 
tions to the spacecraft radio tracking data and consider 
the torques produced by above mentioned forces. The 
Aerospace CHASMP code uses a model for gas leaks that 
can be adjusted to include the effects of the recoil force 
due to emitted radio power and anisotropic thermal ra- 
diation of the spacecraft. 

In principle, one could set up complicated engineering 
models to predict at least some of the effects. However, 
their residual uncertainties might be unacceptable for the 
experiment, in spite of the significant effort required. In 
fact, a constant acceleration produces a linear frequency 
drift that can be accounted for in the data analysis by a 
single unknown parameter. 

The figure against which we compare the effects of non- 
gravitational accelerations on the Pioneers' trajectories is 
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the expected error in the acceleration error estimations. 
This is on the order of 



(To 



2 X 10~^ cm/s^ 



(7) 



where ctq is a single determination accuracy related to ac- 
celeration measurements averaged over number of days. 
This would contribute to our result as (jn cto/VN- 
Thus, if no systematics are involved then will just 
tend to zero as time progresses. 

Therefore, the important thing is to know that these ef- 
fects (systematics) are not too large, thereby overwhelm- 
ing any possibly important signal (such as our anomalous 
acceleration). This will be demonstrated in Sections IVIII 
and lVlTTl 



D. Solar corona model and weighting 

The electron density and density gradient in the so- 
lar atmosphere influence the propagation of radio waves 
through the medium. So, both range and Doppler obser- 
vations at S-band are affected by the electron density in 
the interplanetary medium and outer solar corona. Since, 
throughout the experiment, the closest approach to the 
center of the Sun of a radio ray path was greater than 3.5 
i?Q, the medium may be regarded as coUisionless. The 
one way time delay associated with a plane wave passing 
through the solar corona is obtained 44, 46, 59] by inte- 
grating the group velocity of propagation along the ray's 
path, £: 
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where ne(t,r) is the free electron density in the solar 
plasma, c is the speed of light, and ncrit(i') is the critical 
plasma density for the radio carrier frequency v. The 
plus sign is applied for ranging data and the minus sign 
for Doppler data [60j . 

Therefore, in order to calibrate the plasma contribu- 
tion, one should know the electron density along the 
path. One usually decomposes the electron density, rig, 
into a static, steady-state part, ne(r), plus a fluctuation 
Sne{t,r), i.e., ne{t,r) = ne(r) + 6ne{t,r). The physical 
properties of the second term are hard to quantify. But 
luckily, its effect on Doppler observables and, therefore, 
on our results is small. (We will address this issue in Sec. 
IVIIBh On the contrary, the steady-state corona behav- 
ior is reasonably well known and several plasma models 
can be found in the literature j5gi]-j62i]. 

Consequently, while studying the effect of a system- 
atic error from propagation of the S-band carrier wave 
through the solar plasma, both analyses adopted the fol- 
lowing model for the electron density profile |44| : 



Rr. 



B 



(^) 



2.7 



r is the heliocentric distance to the immediate ray tra- 
jectory and is the helio-latitude normalized by the ref- 
erence latitude of 0o = 10°. The parameters r and 4> 
are determined from the trajectory coordinates on the 
tracking link being modeled. The parameters A, B, C 
are parameters chosen to describe the solar electron den- 
sity. (They are commonly given in two sets of units, 
meters or cm~'^ ,63].) They can be treated as stochas- 
tic parameters, to be determined from the fit. But in 
both analyses we ultimately chose to use the values de- 
termined from the recent solar corona studies done for 
the Cassini mission. These newly obtained values are: 
A = 6.0 X 10^, B = 2.0 X lO'*, C = 0.6 X 10^ aU in meters 
■ [This is what we will refer to as the "Cassini corona 
model."] 

Substitution of Eq. ^ into Eq. © results in the 
following steady-state solar corona contribution to the 
range model that we used in our analysis: 
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and ly are a reference frequency and the actual fre- 
quency of radio-wave [for Pioneer 10 analysis i^o = 2295 
MHz] , p is the impact parameter with respect to the Sun 
and F is a light-time correction factor. For distant space- 
craft this function is given as follows: 
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where and are the heliocentric radial distances to 
the target and to the Earth, respectively. Note that the 
sign of the solar corona range correction is negative for 
Doppler and positive for range. The Doppler correction 
is obtained from Eq. Hl()|l by simple time differentiation. 
Both analyses use the same physical model, Eq. (jlOfl . 
for the steady-state solar corona effect on the radio-wave 
propagation through the solar plasma. Although the ac- 
tual implementation of the model in the two codes is dif- 
ferent, this turns out not to be significant. (See Section 

mH) 

CHASMP can also consider the effect of temporal vari- 
ation in the solar corona by using the recorded history 
of solar activity. The change in solar activity is linked to 
the variation of the total number of sun spots per year 
as observed at a particular wavelength of the solar radia- 
tion, A=10.7 cm. The actual data corresponding to this 
variation is given in Ref. [b^. CHASMP averages this 
data over 81 days and normalizes the value of the flux by 
150. Then it is used as a time-varying scaling factor in 
Eq. ((TUJ). The result is referred to as the "F10.7 model." 

Next we come to corona data weighting. JPL's ODP 
does not apply corona weighting. On the other hand. 
Aerospace's CHASMP can apply corona weighting if de- 
sired. Aerospace uses a standard weight augmented by 
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a weight function that accounts for noise introduced by 
solar plasma and low elevation. The weight values are ad- 
justed so that i) the post-fit weighted sum of the squares 
is close to unity and ii) approximately uniform noise in 
the residuals is observed throughout the fit span. 

Thus, the corresponding solar-corona weight function 

is: 

where, for range data, k is an input constant nominally 
equal to 0.005 light seconds, and v are a reference 
frequency and the actual frequency, p is the trajectory's 
impact parameter with respect to the Sun in km, and Rq 
is the solar radius in km |66j| . The solar-corona weight 
function for Doppler is essentially the same, but obtained 
by numerical time differentiation of Eq. (|12|l . 



by two to three days of gas leakage, large enough to be 
observable in the Doppler data |67| . 

Typically the Doppler data is time averaged over 10 
to 33 minutes, which significantly reduces the high- 
frequency Doppler noise. The residuals represent our fit. 
They are converted from units of Hz to Doppler velocity 
by the formula js^ 

r A 1 r^'^ DSN n\ 
Au DSN = 7: , 13) 

2 l/Q 

where vq is the downlink carrier frequency, ~ 2.29 GHz, 
Av is the Doppler residual in Hz from the fit, and c is 
the speed of light. 

As an illustration, consider the fit to one of the Pio- 
neer 10 maneuvers, # 17, on 22 December 1993, given in 
Figure[Sl This was particularly well covered by low-noise 



E. Modeling of maneuvers 

There were 28 Pioneer 10 maneuvers during our data 
interval from 3 January 1987 to 22 July 1998. Imperfect 
coupling of the hydrazine thrusters used for the spin ori- 
entation maneuvers produced integrated velocity changes 
of a few millimeters per second. The times and durations 
of each maneuver were provided by NASA/Ames. JPL 
used this data as input to ODP. The Aerospace team 
used a slightly different approach. In addition to the orig- 
inal data, CHASMP used the spin-rate data file to help 
determine the times and duration of maneuvers. The 
CHASMP determination mainly agreed with the data 
used by JPL. [There were minor variations in some of 
the times, one maneuver was split in two, and one extra- 
neous maneuver was added to Interval II to account for 
data not analyzed (see below).] 

Because the effect on the spacecraft acceleration could 
not be determined well enough from the engineering 
telemetry, JPL included a single unknown parameter in 
the fitting model for each maneuver. In JPL's ODP anal- 
ysis, the maneuvers were modeled by instantaneous ve- 
locity increments at the beginning time of each maneu- 
ver (instantaneous burn model). [Analyses of individual 
maneuver fits show the residuals to be small.] In the 
CHASMP analysis, a constant acceleration acting over 
the duration of the maneuver was included as a param- 
eter in the fitting model (finite burn model). Analyses 
of individual maneuver fits show the residuals are small. 
Because of the Pioneer spin, these accelerations are im- 
portant only along the Earth-spacecraft line, with the 
other two components averaging out over about 50 revo- 
lutions of the spacecraft over a typical maneuver duration 
of 10 minutes. 

By the time Pioneer 11 reached Saturn, the pattern 
of the thruster firings was understood. Each maneuver 
caused a change in spacecraft spin and a velocity incre- 
ment in the spacecraft trajectory, immediately followed 
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FIG. 5: The Doppler residuals after a fit for maneuver #17 
on 23 December 1993. 

Doppler data near solar opposition. Before the start of 
the maneuver, there is a systematic trend in the residu- 
als which is represented by a cubic polynomial in time. 
The standard error in the residuals is 0.095 mm/s. After 
the maneuver, there is a relatively small velocity discon- 
tinuity of —0.90 ± 0.07 mm/s. The discontinuity arises 
because the model fits the entire data interval. In fact, 
the residuals increase after the maneuver. By 11 January 
1994, 19 days after the maneuver, the residuals are scat- 
tered about their pre-maneuver mean of —0.15 mm/s. 

For purposes of characterizing the gas leak immedi- 
ately after the maneuver, we fit the post-maneuver resid- 
uals by a two-parameter exponential curve. 
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0.15 mm/s. 



(14) 



The best fit yields vq — 0.808 mm/s and the time con- 
stant T is 13.3 days, a reasonable result. The time deriva- 
tive of the exponential curve yields a residual acceleration 
immediately after the maneuver of 7.03 x 10~* cm/s^. 
This is close to the magnitude of the anomalous acceler- 
ation inferred from the Doppler data, but in the opposite 
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direction. However the gas leak rapidly decays and be- 
comes negligible after 20 days or so. 



F. Orbit determination procedure 

Our orbit determination procedure first determines the 
spacecraft's initial position and velocity in a data inter- 
val. For each data interval, we then estimate the mag- 
nitudes of the orientation maneuvers, if any. The anal- 
yses are modeled to include the effects of planetary per- 
turbations, radiation pressure, the interplanetary media, 
general relativity, and bias and drift in the Doppler and 
range (if available). Planetary coordinates and solar sys- 
tem masses are obtained using JPL's Export Planetary 
Ephemeris DE405, where DE stands for the Development 
Ephemeris. [Earlier in the study, DE200 was used. See 
Section ESI] 

We include models of precession, nutation, sidereal ro- 
tation, polar motion, tidal effects, and tectonic plates 
drift. Model values of the tidal deceleration, nonunifor- 
mity of rotation, polar motion, Love numbers, and Chan- 
dler wobble are obtained observationally, by means of 
Lunar and Satellite Laser Ranging (LLR and SLR) tech- 
niques and VLBI. Previously they were combined into 
a common publication by either the International Earth 
Rotation Service (lERS) or by the United States Naval 
Observatory (USNO). Currently this information is pro- 
vided by the ICRF. JPL's Earth Orientation Parameters 
(EOP) is a major source contributor to the ICRF. 

The implementation of the J2000.0 reference coordi- 
nate system in CHASMP involves only rotation from the 
Earth-fixed to the J2000.0 reference frame and the use 
of JPL's DE200 planetary ephemeris [Hi|. The rota- 
tion from J2000.0 to Earth-fixed is computed from 
a series of rotations which include precession, nutation, 
the Greenwich hour angle, and pole wander. Each of 
these general categories is also a multiple rotation and 
is treated separately by most software. Each separate 
rotation matrix is chain multiplied to produce the final 
rotation matrix. 

CHASMP, however, does not separate precession and 
nutation. Rather, it combines them into a single matrix 
operation. This is achieved by using a different set of an- 
gles to describe precession than is used in the ODP. (See 
a description of the standard set of angles in [i^.) These 
angles separate luni-solar precession from planetary pre- 
cession. Luni-solar precession, being the linear term of 
the nutation series for the nutation in longitude, is com- 
bined with the nutation in longitude from the DE200 
ephemeris tape |70| . 

Both JPL's ODP and The Aerospace Corporation's 
CHASMP use the JPL/Earth Orientation Parameters 
(EOP) values. This could be a source of common er- 
ror. However the comparisons between EOP and lERS 
show an insignificant difference. Also, only secular terms, 
such as precession, can contribute errors to the anoma- 
lous acceleration. Errors in short period terms are not 



correlated with the anomalous acceleration. 



G. Parameter estimation strategies 

During the last few decades, the algorithms of orbital 
analysis have been extended to incorporate Kalman-filter 
estimation procedure that is based on the concept of 
"process noise" (i.e., random, non-systematic forces, or 
random- walk effects) . This was motivated by the need to 
respond to the significant improvement in observational 
accuracy and, therefore, to the increasing sensitivity to 
numerous small perturbing factors of a stochastic nature 
that are responsible for observational noise. This ap- 
proach is well justified when one needs to make accurate 
predictions of the spacecraft's future behavior using only 
the spacecraft's past hardware and electronics state his- 
tory as well as the dynamic environment conditions in 
the distant craft's vicinity. Modern navigational software 
often uses Kalman filter estimation since it more easily 
allows determination of the temporal noise history than 
does the weighted least-squares estimation. 

To take advantage of this while obtaining JPL's orig- 
inal results 0, 01 discussed in Section JPL used 
batch-sequential methods with variable batch sizes and 
process noise characteristics. That is, a batch-sequential 
filtering and smoothing algorithm with process noise was 
used with ODP. In this approach any small anomalous 
forces may be treated as stochastic parameters affecting 
the spacecraft trajectory. As such, these parameters are 
also responsible for the stochastic noise in the observa- 
tional data. To better characterize these noise sources, 
one splits the data interval into a number of constant or 
variable size batches and makes assumptions on possible 
statistical properties of these noise factors. One then esti- 
mates the mean values of the unknown parameters within 
the batch and also their second statistical moments. 

Using batches has the advantage of dealing with a 
smaller number of experimental data segments. We ex- 
perimented with a number of different constant batch 
sizes; namely, 0, 5, 30, and 200 day batch sizes. (Later 
we also used 1 and 10 day batch sizes.) In each batch 
one estimates the same number of desired parameters. 
So, one expects that the smaller the batch size the larger 
the resulting statistical errors. This is because a smaller 
number of data points is used to estimate the same num- 
ber of parameters. Using the entire data interval as a 
single batch while changing the process noise a priori val- 
ues is expected in principle (see below) to yield a result 
identical to the least-squares estimation. In the single 
batch case, it would produce only one solution for the 
anomalous acceleration. 

There is another important parameter that was taken 
into account in the statistical data analysis reported here. 
This is the expected correlation time for the underly- 
ing stochastic processes (as well as the process noise) 
that may be responsible for the anomalous acceleration. 
For example, using a zero correlation time is useful in 
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searches for an ap that is generated by a random pro- 
cess. One therefore expects that an ap estimated from 
one batch is statistically independent (uncorrelated) from 
those estimated from other batches. Also, the use of fi- 
nite correlation times indicates one is considering an ap 
that may show a temporal variation within the data in- 
terval. We experimented with a number of possible cor- 
relation times and will discuss the corresponding assump- 
tions when needed. 

In each batch one estimates solutions for the set of de- 
sired parameters at a specified epoch within the batch. 
One usually chooses to report solutions corresponding to 
the beginning, middle, or end of the batch. General co- 
ordinate and time transformations (discussed in Section 



IIVB|I are then used to report the solution in the epoch 
chosen for the entire data interval. One may also adjust 
the solutions among adjacent batches by accounting for 
possible correlations. This process produces a smoothed 
solution for the set of solved-for parameters. More de- 
tails on this so called "batch-sequential algorithm with 
smoothing filter" are available in Refs. |4H-[43*|. 

Even without process noise, the inversion algorithms of 
the Kalman formulation and the weighted least-squares 
method seem radically different. But as shown in [7ll |. 
if one uses a single batch for all the data and if one 
uses certain assumptions about, for instance, the pro- 
cess noise and the smoothing algorithms, then the two 
methods are mathematically identical. When introduc- 
ing process noise, an additional process noise matrix is 
also added into the solution algorithm. The elements 
of this matrix are chosen by the user as prescribed by 
standard statistical techniques used for navigational data 
processing. 

For the recent results reported in Section lVll JPL used 
both the batch-sequential and the weighted least-squares 
estimation approaches. JPL originally implemented only 
the batch-sequential method, which yielded the detec- 
tion (at a level smaller than could be detected with any 
other spacecraft) of an annual oscillatory term smaller in 
size than the anomalous acceleration [13. (This term is 
discussed in Section HX CI ) The recent studies included 
weighted least-squares estimation to see if this annual 
term was a calculational anomaly. 

The Aerospace Corporation uses only the weighted 
least-squares approach with its CHASMP software. A 
test is used as an indicator of the quality of the fit. In this 
case, the anomalous acceleration is treated as a constant 
parameter over the entire data interval. To solve for ap 
one estimates the statistical weights for the data points 
and then uses these in a general weighted least-squares 
fashion. Note that the weighted least-squares method can 
obtain a result similar to that from a batch-sequential ap- 
proach (with smoothing filter, zero correlation time and 
without process noise) by cutting the data interval into 
smaller pieces and then looking at the temporal variation 
among the individual solutions. 

As one will see in the following, in the end, both pro- 
grams yielded very similar results. The differences be- 



tween them can be mainly attributed to (other) system- 
atics. This gives us confidence that both programs and 
their implemented estimation algorithms are correct to 
the accuracy of this investigation. 



V. ORIGINAL DETECTION OF THE 
ANOMALOUS ACCELERATION 

A. Early JPL studies of the anomalous Pioneer 
Doppler residuals 

As mentioned in the introduction, by 1980 Pioneer 10 
was at 20 AU, so the solar radiation pressure accelera- 
tion had decreased to < 5 x 10~* cm/s^. Therefore, a 
search for unmodeled accelerations (at first with the fur- 
ther out Pioneer 10) could begin at this level. With the 
acceptance of a proposal of two of us (JDA and ELL) 
to participate in the Heliospheric Mission on Pioneer 10 
and 11, such a search began in earnest [tJ. 

The JPL analysis of unmodeled accelerations used the 
JPL's Orbit Determination Program (ODP) 41]- 4J. 
Over the years the data continually indicated that the 
largest systematic error in the acceleration residuals is a 
constant bias of ap ^ (8 ± 3) x 10^* cm/s^, directed to- 
ward the Sun (to within the beam-width of the Pioneers' 
antennae [T^I. 

As stated previously, the analyses were modeled to 
include the effects of planetary perturbations, radiation 
pressure, the interplanetary media, general relativity, to- 
gether with bias and drift in the Doppler signal. Plane- 
tary coordinates and the solar system masses were taken 
from JPL's Export Planetary Ephcmeris DE405, refer- 
enced to ICRF. The analyses used the standard space- 
fixed J2000 coordinate system with its associated JPL 
planetary ephemeris DE405 (or earlier, DE200). The 
time-varying Earth orientation in J2000 coordinates is 
defined by a 1998 version of JPL's EOP file, which ac- 
counts for the inertial precession and nutation of the 
Earth's spin axis, the geophysical motion of the Earth's 
pole with respect to its spin axis, and the Earth's time 
varying spin rate. The three-dimensional locations of 
the tracking stations in the Earth's body-fixed coordi- 
nate system (geocentric radius, latitude, longitude) were 
taken from a set recommended by ICRF for JPL's DE405. 

Consider I'obs , the frequency of the re-transmitted sig- 
nal observed by a DSN antennae, and I'modeij the pre- 
dicted frequency of that signal. The observed, two-way 
anomalous effect can be expressed to first order in v/c as 

M 



[f^obs(^) — t'model(^)]BSM ^ 



2ap t 



t'model — t^O 



1 - 



2l^model {t) 



(15) 



Here, vq is the reference frequency, the factor 2 is because 
we use two- and three-way data |3g. Vnodei is the mod- 
eled velocity of the spacecraft due to the gravitational 
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and other large forces discussed in Section Hvl (This ve- 
locity is outwards and hence produces a red shift.) We 
have already included the sign showing that ap is inward. 
(Therefore, ap produces a slight blue shift on top of the 
larger red shift.) By DSN convention 38], the first of Eqs. 

(^nj is [Al^obs — Al'modellusual = — [Al^obs — Aj/modellDSN • 

Over the years the anomaly remained in the data of 
both Pioneer 10 and Pioneer 11 |7^. (See FigureEl) 
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FIG. 6: ODP plots, as a function of distance from the Sun, 
of accelerations on Pioneers 10/11. The accelerations are a) 
the calculated solar radiation acceleration (top line), b) the 
unmodeled acceleration (bottom line), and c) the sum of the 
two above (middle line) \7at . 

In order to model any unknown forces acting on Pi- 
oneer 10, the JPL group introduced a stochastic accel- 
eration, exponentially correlated in time, with a time 
constant that can be varied. This stochastic variable is 
sampled in ten-day batches of data. We found that a 
correlation time of one year produces good results. We 
did, however, experiment with other time constants as 
well, including a zero correlation time (white noise). The 
result of applying this technique to 6.5 years of Pioneer 
10 and 11 data is shown in Figure [7| The plotted points 
represent our determination of the stochastic variable at 
ten-day sample intervals. We plot the stochastic variable 
as a function of heliocentric distance, not time, because 
that is more fundamental in searches for trans-Neptunian 
sources of gravitation. 

As possible "perturbative forces" to explain this bias, 
we considered gravity from the Kuiper belt, gravity from 
the galaxy, spacecraft "gas leaks," errors in the plane- 
tary ephemeris, and errors in the accepted values of the 
Earth's orientation, precession, and nutation. We found 
that none of these mechanisms could explain the appar- 
ent acceleration, and some were three orders of magni- 
tude or more too small. [We also ruled out a number 
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FIG. 7: An ODP plot of the early unmodeled accelerations 
of Pioneer 10 and Pioneer 11, from about 1981 to 1989 and 
1977 to 1989, respectively [t^. 

of specific mechanisms involving heat radiation or "gas 
leaks," even though we feel these are candidates for the 
cause of the anomaly. We will return to this in Sections 

rvnland lvml ] 

We concluded from the JPL-ODP analysis, that 
there is an unmodeled acceleration, ap, towards the Sun 
of (8.09 ± 0.20) X 10"*^ cm/s^ for Pioneer 10 and of 
(8.56 ±0.15) X 10~^ cm/s^ for Pioneer 11. The error was 
determined by use of a five-day batch sequential filter 
with radial acceleration as a stochastic parameter sub- 
ject to white Gaussian noise (~ 500 independent five-day 
samples of radial acceleration) [t^ . No magnitude varia- 
tion of ap with distance was found, within a sensitivity of 
do = 2 X 10^^ cm/s^ over a range of 40 to 60 AU. All our 
errors are taken from the covariance matrices associated 
with the least-squares data analysis. The assumed data 
errors are larger than the standard error on the post-fit 
residuals. [For example, the Pioneer S-band Doppler er- 
ror was set at 1 mm/s at a Doppler integration time of 60 
s, as opposed to a characteristic value of 0.3 mm/s.] 
Consequently, the quoted errors are realistic, not formal, 
and represent our attempt to include systematics and a 
reddening of the noise spectrum by solar plasma. Any 
spectral peaks in the post-fit Pioneer Doppler residuals 
were not significant at a 90% confidence level ^3 • 



B. First Aerospace study of the apparent Pioneer 
acceleration 

With no explanation of this data in hand, our atten- 
tion focused on the possibility that there was some error 
in JPL's ODP. To investigate this, an analysis of the raw 
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data was performed using an independent program, The 
Aerospace Corporation's Compact High Accuracy Satel- 
lite Motion Program (CHASMP) ^ - one of the stan- 
dard Aerospace orbit analysis programs. CHASMP 's or- 
bit determination module is a development of a program 
called POEAS (Planetary Orbiter Error Analysis Study 
program) that was developed at JPL in the early 1970's 
independently of JPL's ODP. As far as we know, not a 
single line of code is common to the two programs [t^ ■ 

Although, by necessity, both ODP and CHASMP use 
the same physical principles, planetary ephemeris, and 
timing and polar motion inputs, the algorithms are oth- 
erwise quite different. If there were an error in either 
program, they would not agree. 

Aerospace analyzed a Pioneer 10 data arc that was 
initialized on 1 January 1987 at 16 hr (the data itself 
started on 3 January) and ended at 14 December 1994, 
hr. The raw data set was averaged to 7560 data points of 
which 6534 points were used. This CHASMP analysis of 
Pioneer 10 data also showed an unmodeled acceleration 
in a direction along the radial toward the Sun [t^ . The 
value is (8.65 ± 0.03) x 10~^ cm/s^, agreeing with JPL's 
result. The smaller error here is because the CHASMP 
analysis used a batch least-squares fit over the whole orbit 
[To. 123 , not looking for a variation of the magnitude of 
ap with distance. 

Without using the apparent acceleration, CHASMP 
shows a steady frequency drift ^38] of about —6 x 10^^ 
Hz/s, or 1.5 Hz over 8 years (one-way only). (See Fig- 
ure IHl) This equates to a clock acceleration, —at, of 
—2.8 X 10"^* s/s^. The identity with the apparent Pio- 
neer acceleration is 



at = ap/c. 



(16) 



The drift in the Doppler residuals (observed minus com- 
puted data) is seen in Figure IHl 

The drift is clear, definite, and cannot be removed 
without either the added acceleration, ap, or the inclu- 
sion in the data itself of a frequency drift, i.e., a "clock 
acceleration" at. If there were a systematic drift in the 
atomic clocks of the DSN or in the time-reference stan- 
dard signals, this would appear like a non- uniformity of 
time; i.e., all clocks would be changing with a constant 
acceleration. We now have been able to rule out this 
possibility. (See Section IXIDI l 

Continuing our search for an explanation, we consid- 
ered the possibilities: i) that the Pioneer 10/11 spacecraft 
had internal systematic properties, undiscovered because 
they are of identical design, and ii) that the accelera- 
tion was due to some not-understood viscous drag force 
(proportional to the approximately constant velocity of 
the Pioneers). Both these possibilities could be investi- 
gated by studying spin-stabilized spacecraft whose spin 
axes are not directed towards the Sun, and whose orbital 
velocity vectors are far from being radially directed. 

Two candidates were Galileo in its Earth-Jupiter mis- 
sion phase and Ulysses in Jupiter-perihelion cruise out 
of the plane of the ecliptic. As well as Doppler, these 
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FIG. 8: CHASMP two-way Doppler residuals (observed 
Doppler velocity minus model Doppler velocity) for Pioneer 
10 vs. time. 1 Hz is equal to 65 mm/s range change per 
second. The model is fuUy-relativistic. The solar system's 
gravitational field is represented by the Sun and its planetary 
systems jS^. 
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FIG. 9: CHASMP best fit for the Pioneer 10 Doppler residuals 
with the anomalous acceleration taken out. After adding one 
more parameter to the model (a constant radial acceleration) 
the residuals are distributed about zero Doppler velocity with 
a systematic variation ~ 3.0 mm/s on a time scale of ~ 3 
months. The outliers on the plot were rejected from the fit. 
[The quality of the fit may be determined by the ratio of 
residuals to the downlink carrier frequency, i/q ~ 2.29 GHz.] 



spacecraft also yielded a considerable quantity of range 
data. By having range data one can tell if a spacecraft 
is accumulating a range effect due to a spacecraft accel- 
eration or if the orbit determination process is fooled by 
a Doppler frequency rate bias. 
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C. Galileo measurement analysis 

We considered the dynamical behavior of GaHleo's tra- 
jectory during its cruise flight from second Earth en- 
counter (on 8 December 1992) to arrival at Jupiter. [This 
period ends just before the Galileo probe release on 13 
July 1995. The probe reached Jupiter on 7 December 
1995.] During this time the spacecraft traversed a dis- 
tance of about 5 AU with an approximately constant ve- 
locity of 7.19(4) km/s. 

A quick JPL look at limited Galileo data (241 days 
from 8 January 1994 to 6 September 1994) demonstrated 
that it was impossible to separate solar radiation effects 
from an anomalous constant acceleration. The Sun was 
simply too close and the radiation cross-section too large. 
The nominal value obtained was ^ 8 x 10~^ cm/s^. 

The Aerospace's analysis of the Galileo data covered 
the same arc as JPL and a second arc from 2 December 
1992 to 24 March 1993. The analysis of Doppler data 
from the first arc resulted in a determination for ap of 
^ (8 ± 3) X 10~* cm/s^, a value similar to that from 
Pioneer 10. But the correlation with solar pressure was 
so high (0.99) that it is impossible to decide whether solar 
pressure is a contributing factor l80l|. 

The second data arc was 113 days long, starting six 
days prior to the second Earth encounter. This solution 
was also too highly correlated with solar pressure, and 
the data analysis was complicated by many mid-course 
maneuvers in the orbit. The uncertainties in the maneu- 
vers were so great, a standard null result could not be 
ruled out. 

However, there was an additional result from the data 
of this second arc. This arc was chosen for study because 
of the availability of ranging data. It had 11596 Doppler 
points of which 10111 were used and 5643 range points 
of which 4863 used. The two-way range change and time 
integrated Doppler are consistent (see Figure [TU)! to ~ 
4 m over a time interval of one day. For comparison, 
note that for a time of t ^ 1 day, (apt^/2) ~ 3 m. For 
the apparent acceleration to be the result of hardware 
problems at the tracking stations, one would need a linear 
frequency drift at all the DSN stations, a drift that is not 
observed. 



D. Ulysses measurement analysis 

1. J PL's analysis 

An analysis of the radiation pressure on Ulysses, in 
its out-of-the-ecliptic journey from 5.4 AU near Jupiter 
in February 1992 to the perihelion at 1.3 AU in Febru- 
ary 1995, found a varying profile with distance jsJl- The 
orbit solution requires a periodic updating of the solar 
radiation pressure. The radio Doppler and ranging data 
can be fit to the noise level with a time- varying solar con- 
stant in the fitting model [s^l- We obtained values for 
the time-varying solar constant determined by Ulysses 
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FIG. 10: Galileo best fit Doppler and range residuals using 
CHASMP. 



navigational data during this south polar pass The 
inferred solar constant is about 40 percent larger at per- 
ihelion (1.3 AU) than at Jupiter (5.2 AU), a physical 
impossibility! 

We sought an alternative explanation. Using physical 
parameters of the Ulysses spacecraft, we first converted 
the time- varying values of the solar constant to a posi- 
tive (i.e., outward) radial spacecraft acceleration, (x^ , as 
a function of heliocentric radius. Then we fit the values 
of ar with the following model: 



/C/o A cos 6>(r) 
cM 



(17) 



where r is the heliocentric distance in AU, M is the to- 
tal mass of the spacecraft, Jq = 1367 W/m^(AU)'^ is 
the (effective-temperature Stefan-Boltzmann) "solar ra- 
diation constant" at 1 AU, A is the cross-sectional area 
of the spacecraft and 9{r) is the angle between the di- 
rection to the Sun at distance r and orientation of the 
antennae. [For the period analyzed 9{r) was almost a 
constant. Therefore its average value was used which 
corresponded to (cos 6'(r)) « 0.82.] Optical parameters 
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defining the reflectivity and emissivity of the spacecraft's 
surface were taken to yield JC ~ 1.8. (See Section rVll Al 
for a discussion on solar radiation pressure.) Finally, the 
parameter ap(jj^ was determined by linear least squares. 
The best-fit value was obtained 

ap(c7) = (12±3) X 10"^ cm/s^ (18) 

where both random and systematic errors are included. 

So, by interpreting this time variation as a true so- 
lar pressure plus a constant radial acceleration, we found 
that Ulysses was subjected to an unmodeled acceleration 
towards the Sun of (12 ± 3) xlO"** cm/s^. 

Note, however, that the determined constant ap([/) is 
highly correlated with solar radiation pressure (0.888). 
This shows that the constant acceleration and the solar- 
radiation acceleration are not independently determined, 
even over a heliocentric distance variation from 5.4 to 1.3 
AU. 



2. Aerospace's analysis 

The next step was to perform a detailed calculation of 
the Ulysses orbit from near Jupiter encounter to Sun per- 
ihelion, using CHASMP to evaluate Doppler and ranging 
data. The data from 30 March 1992 to 11 August 1994 
was processed. It consisted of 50213 Doppler points of 
which 46514 were used and 9851 range points of which 
8465 were used. 

Such a calculation would in principle allow a more pre- 
cise and believable differentiation between an anomalous 
constant acceleration towards the Sun and systematics. 
Solar radiation pressure and radiant heat systematics are 
both larger on Ulysses than on the Pioneers. 

However, this calculation turned out to be a much 
more difficult than imagined. Because of a failed nutation 
damper, an inordinate number of spacecraft maneuvers 
were required (257). Even so, the analysis was completed. 
But even though the Doppler and range residuals were 
consistent as for Galileo, the results were disheartening. 
For an unexpected reason, any fit is not significant. The 
anomaly is dominated by (what appear to be) gas leaks 
|83]. That is, after each maneuver the measured anomaly 
changes. The measured anomalies randomly change sign 
and magnitude. The values go up to about an order of 
magnitude larger than ap. So, although the Ulysses data 
was useful for range/Doppler checks to test models (see 
Section IXI D|) , like Galileo it could not provide a good 
number to compare to ap. 



VI. RECENT RESULTS 

Recent changes to our strategies and orbit determi- 
nation programs, leading to new results, are threefold. 
First, we have added a longer data arc for Pioneer 10, 
extending the data studied up to July 1998. The entire 



data set used (3 Jan. 1987 to 22 July 1998) covers a he- 
liocentric distance interval from 40 AU to 70.5 AU [83| . 
[Pioneer 11 was much closer in (22.42 to 31.7 AU) than 
Pioneer 10 during its data interval (5 January 1987 to 
1 October 1990).] For later use in discussing systemat- 
ics, we here note that in the ODP calculations, masses 
used for the Pioneers were Mpio lo = 251.883 kg and 
Afp„ n = 239.73 kg. CHASMP used 251.883 kg for 
both [13. As the majority of our results are from Pio- 
neer 10, we will make Mq — 251.883 kg to be our nominal 
working mass. 

Second, and as we discuss in the next subsection, we 
have studied the spin histories of the craft. In particular, 
the Pioneer 10 history exhibited a very large anomaly in 
the period 1990.5 to 1992.5. This led us to take a closer 
look at any possible variation of ap among the three time 
intervals: The JPL analysis defined the intervals as I (3 
Jan. 1987 to 17 July 1990); II (17 July 1990 to 12 July 
1992) bounded by 49.5 to 54.8 AU; and III (12 July 1992 
to 22 July 1998). (CHASMP used shghtly different in- 
tervals j85|) The total updated data set now consists of 
20,055 data points for Pioneer 10. (10,616 data points 
were used for Pioneer 11.) This helped us to better un- 
derstand the systematic due to gas leaks, which is taken 
up in Section lynm 

Third, in looking at the detailed measurements of ap 
as a function of time using ODP, we found an anomalous 
oscillatory annual term, smaller in size than the anoma- 
lous acceleration As mentioned in Scction llV Gl and 
as will be discussed in detail in Section llX CI we wanted 
to make sure this annual term was not an artifact of our 
computational method. For the latest results, JPL used 
both the batch-sequential and the least-squares methods. 

All our recent results obtained with both the JPL and 
The Aerospace Corporation software have given us a bet- 
ter understanding of systematic error sources. At the 
same time they have increased our confidence in the de- 
termination of the anomalous acceleration. We present 
a description and summary of the new results in the rest 
of this section. 



A. Analysis of the Pioneer spin history 

Both Pioneers 10 and 11 were spinning down during 
the respective data intervals that determined their ap 
values. Because any changes in spacecraft spin must be 
associated with spacecraft torques (which for lack of a 
plausible external mechanism we assume are internally 
generated), there is also a possibility of a related inter- 
nally generated translational force along the spin axis. 
Therefore, it is important to understand the effects of 
the spin anomalies on the anomalous acceleration. In 
Figures ITTI and IT^ we show the spin histories of the two 
craft during the periods of analysis. 

Consider Pioneer 10 in detail. In time Interval I 
there is a slow spin down at an average rate (slope) of 

(-0.0181±0.0001) rpm/yr. Indeed, a closer look at the 
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FIG. 11: The spin history of Pioneer 10. The vertical lines in- 
dicate the times when precession maneuvers were made. How 
this spin data was obtained is described in Section HllEI The 
final data points were obtained at the times of maneuvers, the 
last being in 1995. 
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FIG. 12: The spin history of Pioneer 11 over the period of 
analysis. The vertical lines indicate the times when precession 
maneuvers were made. This spin calibration was done by 
the DSN until 17 July 1990. At that time the DSN ceased 
doing spin calibrations. From 1990 until the loss of coherent 
Doppler, orbit analysts made estimates of the spin rate. 

curve (either by eye or from an expanded graph) shows 
that the spin down is actually slowing with time (the 
curve is flattening) . This last feature will be discussed in 
Sections IVIII Bl and IVIII PI 

Every time thrusters are used, there tends to be a 
short-term leakage of gas until the valves set (perhaps a 
few days later) . But there can also be long-term leakages 
due to some mechanism which does not quickly correct 
itself. The major Pioneer 10 spin anomaly that marks the 
boundary of Intervals I and II, is a case in point. During 
this interval there was a major factor of ~ 4.5 increase 
in the average spin-rate change to ~ (—0.0861 ± 0.0009) 
rpm/yr. One also notices kinks during the interval. 

Few values of the Pioneer 10 spin rate were obtained 



after mid- 1993, so the long-term spin-rate change is not 
well-determined in Interval III. But from what was mea- 
sured, there was first a short-term transition region of 
about a year where the spin-rate change was ^ —0.0160 
rpm/yr. Then things settled down to a spin-rate change 
of about ^ (—0.0073 ± 0.0015) rpm/yr, which is small 
and less than that of interval I. 

The effects of the maneuvers on the values of ap will 
allow an estimation of the gas leak systematic in Section 
IVIII PI Note, however, that in the time periods studied, 
only orientation maneuvers were made, not trajectory 
maneuvers. 

Shortly after Pioneer 11 was launched on 5 April 1973, 
the spin period was 4.845 s. A spin precession maneuver 
on 18 May 1973 reduced the period to 4.78 s and after- 
wards, because of a series of precession maneuvers, the 
period lengthened until it reached 5.045 s at encounter 
with Jupiter in December 1974. The period was fairly 
constant until 18 December 1976, when a mid-course 
maneuver placed the spacecraft on a Saturn-encounter 
trajectory. Before the maneuver the period was 5.455 s, 
while after the maneuver it was 7.658 s. At Saturn en- 
counter in December 1979 the period was 7.644 s, little 
changed over the three-year post maneuver cruise phase. 
At the start of our data interval on 5 January 1987, the 
period was 7.321 s, while at the end of the data interval 
in October 1990 it was 7.238 s. 

Although the linear fit to the Pioneer 11 spin rate 
shown in Figure 1121 is similar to that for Pioneer 10 in 
Interval I, ~ (— 0.0234 ± 0.0003) rpm/yr, the causes ap- 
pear to be very different. (Remember, although iden- 
tical in design, Pioneers 10 and 11 were not identical 
in quality 14J.) Unlike Pioneer 10, the spin period for 
Pioneer 11 was primarily affected at the time of spin pre- 
cession maneuvers. One sees that at maneuvers the spin 
period decreases very quickly, while in between maneu- 
vers the spin rate actually tends to increase at a rate of 
~ (-(-0.0073 ± 0.0003) rpm/yr (perhaps due to a gas leak 
in the opposite direction). 

All the above observations aid us in the interpretation 
of systematics in the following three sections. 



B. Recent results using JPL software 

The latest results from JPL are based on an upgrade, 
Sigma, to JPL's ODP software 0|. Sigma, developed for 
NASA's Cassini Mission to Saturn, eliminates structural 
restrictions on memory and architecture that were im- 
posed 30 years ago when JPL space navigation depended 
solely on a Univac 1108 mainframe computer. Five ODP 
programs and their interconnecting files have been re- 
placed by the single program Sigma to support filtering, 
smoothing, and mapping functions. 

We used Sigma to reduce the Pioneer 10 (in three time 
intervals) and 11 Doppler of the unmodeled acceleration, 
ap, along the spacecraft spin axis. As mentioned, the 
Pioneer 10 data interval was extended to cover the to- 
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TABLE I: Determinations of ap in units of 10~** cm/s^ from the three time intervals of Pioneer 10 data and from Pioneer 
11. As described in the text, results from various ODP / Sigma and CHASMP calculations are listed. For ODP / Sigma, "WLS" 
signifies a weighted least-squares calculation, which was used with i) no solar corona model and ii) the 'Cassini' solar corona 
model. Also for ODP / Sigma, "BSF" signifies a batch-sequential filter calculation, which was done with iii) the 'Cassini' solar 
corona model. Further (see Section llX CI . a 1-day batch-sequential estimation for the entire data interval of 11.5 years for 
Pioneer 10 yielded a result ap = (7.77 ± 0.16) x 10"** cm/s^. The CHASMP calculations were aU WLS. These calculations were 
done with i) no solar corona model, ii) the 'Cassini' solar corona model, iii) the 'Cassini' solar corona model with corona data 
weighting and F10.7 time variation calibration. Note that the errors given are only formal calculational errors. The much 
larger deviations of the results from each other indicate the sizes of the systematics that are involved. 



Program/Estimation method 


Pio 10 (I) 


Pio 10 (II) 


Pio 10 (III) 


Pio 11 


Sigma, WLS, 

no solar corona model 


8.02 ±0.01 


8.65 ±0.01 


7.83 ±0.01 


8.46 ±0.04 


Sigma, WLS, 

with solar corona model 


8.00 ±0.01 


8.66 ±0.01 


7.84 ±0.01 


8.44 ± 0.04 


Sigma, BSF, 1-day batch, 
with solar corona model 


7.82 ±0.29 


8.16 ±0.40 


7.59 ±0.22 


8.49 ±0.33 


CHASMP, WLS, 

no solar corona model 


8.25 ±0.02 


8.86 ±0.02 


7.85 ±0.01 


8.71 ±0.03 


CHASMP, WLS, 

with solar corona model 


8.22 ±0.02 


8.89 ±0.02 


7.92 ±0.01 


8.69 ±0.03 


CHASMP, WLS, with 
corona, weighting, and F10.7 


8.25 ±0.03 


8.90 ± 0.03 


7.91 ±0.01 


8.91 ±0.04 



tal time interval 3 January 1987 to 22 July 1998. Of 
the total data set of 20,055 Pioneer 10 Doppler points, 
JPL used ~19,403, depending on the initial conditions 
and editing for a particular run. Of the available 10,616 
(mainly shorter time-averaged) Pioneer 11 data points, 
10,252 were used (4919 two-way and 5333 three-way). 

We wanted to produce independent (i.e., uncorrelated) 
solutions for ap in the three Pioneer 10 segments of data. 
The word independent solution in our approach means 
only the fact that data from any of the three segments 
must not have any information (in any form) passed 
onto it from the other two intervals while estimating the 
anomaly. We moved the epoch from the beginning of 
one data interval to the next by numerically integrating 
the equations of motion and not iterating on the data 
to obtain a better initial conditions for this consequent 
segment. Note that this numerical iteration provided us 
only with an a priori estimate for the initial conditions 
for the data interval in question. 

Other parameters included in the fitting model were 
the six spacecraft heliocentric position and velocity co- 
ordinates at the 1987 epoch of 1 January 1987, 01:00:00 
ET, and 84 (i.e., 28 x 3) instantaneous velocity increments 
along the three spacecraft axes for 28 spacecraft attitude 
(or spin orientation) maneuvers. If these orientation ma- 
neuvers had been performed at exactly six month inter- 
vals, there would have been 23 maneuvers over our 11.5 
year data interval. But in fact, five more maneuvers were 
performed than expected over this 11.5 year interval giv- 
ing a total of 28 maneuvers in all. 

As noted previously, in fitting the Pioneer 10 data over 
11.5 years we used the standard space-fixed J2000 coordi- 
nate system with planetary ephemeris DE405, referenced 



to ICRF. The three-dimensional locations of the track- 
ing stations in the Earth's body-fixed coordinate system 
(geocentric radius, latitude, longitude) were taken from a 
set recommended by ICRF for JPL's DE405. The time- 
varying Earth orientation in J2000 coordinates was de- 
fined by a 1998 version of JPL's EOF file. This accounted 
for the geophysical motion of the Earth's pole with re- 
spect to its spin axis and the Earth's time varying spin 
rate. 

JPL used both the weighted least-squares (WLS) and 
the batch-sequential filter (BSF) algorithms for the final 
calculations. In the first three rows of Tableware shown 
the ODP results for i) WLS with no corona, ii) WLS with 
the Cassini corona model, and iii) BSF with the Cassini 
corona model. 

Observe that the WLS acceleration values for Pioneer 
10 in Intervals I, II, and III are larger or smaller, re- 
spectively, just as the spin-rate changes in these inter- 
vals are larger or smaller, respectively. This indicates 
that the small deviations may be due to a correlation 
with the large gas leak/spin anomaly. We will argue this 
quantitatively in Section IVIII Fl For now we just note 
that we therefore expect the number from Interval III, 
ap = 7.83 X 10~®cm/s^, to be close to our basic (least 
perturbed) JPL result for Pioneer 10. We also note that 
the statistical errors and the effect of the solar corona 
are both small for WLS, and will be handled in our error 
budget. 

In FigurelTSIwe show ODP / Sigma WLS Doppler residu- 
als for the entire Pioneer 10 data set. The residuals were 
obtained by first solving for ap with no corona in each of 
the three Now look at the batch-sequential results in row 
3 of Table m First, note that the statistical Intervals in- 



25 



dependently and then subtracting these solutions (given 
in Table nj from the fits within the corresponding data 
intervals. 




FIG. 13: ODP Doppler residuals in Hz for the entire Pioneer 

10 data span. The two solid vertical lines in the upper part of 
the plot indicate the boundaries between data Intervals I/II 
and II/III, respectively. Maneuver times are indicated by the 
vertical dashed lines in the lower part of the plot. 

One can easily see the very close agreement with the 
CHASMP residuals of Figure which go up to 14 De- 
cember 1994. 

The Pioneer 11 number is significantly higher. A de- 
viation is not totally unexpected since the data was rela- 
tively noisy, was from much closer in to the Sun, and was 
taken during a period of high solar activity. We also do 
not have the same handle on spin-rate change effects as 
we did for Pioneer 10. We must simply take the number 
for what it is, and give the basic JPL result for Pioneer 

11 as ap = 8.46 x 10"^ cm/s^. 

Now look at the batch-sequential results in row 3 of Ta- 
ble ^ First, note that the statistical errors are an order 
of magnitude larger than for WLS. This is not surprising 
since: i) the process noise significantly affects the preci- 
sion, ii) BSF smoothes the data and the data from the 
various intervals is more correlated than in WLS. The ef- 
fects of all this are that all four numbers change so as to 
make them all closer to each other, but yet all the num- 
bers vary by less than 2a from their WLS counterparts. 

Finally, there is the annual term. It remains in the data 
(for both Pioneers 10 and 11). A representation of it can 
be seen in a 1-day batch-sequential averaged over all 11.5 
years. It yielded a result ap = (7.77±0.16) x 10"*^ cm/s^, 
consistent with the other numbers/errors, but with an 
added annual oscillation. In the following subsection we 
will compare JPL results showing the annual term with 
the counterpart Aerospace results. 

We will argue in Section IIX CI that this annual term 
is due to the inability to model the angles of the Pio- 
neers' orbits accurately enough. [Note that this annual 
term is not to be confused with a small oscillation seen 
in Figure IHlthat can be caused by mispointing towards 
the spacecraft by the fit programs.] 



C. Recent results using The Aerospace 
Corporation software 

As part of an ongoing upgrade to CHASMP 's accuracy. 
Aerospace has used Pioneer 10 and 11 as a test bed to 
confirm the revision's improvement. In accordance with 
the JPL results of Section IVl Bl we used the new version 
of CHASMP to concentrate on the Pioneer 10 and 11 
data. The physical models are basically the same ones 
that JPL used, but the techniques and methods used are 
largely different. (See Section liXBI ) 

The new results from the Aerospace Corporation's 
software are based on first improving the Planetary 
Ephemeris and Earth orientation and spacecraft spin 
models required by the program. That is: i) the spin 
data file has been included with full detail; ii) a newer 
JPL Earth Orientation Parameters file was used; iii) all 
lERS tidal terms were included; iv) plate tectonics were 
included; v) DE405 was used; vi) no a priori information 
on the solved for parameters was included in the fit; vii) 
Pioneer 11 was considered, viii) the Pioneer 10 data set 
used was extended to 14 Feb. 1998. Then the Doppler 
data was refitted. 

Beginning with this last point: CHASMP uses the 
same original data file, but it performs an additional data 
compression. This compression combines the longest con- 
tiguous data composed of adjacent data intervals or data 
spans with duration > 600 s (effectively it prefers 600 
and 1980 second data intervals). It ignores short-time 
data points. Also, Aerospace uses an N-cr/fixed bound- 
ary rejection criteria that rejects all data in the fit with a 
residual greater than ±0.025 Hz. These rejection criteria 
resulted in the loss of about 10 % of the original data 
for both Pioneers 10 and 11. In particular, the last five 
months of Pioneer 10 data, which was all of data-lengths 
less than 600 s, was ignored. Once these data compres- 
sion/cuts were made, CHASMP used 10,499 of its 11,610 
data points for Pioneer 10 and 4,380 of its 5,137 data 
points for Pioneer 11. 

Because of the spin-anomaly in the Pioneer 10 data, 
the data arc was also divided into three time intervals 
(although the I/II boundary was taken as 31 August 
1990 85] ). In what was especially useful, the Aerospace 
analysis uses direct propagation of the trajectory data 
and solves for the parameter of interest only for the data 
within a particular data interval. That means the three 
interval results were truly independent. Pioneer 11 was 
fit as a single arc. 

Three types of runs are listed, with: i) no corona; ii) 
with Cassini corona model of Sections IIV Dl and IVII CI 
and iii) with the Cassini corona model, but added are 
corona data weighting (Section IIVD|I and the time- 
variation called "F10.7" ;65]. (The number 10.7 labels 
the wavelength of solar radiation, A=10.7 cm, that, in 
our analysis, is averaged over 81 days.) 

The results are given in rows 4-6 of Table ^ The no 
corona results (row 4) are in good agreement with the 
Sigma results of the first row. This is especially true 



for the extended-time Interval III values for Pioneer 10, 
which interval had clean data. However there is more 
disagreement with the values for Pioneer 10 in Intervals 
I and II and for Pioneer 11. These three data sets all 
were noisy and underwent more data-editing. Therefore, 
it is significant that the deviations between Sigma and 
CHASMP in these arcs are all similar, but small, be- 
tween 0.20 to 0.25 of our units. As before, the effect 
of the solar corona is small, even with the various model 
variations. But most important, the numbers from Sigma 
and CHASMP for Pioneer 10 Interval III are in excellent 
agreement. 

Further, CHASMP also found the annual term. (Re- 
call that CHASMP can also look for a temporal varia- 
tion by calculating short time averages.) Results on the 
time variation in ap can be seen in Figure [HI Although 
there could possibly be ap variations of ±2 x lO^'* cm/s^ 
on a 200-day time scale, a comparison of the variations 
with the error limits shown in Figure lHI indicate that our 
measurements of these variations are not statistically sig- 
nificant. The 5-day averages of ap from ODP (using the 
batch-sequential method) are not reliable at solar con- 
junction in the middle (June) of each year, and hence 
should be ignored there. The CHASMP 200-day aver- 
ages suppress the solar conjunction bias inherent in the 
ODP 5-day averages, and they reliably indicate a con- 
stant value of flp. Most encouraging, these results clearly 
indicate that the obtained solution is consistent, stable, 
and its mean value does not strongly depend on the esti- 
mation procedure used. The presence of the small annual 
term on top of the obtained solution is apparent. 



D. Our solution, before systematics, for the 
anomalous acceleration 

From Table ID we can intuitively draw a number of con- 
clusions: 

A) The effect of the corona is small. This systematic 
will be analyzed in Section FV 11 HI 

B) The numerical error is small. This systematic will 
be analyzed in Section llX Al 

C) The differences between the Sigma and CHASMP 
Pioneer 10 results for Interval I and Interval II, respec- 
tively, we attribute to two main causes: especially i) the 
different data rejection techniques of the two analyses but 
also ii) the different maneuver simulations. Both of these 
effects were especially significant in Interval II, where the 
data arc was small and a large amount of noisy data was 
present. Also, to account for the discontinuity in the spin 
data that occurred on 28 January 1992 (see Figure [TT|l . 
Aerospace introduced a fictitious maneuver for this in- 
terval. Even so, the deviation in the two values of ap 
was relatively small, namely 0.23 and 0.21, respectively, 
xlO~^ cm/s^. 

D) The changes in ap in the different Intervals, corre- 
lated with the changes in spin-rate change, are likely (at 
least partially) due to gas leakage. This will be discussed 
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Anomaious Acceleration From ODP and CHASMP 
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FIG. 14: Consistency of the ODP /Stgma and CHASMP time- 
variation signals. The dots show 5-day sample averages of 
the anomalous acceleration of Pioneer 10 from ODP /Sigma 
using BSF with a 200-day correlation time. From this data, 
the solid lines show the mean values of ap in the three inter- 
vals corresponding to the three separate spin down histories. 
The dashed lines represent the large batch-sequential compu- 
tational error bounds on the three values of ap. The 200- 
day acceleration values using CHASMP are the solid squares. 
At the time positions where there are CHASMP results, the 
agreement between the CHASMP and the ODP/ Sigma results 
is clear. 

in Section rvTHFl 

But independent of the origin, this last correlation be- 
tween shifts in ap and changes in spin rate actually allows 
us to calculate the best "experimental" base number for 
Pioneer 10. To do this, assume that the spin-rate change 
is directly contributing to an anomalous acceleration off- 
set. Mathematically, this is saying that in any interval 
i — I, II, HI, for which the spin-rate change is an approx- 
imate constant, one has 

ap(^') = ap(o) - K (19) 

where k is a constant with units of length and ap(o) = 
ap{9 = 0) is the Pioneer acceleration without any spin- 
rate change. 

One now can fit the data to Eq. (|19|l to obtain solutions 
for K and ap(o). The three intervals i — I, II, HI provide 
three data combinations {ap(i)(6'), di}. We take our base 
number, with which to reference systematics, to be the 
weighted average of the Sigma and CHASMP results for 
a p(o) when no corona model was used. Start first with the 
Sigma Pioneer 10 solutions in row one of Tableland the 
Pioneer 10 spin-down rates given in Section rVI Bl and Fig- 
ureCH ap(^^ = (8.02 ±0.01, 8.65 ±0.01, 7.83 ±0.01) in 

units of 10-8 cm/s2 and 6, = -(0.0181±0.0001, 0.0861± 
0.0009, 0.0073 ± 0.0015) in units of rpm/yr, where 

1 rpm/ycar = 5.281 x 10"^" rev/s^ 

= 3.318 X 10^^ radians/s^ (20) 



27 



With these data we use the maximum hkehhood 
and minimum variance approach to find the optimally 
weighted least-squares solution for ap(o): 



a 



Sigma 
P(0) 



(7.82 ±0.01) X 10~^ cm/s^ 



(21) 



with solution for the parameter k obtained as K^^sma _ 
(29.2 ± 0.7) cm. Similarly, for CHASMP one takes the 
values for ap from row four of TableHl afl^^^ = (8.25 ± 
0.02, 8.86 ± 0.02, 7.85 ± 0.01) and uses them with the 
same 9i as above. The solution for ap(o) in this case is 



^P(O) 



(7.89 ±0.02) X 10"» cm/s 



(22) 



together with k™^"^ = (34.7 ± 1.1) cm. The solutions 
for Sigma and CHASMP are similar, 7.82 and 7.89 in our 
units. We take the weighted average of these two to yield 
our base line "experimental" number for ap: 



•'P(expe 



(7.84 ±0.01) X 10"« cm/s\ (23) 



[The weighted average constant k is kq — (30.7 ± 0.6) 
cm.] 

For Pioneer 11, we only have the one 3| year data arc. 
The weighted average of the two programs' no corona re- 
sults is (8.62±0.02)x 10^* cm/s^. We observed in Section 
IVI Al that between maneuvers (which are accounted for 
- see Section llV Ep there is actually a spin rate increase 
of ^ (+0.0073 ± 0.0003) rpm/yr. If one uses this spin-up 
rate and the Pioneer 10 value for kq = 30.7 cm given 
above, one obtains a spin-rate change corrected value for 
ap. We take this as the experimental value for Pioneer 
11: 



a 



P(exper) 



.55 ± 0.02) X 10"^ cm/s . (24) 



VII. SOURCES OF SYSTEMATIC ERROR 
EXTERNAL TO THE SPACECRAFT 

We are concerned with possible systematic acceleration 
errors that could account for the unexplained anomalous 
acceleration directed toward the Sun. There exist de- 
tailed publications describing analytic recipes developed 
to account for non-gravitational accelerations acting on 
spacecraft. (For a summary see Milani et al. H^].) With 
regard to the specific Pioneer spacecraft, possible sources 
of systematic acceleration have been discussed before for 
Pioneer 10 and 11 at Jupiter [s^l and Pioneer 11 at Sat- 
urn (67j . 

External forces can produce three vector components 
of spacecraft acceleration, unlike forces generated on 
board the spacecraft, where the two non-radial compo- 
nents (i.e., those that are eifectively perpendicular to 
the spacecraft spin) are canceled out by spacecraft ro- 
tation. However, non-radial spacecraft accelerations are 
difficult to observe by the Doppler technique, which mea- 
sures spacecraft velocity along the Earth-spacecraft line 



of sight. But with several years of Doppler data, it is in 
principle possible to detect systematic non-radial accel- 
eration components |73j . 

With our present analysis [t^ wc find that the Doppler 
data yields only one significant component of unmodeled 
acceleration, and that any acceleration components per- 
pendicular to the spin axis are small. This is because 
in the fitting we tried including three unmodeled accel- 
eration constants along the three spacecraft axes (spin 
axis and two orthogonal axes perpendicular to the spin 
axis). The components perpendicular to the spin axis 
had values consistent with zero to a 1-cr accuracy of 2 
X 10"^ cm/s^ and the radial component was equal to 
the reported anomalous acceleration. Further, the ra- 
dial acceleration was not correlated with the other two 
unmodeled acceleration components. 

Although one could in principle set up complicated en- 
gineering models to predict all or each of the systemat- 
ics, often the uncertainty of the models is too large to 
make them useful, despite the significant efi^ort required. 
A different approach is to accept our ignorance about a 
non-gravitational acceleration and assess to what extent 
these can be assumed a constant bias over the time scale 
of all or part of the mission. (In fact, a constant ac- 
celeration produces a linear frequency drift that can be 
accounted for in the data analysis by a single unknown 
parameter.) In fact, we will use both approaches. 

In most orbit determination programs some effects, like 
the solar radiation pressure, are included in the set of 
routinely estimated parameters. Nevertheless we want 
to demonstrate their influence on Pioneer's navigation 
from the general physics standpoint. This is not only to 
validate our results, but also to be a model as to how 
to study the influence of the other physical phenomena 
that are not yet included in the standard navigational 
packages for future more demanding missions. Such mis- 
sions will involve either spacecraft that will be distant 
or spacecraft at shorter distances where high-precision 
spacecraft navigation will be required. 

In this section we will discuss possible systematics (in- 
cluding forces) generated external to the spacecraft which 
might signiflcantly affect our results. These start with 
true forces due to (1) solar-radiation pressure and (2) so- 
lar wind pressure. We go on to discuss (3) the effect of the 
solar corona and its mismodeling, (4) electro-magnetic 
Lorentz forces, (5) the influence of the Kuiper belt, (6) 
the phase stability of the reference atomic clocks, and 
(7) the mechanical and phase stability of the DSN an- 
tennae, together with influence of the station locations 
and troposphere and ionosphere contributions. 



A. Direct solar radiation pressure and mass 

There is an exchange of momentum when solar pho- 
tons impact the spacecraft and are either absorbed or 
reflected. Models for this solar pressure effect were de- 
veloped before either Pioneer 10 or 11 were launched [s^ 
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and have been refined since then. The models take into 
account various parts of the spacecraft exposed to solar 
radiation, primarily the high-gain antenna. It computes 
an acceleration directed away from the Sun as a function 
of spacecraft orientation and solar distance. 

The models for the acceleration due to solar radiation 
can be formulated as 



as.p.(r) 



ICfQAcose{r) 
c M r2 ' 



(25) 



/© = 1367 W/m^(AU)^ is the (effective-temperature 
Stefan-Boltzmann) "solar radiation constant" at 1 AU 
from the Sun and A is the effective size of the craft as 
seen by the Sun (For Pioneer the area was taken 

to be the antenna dish of radius 1.73 m.) 6 is the an- 
gle between the axis of the antenna and the direction of 
the Sun, c is the speed of light, M is the mass of the 
spacecraft (taken to be 251.883 for Pioneer 10), and r 
is the distance from the Sun to the spacecraft in AU. JC 
[90l | is the effective absorption/reflection coefficient. 
For Pioneer 10 the simplest approximately correct model 
yields /Co — 1.71 joj- Eq. provides a good model 

for analysis of the effect of solar radiation pressure on 
the motion of distant spacecraft and is accounted for by 
most of the programs used for orbit determination. 

However, in reality the absorptivities, emissivities, and 
effective areas of spacecraft parts parameters which, al- 
though modeled by design, are determined by calibration 
early in the mission |92| . One determines the magnitude 
of the solar-pressure acceleration at various orientations 
using Doppler data. (The solar pressure effect can be 
distinguished from gravity's 1 /r^ law because cos 9 varies 
pl5j.) The complicated set of program input parameters 
that yield the parameters in Eq. H25|l are then set for 
later use 592j. Such a determination of the parameters 
for Pioneer 10 was done, soon after launch and later. 
When applied to the solar radiation acceleration in the 
region of Jupiter, this yields (from a 5 % uncertainty in 
a..p. 113) 



.(r = 5.2AU) = (70.0±3.5) X 10"^ cm/s , 
/Cs.2 = 1.77. 



(26) 



The second of Eqs. (|26|l comes from putting the first into 
Eq. H25|l . Note, specifically, that in a fit a too high input 
mass will be compensated for by a higher effective K.. 
Because of the 1/r^ law, by the time the craft reached 

10 AU the solar radiation acceleration was 18.9 x 10~^ 
cm/s^ going down to 0.39 of those units by 70 AU. Since 
this systematic falls off as r~^, it can bias the Doppler de- 
termination of a constant acceleration at some level, even 
though most of the systematic is correctly modeled by 
the program itself. By taking the average of the ac- 
celeration curves over the Pioneer distance intervals, we 
estimate that the systematic error from solar-radiation 
pressure in units of 10^* cm/s^ is 0.001 for Pioneer 10 
over an interval from 40 to 70 AU, and 0.006 for Pioneer 

11 over an interval from 22 to 32 AU. 



However, this small uncertainty is not our main prob- 
lem. In actuality, since the parameters were fit the mass 
has decreased with the consumption of propellant. Effec- 
tively, the 1/r^ systematic has changed its normalization 
with time. If not corrected for, the difference between 
the original l/r^ and the corrected 1/r^ will be inter- 
preted as a bias in ap. Unfortunately, exact information 
on gas usage is unavailable 16] . Therefore, in dealing 
with the effect of the temporal mass variation during the 
entire data span (i.e. nominal input mass vs. actual mass 
history J£] ) we have to address two effects on the so- 
lutions for the anomalous acceleration ap. They are i) 
the effect of mass variation from gas consumption and ii) 
the effect of an incorrect input mass 0, 0| . 

To resolve the issue of mass variation uncertainty we 
performed a sensitivity analysis of our solutions to dif- 
ferent spacecraft input masses. We simply re-did the no- 
corona, WLS runs of Table U with a range of different 
masses. The initial wet weight of the package was 259 kg 
with about 36 kg of consumable propellant. For Pioneer 
10, the input mass in the program fit was 251.883 kg, 
roughly corresponding to the mass after spin-down. By 
our data period, roughly half the fuel (18 kg) was gone 
so we take 241 kg as our nominal Pioneer 10 mass. Thus, 
the effect of going from 251.883 kg to 241 kg we take to 
be our bias correction for Pioneer 10. We take the un- 
certainty to be given by one half the effect of going from 
plus to minus 9 kg (plus or minus a quarter tank) from 
the nominal mass of 241 kg. 

For the three intervals of Pioneer 10 data, using 
ODP / Sigma yields the following changes in the accel- 
erations: 



= [(0.040 ± 0.035), (0.029 ±0.025), 

(0.020 ±0.017)] X 10"* cm/s^ 



As expected, these results make ap larger. For our sys- 
tematic bias we take the weighted average of 6a'^^^ for 
the three intervals of Pioneer 10. The end result is 



,.p. = (0.03 ± 0.01) X 10"* cm/s 



(27) 



For Pioneer 11 we did the same except our bias point 
was 3/4 of the fuel gone (232 kg). Therefore the bias 
results by going from the input mass of 239.73 to 232 kg. 
The uncertainty is again defined by ± 9 kg. The result 
for Pioneer 11 is more sensitive to mass changes, and we 
find using OT)P / Sigma 



(0.09 ± 0.21) X 10"* cm/s^ 



(28) 



The bias number is three times larger than the similar 
number for Pioneer 10, and the uncertainty much larger. 
We return to this difference in Section IVlII Gl 

The previous analysis also allowed us to perform 
consistency checks on the effective values of K, which 
the programs were using. By taking [rminTmax]"^ — 
[/ {dr jr^ )/ / dr] for the inverse distance squared of a data 
set, varying the masses, and determining the shifts in ap 
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'^Pio-lO(I) " 

and ^-f_- 



we could determine the values of K. implied, We found: 
1.72; IC°Z-n « 1-82; /Cp=ff„To(i) « 1-74; 
1.84. [The hat over the last /C indicates 
it was multiplied by (237.73/251.883) because CHASMP 
uses 259.883 kg instead of 239.73 kg for the input mass.] 
All these values of /C are in the region expected and are 
clustered around the value /C5.2 in Eq. (|26|l . 

Finally, if you take the average values of /C for Pio- 
neers 10 and 11 (1.73, 1.83), multiply these numbers by 
the input masses (251.883, 239.73) kg, and divide them 
by our nominal masses (241, 232) kg, you obtain (1.87, 
1.89), indicating our choice of nominal masses was well 
motivated. 



B. The solar wind 

The acceleration caused by the solar wind has the same 
form as Eq. H25II . with /q replaced by nipV^n, where 
n « 5 cm"'^ is the proton density at 1 AU and v « 400 
km/s is the speed of the wind. Thus, 



JCs 



A cos t 



cM r2 



1.24 X 10"" 



20 AU 



cm/s^ (29) 



Because the density can change by as much as 100%, 
the exact acceleration is unpredictable. But there are 
measurements |s3 showing that it is about 10"'"' times 
smaller than the direct solar radiation pressure. Even if 
we make the very conservative assumption that the solar 
wind contributes only 100 times less force than the solar 
radiation, its smaller contribution is completely negligi- 
ble. 



C. The effects of the solar corona and models of it 



As we saw in the previous Section FVII Bl the effect of 
the solar wind pressure is negligible for distant spacecraft 
motion in the solar system. However, the solar corona 
effect on propagation of radio waves between the Earth 
and the spacecraft needs to be analyzed in more detail. 

Initially, to study the sensitivity of ap to the solar 
corona model, we were also solving for the solar corona 
parameters A, B, and C of Eq. © in addition to ap. 
However, we realized that the Pioneer Doppler data is 
not precise enough to produce credible results for these 
physical parameters. In particular, we found that solu- 
tions could yield a value of ap which was changed by 
of order 10 % even though it gave unphysical values of 
the parameters (especially B, which previously had been 
poorly defined even by the Ulysses mission [s^). [By 
"unphysical" we mean electron densities that were either 
negative or positive with values that are vastly different 
from what would be expected.] 



Therefore, as noted in Sect ion 11 VDI we decided to use 
the newly obtained values for A, B, and C from the 
Cassini mission and use them as inputs for our analy- 
ses: A = 6.0 X 10^,5 = 2.0 X 10'',C = 0.6 x 10^ all in 
meters 64']. This is the "Cassini corona model." 

The effect of the solar corona is expected to be small 
for Doppler and large for range. Indeed it is small for 
Sigma. For ODP / Sigma, the time-averaged effect of the 
corona was small, of order 



= ±0.02 X 10"* cm/s , 



(30) 



as might be expected. We take this number to be the 
error due to the corona. 

What about the results from CHASMP. Both analyses 
use the same physical model for the effect of the steady- 
state solar corona on radio- wave propagation through the 
solar plasma (that is given by Eq. H10|l ). However, there 
is a slight difference in the actual implementation of the 
model in the two codes. 

ODP calculates the corona effect only when the Sun- 
spacecraft separation angle as seen from the Earth (or 
Sun-Earth-spacecraft angle) is less then tt/2. It sets the 
corona contribution to zero in all other cases. Earlier 
CHASMP used the same model and got a small corona 
effect. Presently CHASMP calculates an approximate 
corona contribution for all the trajectory. Specific atten- 
tion is given to the region when the spacecraft is at oppo- 
sition from the Sun and the Sun-Earth-spacecraft angle 
~ TT. There CHASMP 's implementation truncates the 
code approximation to the scaling factor F in Eq. H10|l . 
This is specifically done to remove the fictitious diver- 
gence in the region where "impact parameter" is small. 

However, both this and also the more complicated 
corona models (with data- weighting and/or "F10.7" time 
variation) used by CHASMP produce small deviations 
from the no-corona results. Our decision was to incorpo- 
rate these small deviations between the two results due 
to corona modeling into our overall error budget as a 
separate item: 



^corona_niodel 



±0.02 X 10"® cm/s^ 



(31) 



This number could be discussed in Section lIXI on com- 
putational systematics. Indeed, that is where it will be 
listed in our error budget. 



D. Electro-magnetic Lorentz forces 

The possibility that the spacecraft could hold a charge, 
and be deflected in its trajectory by Lorentz forces, was 
a concern for the magnetic field strengths at Jupiter and 
Saturn. However, the magnetic field strength in the outer 
solar system is on the order of < 1 7 (7 = 10"^ Gauss). 
This is about a factor of 10^ times smaller than the mag- 
netic field strengths measured by the Pioneers at their 
nearest approaches to Jupiter: 0.185 Gauss for Pioneer 
10 and 1.135 Gauss for the closer in Pioneer 11 l93l. 
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Also, there is an upper limit to the charge that a space- 
craft can hold. For the Pioneers that limit produced an 
upper bound on the Lorentz acceleration at closest ap- 
proach to Jupiter of 20 x 10"^ cm/s^ [13 ■ With the 
interplanetary field being so much lower than at Jupiter, 
we conclude that the electro-magnetic force on the Pio- 
neer spacecraft in the outer solar system is at worst on 
the order of 10~^^ cm/s^, completely negligible [93 |. 

Similarly, the magnetic torques acting on the space- 
craft were about a factor of 10~^ times smaller than 
those acting on Earth satellites, where they are a con- 
cern. Therefore, for the Pioneers any observed changes 
in spacecraft spin cannot be caused by magnetic torques. 

E. The Kuiper belt's gravity 



effect. (See Figure El) Further, the accelerations are 
not constant across the data range. Rather, they show 
an increasing effect as Pioneer 10 approaches the belt 
and a decreasing effect as Pioneer 10 recedes from the 
belt, even with a uniform density model. For these two 
reasons, we excluded the dust belt as a source for the 
Pioneer effect. 

More recent infrared observations have ruled out more 
than 0.3 Earth mass of Kuiper Belt dust in the trans- 
Neptunian region jol, [o^ • Therefore, we can now place 
a limit of ±3 x 10~^° cm/s^ for the contribution of the 
Kuiper belt. 

Finally, we note that searches for gravitational encoun- 
ters of Pioneer with large Kuiper-belt objects have so far 
not been successful EM. 



From the study of the resonance effect of Neptune upon 
Pluto, two primary mass concentration resonances of 3:2 
and 2:1 were discovered [o^, corresponding to 39.4 AU 
and 47.8 AU, respectively. Previously, Boss and Peale 
had derived a model for a non-uniform density distribu- 
tion in the form of an infinitesimally thin disc extending 
from 30 AU to 100 AU in the eclipticplane We 
combined the results of Refs. and [o^ to determine 
if the matter in the Kuiper belt could be the source of 
the anomalous acceleration of Pioneer 10 [97l |. 

We specifically studied three distributions, namely: i) 
a uniform distribution, ii) a 2:1 resonance distribution 
with a peak at 47.8 AU, and iii) a 3:2 resonance distri- 
bution with a peak at 39.4 AU. Figure [T51 exhibits the 
resulting acceleration felt by Pioneer 10, from 30 to 65 
AU which encompassed our data set at the time. 
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FIG. 15: Possible acceleration caused by dust in the Kuiper 
belt. 

We assumed a total mass of one Earth mass, which is 
significantly larger than standard estimates. Even so, the 
accelerations are only on the order of 10~^ cm/s^, which 
is two orders of magnitude smaller than the observed 



F. Phase and frequency stability of clocks 

After traversing the mechanical components of the an- 
tenna, the radio signal enters the DSN antenna feed and 
passes through a series of amplifiers, filters, and cables. 
Averaged over many experiments, the net effect of this 
on the calculated dynamical parameters of a spacecraft 
should be very small. We expect instrumental calibra- 
tion instabilities to contribute 0.2 x 10"^ cm/s^ to the 
anomalous acceleration on a 6 s t ime interval. Thus, 
in order for the atomic clocks |lOlj | to have caused the 
Pioneer effect, all the atomic clocks used for signal refer- 
encing clocks would have had to have drifted in the same 
manner as the local DSN clocks. 

In Section^we observed that without using the appar- 
ent anomalous acceleration^ the CHASMP residuals show 
a steady frequency drift |30j of about —6 x 10~^ Hz/s, 
or 1.5 Hz over 8 years (one-way only). This equates to 
a clock acceleration, —at, of —2.8 x 10~^* s/s^. (See Eq. 
(|16|1 and Figure |H1) To verify that it is actually not the 
clocks that are drifting, we analyzed the calibration of 
the frequency standards used in the DSN complex. 

The calibration system itself is referenced to Hydro- 
gen maser atomic clocks. Instabilities in these clocks 
are another source of instrumental error which needs 
to be addressed. The local reference is synchronized 
to the frequency standards generated either at the Na- 
tional Institute of Standards and Technology (NIST), lo- 
cated in Boulder, Colorado or at the U. S. Naval Obser- 
vatory (USNO), Washington, DC. These standards are 
presently distributed to local stations by the Global Po- 
sitioning System (GPS) satellites. [During the pre-GPS 
era, the station clocks used signals from WWV to set the 
Cesium or Hydrogen masers. WWV, the radio station 
which broadcasts time and frequency services, is located 
in Fort Collins, CO.] While on a track, the station is 
"free- running," i.e., the frequency and timing data are 
generated locally at the station. The Allan variances are 
about 10"^'^ for Cesium and 10"^^ for Hydrogen masers. 
Therefore, over the data-pass time interval, the data ac- 
curacy is on the order of one part in 1000 GHz or better. 
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Long-term frequency stability tests are conducted with 
the exciter/transmitter subsystems and the DSN's radio- 
science open-loop subsystem. An uplink signal gener- 
ated by the exciter is translated at the antenna by a test 
translator to a downlink frequency. (See Section UTTl ') 
The downlink signal is then passed through the RF-IF 
downconverter present at the antenna and into the ra- 
dio science receiver chain [3ll | . This technique allows the 
processes to be synchronized in the DSN complex based 
on the frequency standards whose Allan variances are of 
the order of ay ~ 10~^^— 10~^^ for integration time in 
the range from 10 s to 10'^ s. For the S-band frequen- 
cies of the Pioneers, the corresponding Allan variances 
are 1.3 x 10^^^ and 1.0 x 10^^^, respectively, for a 10'^ 
s Doppler integration time. 

Phase-stability testing characterizes stability over very 
short integration times; that is, spurious signals whose 
frequencies are very close to the carrier (frequency). The 
phase noise region is defined to be frequencies within 100 
kHz of the carrier. Both amplitude and phase variations 
appear as phase noise. Phase noise is quoted in dB rela- 
tive to the carrier, in a 1 Hz band at a specified deviation 
from the carrier; for example, dBc-Hz at 10 Hz. Thus, 
for the frequency 1 Hz, the noise level is at —51 dBc and 
10 Hz corresponds to —60 dBc. This was not significant 
for our study. 

Finally, the influence of the clock stability on the de- 
tected acceleration, ap, may be estimated based on the 
reported Allan variances for the clocks, ay. Thus, the 
standard 'single measurement' error on acceleration as 
derived by the time derivative of the Doppler frequency 
data is {cay)/T, where the Allan variance, ay, is calcu- 
lated for 1000 s Doppler integration time, and t is the 
signal averaging time. This formula provides a good rule 
of thumb when the Doppler power spectral density func- 
tion obeys a 1// flicker-noise law, which is approximately 
the case when plasma noise dominates the Doppler error 
budget. Assume a worst case scenario, where only one 
clock was used for the whole 11 years study. (In reality 
each DSN station has its own atomic clock.) To estimate 
the influence of that one clock on the reported accuracy 
of the detected anomaly ap, combine ay — I^v/vq, the 
fractional Doppler frequency shift from the reference fre- 
quency of vq =^ 2.29 GHz, with the estimate for the Al- 
lan variance, ay = 1.3 x 10^^^. This yields a number that 
characterizes the upper limit for a frequency uncertainty 
introduced in a single measurement by the instabilities 
in the atomic clock: cr^ = fof y = 2.98 x lO^'^ Hz for a 
lO'^ Doppler integration time. 

In order to derive an estimate for the total effect, re- 
call that the Doppler observation technique is essentially 
a continuous count of the total number of complete fre- 
quency circles during observational time. Within a year 
one can have as many as N « 3.156 x lO'^ indepen- 
dent single measurements of the clock with duration 10'^ 
seconds. This yields an upper limit for the contribu- 
tion of atomic clock instability on the frequency drift of 
O'clock = ay/^/N w 5.3 X 10^^ Hz/year. But in Section 



IVBI we noted that the observed ap corresponds to a fre- 
quency drift of about 0.2 Hz/year, so the error in ap is 
about 0.0003 x 10^^ cm/s^. Since all data is not inte- 
grated over 1,000 seconds and is data is not available for 
all time, we increase the numerical factor to 0.001, which 
is still negligible to us. [But further, this upper limit 
for the error becomes even smaller if one accounts for 
the number of DSN stations and corresponding atomic 
clocks that were used for the study.] 

Therefore, we conclude that the clocks are not a con- 
tributing factor to the anomalous acceleration at a mean- 
ingfully level. We will return to this issue in Section lXI Dl 
where we will discuss a number of phenomenological time 
models that were used to fit the data. 



G. DSN antennae complex 

The mechanical structures which support the refiecting 
surfaces of the antenna are not perfectly stable. Among 
the numerous effects influencing the DSN antennae per- 
formance, we are only interested in those whose behav- 
ior might contribute to the estimated solutions for ap. 
The largest systematic instability over a long period is 
due to gravit y loa ds and the aging of the structure. As 
discussed in |l02l |. antenna deformations due to gravity 
loads should be absorbed almost entirely into biases of 
the estimated station locations and clock offsets. There- 
fore, they will have little effect on the derived solutions 
for the purposes of spacecraft navigation. 

One can also consider ocean loading, wind loading, 
thermal expansion, and aging of the structure. We 
found none of these can produce the constant drift in 
the Doppler frequency on a time scale comparable to the 
Pioneer data. Also, routine tests are performed by DSN 
personnel on a regular basis to access all the effects that 
may contribute to the overall performance of the DSN 
complex. The information is available and it shows all pa- 
rameters are in the required ranges. Detailed assessments 
of all these effect on th e astrometric VLBI solutions were 
published in Il02j . The results for the astrometric 
errors introduced by the above factors may be directly 
translated to the error budget for the Pioneers, scaled by 
the number of years. It yields a negligible contribution. 

Our analyses also estimated errors introduced by a 
number of station-specific parameters. These include the 
error due to imperfect knowledge in a DSN station lo- 
cation, errors due to troposphere and ionosphere models 
for different stations, and errors due to the Faraday rota- 
tion effects in the Earth's atmosphere. Our analysis indi- 
cates that at most these effects would produce a distance- 
and/or time-dependent drifts that would be easily notice- 
able in the radio Doppler data. What is more important 
is that none of the effects would be able to produce a con- 
stant drift in the Doppler residuals of Pioneers over such a 
long time scale. The updated version of the ODP, Sigma, 
routinely accounts for these error factors. Thus, we run 
covariance analysis for the whole set of these parameters 
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using both Sigma and CHASMP. Based on these stud- 
ies we conclude that mechanical and phase stability of 
the DSN antennae together with geographical locations 
of the antennae, geophysical and atmospheric conditions 
on the antennae site have negligible effects on our solu- 
tions for ap. At most their contributions are at the level 
of (Jdsn < lO^^flp. 



VIII. SOURCES OF SYSTEMATIC ERROR 
INTERNAL TO THE SPACECRAFT 

In this section we will discuss the forces that may be 
generated by spacecraft systems. The mechanisms we 
consider that may contribute to the found constant ac- 
celeration, ap, and that may be caused by the on-board 
mechanisms include: (1) the radio beam reaction force, 
(2) RTG heat reflecting off the spacecraft, (3) differen- 
tial emissivity of the RTGs, (4) non-isotropic radiative 
cooling of the spacecraft, (5) expelled Hehum produced 
within the RTG, (6) thruster gas leakage, and (7) the dif- 
ference in experimental results from the two spacecraft. 



A. Radio beam reaction force 

The Pioneer navigation does not require that the 
spacecraft constantly beam its radio signal, but instead 
it does so only when it is requested to do so from the 
ground control. Nevertheless, the recoil force due to the 
emitted radio-power must also be analyzed. 

The Pioneers have a total nominal emitted radio power 
of eight Watts. It is parameterized as 



de sine vi0), 



(32) 



V^O) being the antenna power distribution. The radiated 
power has been kept constant in time, independent of 
the coverage from ground stations. That is, the radio 
transmitter is always on, even when not received by a 
ground station. 

The recoil from this emitted radiation produces an 
acceleration bias, &rp, on the spacecraft away from the 
Earth of 
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Mc 



(33) 



M is taken to be the Pioneer mass when half the fuel is 
gone [3 • P is the fractional component of the radiation 
momentum that is going in a direction opposite to ap: 



de sine cose v{e). 



(34) 



Ref |j| describes the HGA and shows its downlink an- 
tenna pattern in Fig. 3.6-13. (Thermal antenna expan- 
sion mismodeling is thought to be negligible.) The gain 



is given as (33.3 ± 0.4) dB at zero (peak) degrees. The 
intensity is down by a factor of two (—3 dB) at 1.8 de- 
grees. It is down a factor of 10 (—10 dB) at 2.7 degrees 
and down by a factor of 100 (—20 dB) at 3.75 degrees. 
[The first diffraction minimum is at a little over four de- 
grees.] Therefore, the pattern is a very good conical 
beam. Further, since cos[3.75°] = 0.9978, we can take 
13 = (0.99 ± 0.01), yielding b^p = 1-10. 

Finally, taking the error for the nominal 8 Watts power 
to be given by the 0.4 dB antenna error (0.10) and the 
error due to the uncertainty in our nominal mass (0.04), 
we arrive at the result 



^rp 



= &rp ± cTrp = (1-10 ± 0.11) X 10"^ cm/s^ 



(35) 



B. RTG heat reflecting off" the spacecraft 



It has been argued that the anomalous acceleration 
seen in the Pioneer spacecraft is due to anisotropic heat 
reflection off of the back of the spacecraft high-gain an- 
tennae, the heat coming from the RTGs | 103| . Before 
launch, the four RTGs had a total thermal fuel inventory 
of 2580 W (now - 2070 W). They produced a total elec- 
trical power of 160 W (now ~ 65 W). Presently ~ 2000 
W of RTG heat must be dissipated. Only - 63 W of 
directed power could explain the anomaly. Therefore, in 
principle there is enough power to explain the anomaly 
this way. However, there are two reasons that preclude 
such a mechanism, namely: 

i) The spacecraft geometry: The RTGs are located at 
the end of booms, and rotate about the spacecraft in a 
plane that contains the approximate base of the antenna. 
From the closest axial center point of the RTGs, the an- 
tenna is seen nearly "edge on" (the longitudinal angular 
width is 24.5°). The tota l solid angle subtended is ^ 1- 
2% of 47r steradians |104| . Even though a more detailed 
calculation yields a value of 1.5% |l05j . even taking the 
higher bound of 2% means this pr opos al could provide at 
most ^ 40 W. But there is more |106| |. 

ii) The RTGs' radiation pattern: The above estimate 
was based on the assumption that the RTGs are spher- 
ical black bodies. But they are not. The main bodies 
of the RTGs are cylinders and they are grouped in two 
packages of two. Each package has the two cylinders end 
to end extending away from the antenna. Every RTG has 
six fins separated by equal angles of 60 degrees that go 
radially out from the cylinder. Presumably this results 
in a symmetrical radiation of thermal power into space. 

Thus, the fins are "edge on" to the antenna (the fins 
point perpendicular to the cylinder axes). The largest 
opening angle of the fins is seen only by the narrow-angle 
parts of the antenna's outer edges. Ignoring these edge 
effects, only ^^2.5% of the surface area of the RTGs is 
facing the antenna. This is a factor 10 less than that from 
integrating the directional intensity from a hemisphere: 
[(/''■'P'' df]cos0)/(47r)] = 1/4. So, one has only 4 W 
of directed power. This suggests a systematic bias of 
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^ 0.55 X 10 ^ cm/s^. Even adding an uncertainty of the 
same size yields a systematic for heat reflection of 

flh.r. = (-0.55 ± 0.55) X 10"* cm/s^ (36) 

But there are reasons to consider this an upper bound. 
The Pioneer SNAP 19 RTGs have larger fins than the 
earlier test models and the packages were insulated so 
that the end caps have lower temperatures. This results 
in lower radiation from the end caps than from the cylin- 
der/fins pol I2H . As a result, even though this is not 
exact, we can argue that the vast majority of the heat 
radiated by the RTGs is symmetrically directed to space 
unobscured by the antenna. Further, for this mechanism 
to work one still has to assume that the energy hitting 
the antenna i s com pletely reradiated in the direction of 
the spin axis |l06j . 

Finally, if this mechanism were the cause, ultimately 
an unambiguous decrease in the size of ap should be seen 
because the RTGs' radioactively produced radiant heat 
is decreasing. As noted previously, the heat produced is 
now about 80% of the original magnitude. In fact, one 
would similarly expect a decrease of about 0.75 x 10"* 
cm/s^ in ap over the 11.5 year Pioneer 10 data interval 
if this mechanism were the origin of ap. 

So, even though a complete thermal/physical model of 
the spacecraft might be able to ascertain if there are any 
other unsuspected heat systematics, we conclude that 
this particular mechanism does n ot pr ovide enough power 
to explain the Pioneer anomaly |107| . 

In addition to the observed constancy of the anomalous 
acceleration, any explanation involving thermal radiation 
must also discuss the absence of a disturbance to the spin 
of the spacecraft. There may be a small correlation of 
the spin angular acceleration with the anomalous linear 
acceleration. However, as described in Section IVII the 
linear acceleration is much more constant than the spin. 
This suggests that most of the linear acceleration is not 
caused by whatever disturbs the spin, thermal or not. 

However, a careful look at the Interval I results of Fig- 
ure llll shows that the nearly steady, background spin-rate 
change of about 6 x 10~^ rpm/day is slowly decreasing. 

In principle this could be caused by heat. 

The spin-rate change produced by the torque of radiant 
power directed against the rotation with a lever arm d is 



where Xz is the moment of inertia, 588.3 kg m^ We 
take a base unit of 9o for a power of one Watt and a lever 
arm of one meter. This is 

00 ^ 5.63 X 10"^^ rad/s^ = 4.65 x 10"** rpm/day = 
= 1.71 X 10"^ rpm/yr. (38) 

So, about 13 Watt-meters of directed power could cause 
the base spin-rate change. 



It turns out that such sources could, in principle, be 
available. There are 3x3 = 9 radioisotope heater units 
(RHUs) with one Watt power to heat the Thruster Clus- 
ter Assembly (TCA). (See pages 3.4-4 and 3.8-1-3.8-17 
of Ref. 4].) The units are on the edge of the antenna of 
radius 1.37 m, in the housings of the TCAs which are ap- 
proximately 180° apart from each other. At one position 
there are six RHUs and at the other position there are 
three. An additional RHU is near the sun sensor which 
is located near the second assembly. The final RHU is 
located at the magnetometer, 6.6 meters out from the 
center of the spacecraft. 

The placement gives an "ideal" rotational asymmetry 
of two Watts. But note, the real asymmetry should be 
less, since these RHUs do not radiate only in one direc- 
tion. Even one Watt unidirected at the magnetometer, is 
not enough to cause the baseline spin rate decrease. Fur- 
ther, since the base line is decreasing faster than what 
would come from the change cause by radioactive decay 
decrease, one cannot look for this effect or some com- 
plicated RTG source as the entire origin of the baseline 
change. One would suspect a very small gas leak or a 
combination of this and heat from the powered bus. (See 
Section rvTHDl ) Indeed, the factor 1/c in Eq. (|S7|) is 
a manifestation of the energy-momentum conservation 
power needed to produce 9 by heat vs. massive particles. 

But in any event, this baseline spin-rate change is not 
significantly correlated with the anomalous acceleration, 
so we do not have to pursue it further. 



C. Differential emissivity of the RTGs 

A noth er suggestion related to the RTGs is the follow- 
ing |108| |: during the early parts of the missions, there 
might have been a differential change of the radiant emis- 
sivity of the solar-pointing sides of the RTGs with re- 
spect to the deep-space facing sides. Note that, especially 
closer in the Sun, the inner sides were subjected to the 
solar wind. Contrariwise, the outer sides were sweeping 
through the solar-system dust cloud. Therefore, it can 
be argued that these two processes could have caused 
the effect. However, other information seems to make it 
difhcult for this explanation to work. 

The six fins of each RTG, designed to "provide the 
bulk of the heat rejection capacity," were fabricated of 
HM21A-T8 magnesium alloy plate The metal, after 
being specially prepared, was coated with two to three 
mils of zirconia in a sodium silicate binder to provide 
a high emissivity {^-^ 0.9) and low absorptivity 0.2). 
Depending on how symmetrically fore-and-aft they ra- 
diated, the relative fore-and-aft emissivity of the alloy 
would have had to have changed by ~ 10% to account 
for ap (see below). Given our knowledge of the solar 
wind and the interplanetary dust (see Section IXI Al) . we 
find that this amount of a radiant change would be diffi- 
cult to explain, even if it were of the right sign. (In fact, 
even the brace bars holding the RTGs were built such 



34 



that radiation is roughly fore/aft symmetric,) 

We also have "visual" evidence from the Voyager 
spacecraft. As mentioned, the Voyagers are not spin- 
stabilized. They have imaging video cameras attached 
flOQj . The cameras are mounted on a scan platform that 
is pointed u nder both celestial and inertial attitude con- 
trol modes The cameras do not have lens covers 
During the outward cruise calibrations, the cam- 
eras were sometimes pointed towards an imaging tar- 
get plate mounted at the rear of the spacecraft. But 
most often they were pointed all over the sky at specific 
star fields in support of ultraviolet spectrometer observa- 
tions. Meanwhile, the spacecraft antennae were pointed 
towards Earth. Therefore, at an angle, the lenses were 
sometimes hit by the solar wind and sometimes by the 
interplanetary dust. Even so, there was no noticeable 
deterioration of the imag es received, even when Voyager 
2 reached Neptune '112]. We infer, therefore, that this 
mechanism can not explain the Pioneer effect. 

It turned out that the greatest radiation damage oc- 
curred during the flybys. The peak Pioneer 10 radiation 
flux near Jupiter was about 10,000 times that of Earth 
for electrons (1,000 times for protons). Pioneer 11 expe- 
rienced an even higher radiation flux and also went by 
Saturn (We return to this in Section IVlTmi ) There- 
fore, if radiation damage was a problem, one should have 
seen an approximately uniform change in emissivity dur- 
ing flyby. Since the total heat flux, J^, from the RTGs was 
a constant over a flyby, there would have been a change 
in the RTG surface temperature manifested by the radi- 
ation formula !F oc eirf = e2T2 , the being the emissiv- 
ities of the fin material. There are several temperature 
sensors mounted at RTG fin bases. They measured aver- 
age temperatures of approximately 330 F, roughly 440 K. 
Therefore, a 10% change in the total average emissivity 
would have produced a temperature change of ~12.2 K 
= 22 F. Such a change would have been noticed. (Mea- 
surements would be compared from, say, 30 days before 
and after flyby to eliminate the flyby power/thermal dis- 
tortions.) Since (see below) a 10% differential fore/aft 
emissivity could cause the Pioneer effect, the lack of ob- 
servation of a 10% total average emissivity change limits 
the size of the differential emissivity systematic. 

To obtain a reasonable estimate of the uncertainty, 
consider if one side (fore or aft) of the RTGs had its 
emissivity changed by 1% with respect to the other side. 
In a simple cylindrical model of the RTGs, with 2000 
W power (here we presume only radial emission with no 
loss out the sides), the ratio of power emitted by the 
two sides would be 0.99 = 995/1005, or a differential 
emission between the half cylinders of 10 W. Therefore, 
the fore/aft asymmetry towards the normal would be 
[10 W] X [sin (I)]d4> / n w 6.37 W. If one does a more 
sophisticated fin model, with 4 of the 12 fins facing the 
normal (two flat and two at 30°), one gets a number of 
6.12 W. We take this to yield our uncertainty, 

CTd.e. = 0.85 X 10"* cm/s^ (39) 



Note that 10 (Td.e. almost equals our final ap. This is the 
origin of our previous statement that ~ 10% differential 
emissivity (in the correct direction) would be needed to 
explain ap. 

Finally, we want to comment on the significance of ra- 
dioactive decay for this mechanism. Even acknowledging 
the Interval jumps due to gas leaks (see below), we re- 
ported a one-day batch-sequential value (before system- 
atics) for ap, averaged over the entire 11.5 year interval, 
of ap = (7.77 ± 0.16) x 10"* cm/s^. From radioactive 
decay, the value of ap should have decreased by 0.75 of 
these units over 11.5 years. This is 5 times the above 
variance, which is very large with batch sequential. Even 
more stringently, this bound is good for all radioactive 
heat sources. So, what if one were to argue that emissiv- 
ity changes occurring before 1987 were the cause of the 
Pioneer effect? There still should have been a decrease 
in ap with time since then, which has not been observed. 

We will return to these points in Section IVIII Gl 



D. Non-isotropic radiative cooling of the spacecraft 

It has also been suggested that the anomalous acceler- 
ation seen in the Pioneer 10/11 spacecraft can be, "ex- 
plained, at least in p art, b y non-isotropic radiative cool- 
ing of the spacecraft |ll3j ." So, the question is, does "at 
least in part" mean this effect come s nea r to explaining 
the anomaly? We argue it does not |114| . 

Consider radiation of the main-bus electrical systems 
power from the spacecraft rear. For the Pioneers, the 
aft has a louver system, and "the louver system acts to 
control the heat rejection of the radiating platform. A 
bimetallic spring, thermally coupled radiatively to the 
platform, provides the motive force for altering the angle 
of each blade. In a closed position (below 40 F) the heat 
rejection of the platform is minimized by virtue of the 
blockage of the blades while open fin louvers provide the 
platform with a nearly unobstructed view of space 0." 

If these louvers were open (above ^ 88 F) and all the 
diminishing electrical-power heat was radiated only out 
of the louvers, this mechanism could produce a signif- 
icant effect. However, by nine AU the actuator spring 
temperature had already reached ~40 F Q . This means 
the louver doors were closed (i.e., the louver angle was 
zero) from where we obtained our data. Thus, from that 
time on of the radiation properties, the contribution of 
the thermal radiation to the Pioneer anomalous accelera- 
tion should be small. Although one might speculate that 
a louver stuck, there are 30 louvers on each craft. They 
clearly worked as designed, or else the temperature of the 
crafts' interiors would have fallen to disastrous levels. 

As shown in Figure lTCl in 1984 Pioneer 10 was at about 
33 AU and the power was about 105 W. (Always reduce 
the total power numbers by 8 W to account for the radio 
beam power.) In (1987, 1992, 1996) the spacecraft was 
at ~(41, 55, 65) AU and the power was ~(95, 82, 73) 
W. The louvers were inactive, and no decrease in ap was 



35 



seen. 

190 r- 



lao 




1000 2000 3000 4000 5000 6000 7000 BOOO 9000 

Da^s Ifom Launch 



FIG. 16: The Pioneer 10 electrical power generated at the 
RTGs as a function of time from launch to near the end of 
1994. By 1998.5, only ~68 W was generated. 



three times this error to be our systematic uncertainty 
for radiative cooling of the craft, CTr.c. = ±0.48 x 10~* 
cm/s^. 

Although doubtful, one can also speculate that some 
mechanism like this might be involved with the baseline 
spin-rate change discussed in Section FVIII Bl In 1986-7, 
Pioneer 10 power was about 97 W, decreasing at about 
2.5-3.0 W/yr. If you take a lever arm of 0.71 meters (the 
hexagonal bus size), this is more than enough to pro- 
vide the 13 W-meters necessary to produce the baseline 
spin-rate change of Figure [TTl Further for the first three 
years the decrease about matches the bus power loss rate. 
Then after the complex changes associated with the end 
of 1989 to 1990, there is a decrease in the base rate with 
a continued similar slope. 

Perhaps the "baseline" rate is indeed from the heat of 
the bus being vented to the side. But the much larger 
gas leaks would be on top of the baseline. 



In fact, during the entire 11.5 year period from 1987 to 
1998 the electrical power decreased from around 95 W to 
around 68 W, a change of 27 W. Since we already have 
noted that about ~ 65 W is needed to cause our effect, 
such a large decrease in the "source" of the acceleration 
would have been seen. But as shown in Sect ion IvTI it was 
not. Even the small differences in the three intervals are 
most likely to be from gas leaks (as will be demonstrated 

in Section rvTrmi . 

Later a double modification of this idea was given. It 
was first suggested that "most, if not all, of the unmod- 
eled acceleration" of Pioneer 10 and 11 is due to an es- 
sentially constant supply of heat coming from the central 
compartment, dire cted out the front of the craft through 
the closed louvers |ll5l |-faV However, when one studies 
the electrical power history in both parts (instruments 
and experimental) of the central compartment, there is 
no constancy of heat. (See the details in |ll6j .') Indeed 
during our data period the heat from this compartment 
decreased from about 73 W to about 57 Watts, or a fac- 
tor of 1.26. This is inconsistent with the constancy of our 
result. F urthe r, if one looks at the earlier, very roughly 
analyzed Sll'^ data in Figure [3 one sees nothing close to 
the i nternal power change of 93 to 57 W (a factor of 1.6) 
[Til . 

To address this inconsistency a second modification 
|ll5| -b.c was made. It was arbitrarily argued that 
there was an incorrect determination of the reflec- 
tion/absorption coefRcients by a large factor. But these 
coefficients are known to 5%. If they were as poorly de- 
termined as speculated, the missi on w ould have failed 
early on. (Further discussion is in |ll6l |.) 

We conclude that neither the original proposal '1135 
nor the modification .115j can explain the anomalous 
Pioneer acceleration A bound on the con- 

stancy of ap comes from first noting the 11.5 year 1- 
day batch-sequential result, sensitive to time variation: 
ap = (7.77±0.16)xl0"^cm/s^. Also given the constancy 
of the earlier imprecise date, it is conservative to take 



E. Expelled Helium produced within the RTGs 

Another possible on-board systematic is from the ex- 
pulsion of the He being created in the RTGs from the 
a-decay of ^^^Pu. To make this mechanism work, one 
would need that the He leakage from the RTGs be pref- 
erentially directed away from the Sun, with a velocity 
large enough to cause the acceleration. 

The SNAP- 19 Pioneer RTGs were designed in a such a 
way that the He pressure has not been totally contained 
within the Pioneer heat source over the life of RTGs . 
Instead, the Pioneer heat source contains a pressure re- 
lief device which allows the generated He to vent out 
of the heat source and into the thermoelectric converter. 
(The strength member and the capsule clad contain small 
holes to permit He to escape into the thermoelectric con- 
verter.) The thermoelectric converter housing-to-power 
output receptacle interface is sealed with a viton 0-ring. 
The 0-ring allows the helium gas within the converter 
to be released by permeation to the space environment 
throughout the mission life of the Pioneer RTGs. 

Information on the fuel pucks |ll8l | shows that they 
each have heights of 0.212 inches with diameters of 2.145 
inches. With 18 in each RTG and four RTGs per mission, 
this gives a total volume of fuel of about 904 cm'^. The 
fuel is PMC Pu conglomerate. The amount of ^^^Pu in 
this fuel is about 5.8 kg. With a half life of 87.74 years, 
that means the rate of He production (from Pu decay) 
is about 0.77 gm/year, assuming it all leaves the cermet. 
Taking on operational temperature on the RTG surface 
of 320 F 433 K, imphes a 3fcr/2 hehum velocity of 
1.22 km/s. (The possible energy loss coming out of the 
viton is neglected for helium.) Using this in the rocket 
equation. 



a{t) = \nM{t) 



(40) 



with our nominal Pioneer mass with half the fuel gone 
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and the assumption that the gas is all unidirected, yields a 
maximal bound on the possible acceleration of 1.16x 10^^ 
cm/s^. So, we can rule out helium permeating through 
the 0-rings as the cause of ap although it is a systematic 
to be dealt with. 

Of course, the gas is not totally unidirected. As one 
can see by looking at Figures El and III-2 of '2^ : the 
connectors with the 0-rings are on the RTG cylinder 
surfaces, on the ends of the cylinders where the fins are 
notched. They are equidistant (30 degrees) from two of 
the fins. The placement is exactly at the "rear" direction 
of the RTG cylinders, i.e., at the position closest to the 
Sun/Earth. The axis through the 0-rings is parallel to 
the spin-axis. The 0-rings, sandwiched by the receptacle 
and connector plates, "see" th e out side world through an 
angle of about 90° in latitude [iTflj . (Overhead of the O- 
rings is towards the Sun.) In longitude the 0-rings see 
the direction of the bus and space through about 90°, and 
"see" the fins through most of the rest of the longitudinal 
angle. 

If one assumes a single elastic reflection, one can es- 
timate the fraction of the bias away from the Sun. (In- 
deed, multiple and back reflections will produce an even 
greater bias. Therefore, we feel this approximation is 
justified.) This estimate is (3/4) sin 30° times the av- 
erage of the heat momentum component parallel to the 
shortest distance to the RTG fin. Using this, we find the 
bias would be 0.31 x 10^^ cm/s^. This bias effectively 
increases the value of our solution for a p , which we hesi- 
tate to accept given all the true complications of the real 
system. Therefore we take the systematic expulsion to 
be ane = (0.15 ± 0.16) x 10"^ cm/s^. 



F. Propulsive mass expulsion due to gas leakage 

The eff'ect of propulsive mass expulsion due to gas leak- 
age has to be assessed. Although this effect is largely un- 
predictable, many spacecraft have experienced gas leaks 
producing accelerations on the order of 10~^ cm/s^. [The 
reader will recall the even higher figure for Ulysses found 
in Section FV D 21 ] As noted previously, gas leaks gener- 
ally behave differently after each maneuver. The leakage 
often decreases with time and becomes negligibly small. 

Gas leaks can originate from Pioneer's propulsion sys- 
tem, which is used for mid-course trajectory maneu- 
vers, for spinning-up or -down the spacecraft, and for 
orientation of the spinning spacecraft. The Pioneers 
are equipped with three pairs of hydrazine thrusters 
which are mounted on the circumference of the Earth- 
pointing high gain antenna. Each pair of thrusters forms 
a Thruster Cluster Assembly (TCA) with two nozzles 
aligned in opposition to each other. For attitude control, 
two pairs of thrusters can be fired forward or aft and are 
used to precess the spinning antenna (See Section HlBI l 
The other pair of thrusters is aligned parallel to the rim 
of the antenna with nozzles oriented in co- and contra- 
rotation directions for spin/despin maneuvers. 



During both observing intervals for the two Pioneers, 
there were no trajectory or spin/despin maneuvers. So, 
in this analysis we are mainly concerned with preces- 
sion (i.e., orientation or attitude control) maneuvers only. 
(See Section FlI Bl ) Since the valve seals in the thrusters 
can never be perfect, one can ask if the leakages through 
the hydrazine thrusters could be the cause of the anoma- 
lous acceleration, ap. 

However, when we investigate the total computational 
accuracy of our solution in Section IIX Bl we will show 
that the currently implemented models of propulsion ma- 
neuvers may be responsible for an uncertainty in ap only 
at the level of ±0.01 x 10~^ cm/s^. Therefore, the ma- 
neuvers themselves are the main contributors neither to 
the total error budget nor to the gas leak uncertainty, as 
we now detail 

The serious uncertainty comes from the possibility of 
undetected gas leaks. We will address this issue in some 
detail. First consider the possible action of gas leaks 
originating from the spin/despin TCA. Each nozzle from 
this pair of thrusters is subject to a certain amount of 
gas leakage. But only a differential leakage from the two 
nozzles would produce an observable effec t cau sing the 
spacecraft to either spin-down or spin- up |120| . So, to 
obtain a gas leak uncertainty (and we emphasize "uncer- 
tainty" vs. "error" because we have no other evidence) 
let us ask how large a differential force is needed to cause 
the spin-down or spin-up effects observed? 

Using the moment of inertia about the spin axis, Iz =~ 
588.3 kg-m^ (T^J, and the antenna radius, TZ = 1.37 m, as 
the lever arm, one can calculate that the differential force 
needed to torque the spin-rate change, in Intervals 
i =1,11,111 is 

X 6- 

Fg^ = = (2.57, 12.24, 1.03) X 10-3 dynes. (41) 

It is possible that a similar mechanism of undetected 
gas leakage could be responsible for the net differential 
force acting in the direction along the line of sight. In 
other words, what if there were some undetected gas leak- 
age from the thrusters oriented along the spin axis of the 
spacecraft that is causing ap? How large would this have 
to be? A force of (Af = 241 kg) 

= Af ap = 21.11 X IQ-^ dynes (42) 

would be needed to produce our final unbiased value of 
ap. (See Section |3 That is, one would need even more 
force than is needed to produce the anomalously high 
rotational gas leak of Interval II. Furthermore, the dif- 
ferential leakage to produce this ap would have had to 
have been constant over many years and in the same di- 
rection for both spacecraft, without being detected as 
a spin-rate change. That is possible, but certainly not 
demonstrated. Furthermore if the gas leaks hypothesis 
were true, one would expect to see a dramatic difference 
in ap during the three Intervals of Pioneer 10 data. In- 
stead an almost 500 % spin-down rate change between 
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Intervals I and II resulted only in a less than 8% change 
in ap. 

Given the small amount of information, we propose 
to conservatively take as our gas leak uncertainties the 
acceleration values that would be produced by differential 
forces equal to 

F,^(,)g.i. 0, ±V2F^^ = (43) 
= ( ± 3.64, ± 17.31, ± 1.46) x 10~^ dynes. 

The argument for this is that, in the root sum of squares 
sense, one is accounting for the differential leakages from 
the two pairs of thrusters with their nozzles oriented 
along the line of sight direction. This directly trans- 
lates into the acceleration errors introduced by the leak- 
age during the three intervals of Pioneer 10 data, 

CT(ap(.Og.i.) = ±^ap (Og.i.A^ = (44) 
= (±1.51, ±7.18, ±0.61) X 10"* cm/s^ 

Assuming that these errors are uncorrelated and are nor- 
mally distributed around zero mean, we find the gas leak 
uncertainty for the entire Pioneer 10 data span to be 

o-g.i. = ±0.56 X 10"** cm/s^ (45) 

This is one of our largest uncertainties. 

The data set from Pioneer 11 is over a much smaller 
time span, taken when Pioneer 11 was much closer to 
the Sun (off the plane of the ecliptic) , and during a max- 
imum of solar activity. For Pioneer 11 the main effects 
of gas leaks occurred at the maneuvers, when there were 
impulsive lowerings of the spin-down rate. These dom- 
inated the over-all spin rate change of 6*11 = —0.0234 
rpm/yr. (See Figure IT^) But in between maneuvers the 
spin rate was actually increasing. One can argue that 
this explains the higher value for ap^n) in Table as 
compared to ap(^iQy Unfortunately, one has no a priori 
way of predicting the effect here. We do not know that 
the same specific gas leak mechanism applied here as did 
in the case of Pioneer 10 and there is no well-defined in- 
terval set as there is for Pioneer 10. Therefore, although 
we feel this "spin up" may be part of the explanation 
of the higher value of ap for Pioneer 11, we leave the 
different numbers as a separate systematic for the next 
subsection. 

At this point, we must conclude that the gas leak mech- 
anism for explaining the anomalous acceleration seems 
very unlikely, because it is hard to understand why it 
would affect Pioneer 10 and 11 at the same level (given 
that both spacecraft had different quality of propulsion 
systems, see Section FlI Bll . One also expects a gas leak 
would obey the rules of a Poisson distribution. That 
clearly is not true. Instead, our analyses of different data 
sets indicate that ap behaves as a constant bias rather 
than as a random variable. (This is clearly seen in the 
time history of ap obtained with batch-sequential esti- 
mation.) 



G. Variation between determinations from the two 
spacecraft 

Finally there is the important point that we have two 
"experimental" results from the two spacecraft, given in 
Eqs. and (j^SJ: 7.84 and 8.55, respectively, in units 
of 10"* cm/s^. If the Pioneer effect is real, and not a sys- 
tematic, these numbers should be approximately equal. 

The first number, 7.84, is for Pioneer 10. In Section 
IVIDI we obtained this number by correlating the values 
of ap in the three data Intervals with the different spin- 
down rates in these Intervals. The weighted correlation 
between a shift in ap and the spin-down rate is kq = 
(30.7±0.6) cm. (We argued in the previous Section lVlII 11 
that this correlation is the manifestation of the rotational 
gas leak systematic.) Therefore, this number represents 
the entire 11.5 year data arc of Pioneer 10. Similarly, 
Pioneer ll's number, 8.55, represents a 3| year data arc. 

Even though the Pioneer 11 number may be less re- 
liable since the craft was so much closer to the Sun, 
we calculate the time-weighted average of the exper- 
imental results from the two craft: [(11.5)(7.84) + 
(3.75)(8.55)]/(15.25) = 8.01 in units of 10"* cm/s^. This 
implies a bias of b2_cza.1t = +0.17 x 10~* cm/s^ with re- 
spect to the Pioneer 10 experimental result ap(exper)- We 
also take this number to be our two spacecraft uncer- 
tainty. This means 

a2-craft = ^2-craft ± C''2_craft = 

= (0.17 ±0.17) X 10"* cm/s^ (46) 



The difference between the two craft could be due to 
different gas leakage. But it also could be due to heat 
emitted from the RTGs. In particular, the two sets of 
RTGs have had different histories and so might have dif- 
ferent emissivities. Pioneer 11 spent more time in the 
inner solar system (absorbing radiation). Pioneer 10 has 
swept out more dust in deep space. Further, Pioneer 11 
experienced about twice as much Jupiter/Saturn radia- 
tion as Pioneer 10. 

Further, note that [ap^exper) ~ °^p(°ex°per)] ^'^'^ uncer- 
tainty from differential emissivity of the RTGs, cr^.e., are 
of the same size: 0.71 and 0.85 xlO"* cm/s^. It could 
therefore be argued that Pioneer ll's offset from Pioneer 
10 comes from Pioneer 11 having obtained twice as large 
a differential emissivity bias as Pioneer 10. Then our fi- 
nal value of ap, given in Section ^ would be reduced 
by about 0.7 of our units since dd.e. would have become 
mainly a negative bias, b^.e.- This would make the final 
number closer to 8 x 10~* cm/s^. Because this model 
and our final number are consistent, we present this ob- 
servation only for completeness and as a possible reason 
for the different results of the two spacecraft. 
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IX. COMPUTATIONAL SYSTEMATICS 

Given the very large number of observations for the 
same spacecraft, the error contribution from observa- 
tional noise is very small and not a meaningful measure of 
uncertainty. It is therefore necessary to consider several 
other effects in order to assign realistic errors. Our first 
consideration is the statistical and numerical stability of 
of the calculations. We then go on to the cumulative 
influence of all modeling errors and editing decisions. Fi- 
nally we discuss the reasons for and significance of the 
annual term. 

Besides the factors mentioned above, we will discuss in 
this section errors that may be attributed to the specific 
hardware used to run the orbit determination computer 
codes, together with computational algorithms and sta- 
tistical methods used to derive the solution. 



A. Numerical stability of least-squares estimation 

Having presented estimated solutions along with their 
formal statistics, we should now attempt to characterize 
the true accuracy of these results. Of course, the signifi- 
cance of the results must be assessed on the basis of the 
expected measurement errors. These expected errors are 
used to weight a least-squares adjustment to parameters 
which describe the theoretical model. [Examination of 
experimental systematics from sources both external to 
and also internal to the spacecraft was covered in Sections 
IViilVllll l 

First we look at the numerical stability of the least 
squares estimation algorithm and the derived solution. 
The leading computational error source turns out to be 
subtraction of similar numbers. Due to the nature of 
floating point arithmetic, two numbers with high order 
digits the same are subtracted one from the other re- 
sults in the low order digits being lost. This situation 
occurs with time tags on the data. Time tags are refer- 
enced to some epoch, such as say 1 January 1 1950 which 
is used by CHASMP. As more than one billion seconds 
have passed since 1950, time tags on the Doppler data 
have a start and end time that have five or six common 
leading digits. Doppler signal is computed by a differ- 
enced range formulation (see Section Fill Bp . This noise 
in the time tags causes noise in the computed Doppler 
at the 0.0006 Hz level for both Pioneers. This noise can 
be reduced by shifting the reference epoch closer to the 
data or increasing the word length of the computation, 
however, it is not a significant error source for this anal- 
ysis. 

In order to guard against possible computer com- 
piler and/or hardware errors we ran orbit determina- 
tion programs on different computer platforms. JPL's 
ODP resides on an HP workstation. The Aerospace 
Corporation ran the analysis on three different com- 
puter architectures: (i) Aerospace's DEC 64-bit RISC 
architecture workstation (Alphastation 500/266), (ii) 



Aerospace's DEC 32-bit CISC architecture workstation 
(VAX 4000/60), and (iii) Pentium Pro PC. Comparisons 
of computations performed for CHASMP in the three 
machine show consistency to 15 digits which is just suffi- 
cient to represent the data. While this comparison does 
not eliminate the possibility of systematic errors that are 
common to both systems, it does test the numerical sta- 
bility of the analysis on three very different computer 
architectures. 

The results of the individual programs were given 
in Sections Ivknd IVII In a test we took the JPL re- 
sults for a batch-sequential Sigma run with 50-day av- 
erages of the anomalous acceleration of Pioneer 10, ap. 
The data interval was from January 1987 to July 1998. 
We compared this to an Aerospace determination using 
CHASMP, where the was split into 200 day intervals, 
over a shorter data interval ending in 1994. As seen in 
Figure ITU the results basically agree. 

Given the excellent agreement in these implementa- 
tions of the modeling software, we conclude that differ- 
ences in analyst choices (parameterization of clocks, data 
editing, modeling options, etc.) give rise to coordinate 
discrepancies only at the level of 0.3 cm. This number 
corresponds to an uncertainty in estimating the anoma- 
lous acceleration on the order of 8 x 10^^^ cm/s^. 

But there is a slightly larger error to contend with. In 
principle the STRIPPER can give output to 16 significant 
figures. From the beginning the output was-rounded off 
to 15 and later to 14 significant figures. When Block 
5 came on near the beginning of 1995, the output was 
rounded off to 13 significant figures. Since the Doppler 
residuals are 1.12 mm/s this last truncation means an 
error of order 0.01 mm/s. If we divide this number by 2 
for an average round off, this translates to ±0.04 x 10^^ 
cm/s^. The roundoff occurred in approximately all the 
data we added for this paper. This is the cleanest 1/3 
of the Pioneer 10 data. Considering this we take the 
uncertainty to be 

CTnum ± 0.02 X 10"^ cm/s^ (47) 

It needs to be stressed that such tests examine only 
the accuracy of implementing a given set of model codes, 
without consideration of the inherent accuracy of the 
models themselves. Numerous external tests, which we 
have been discussing in the previous three sections, are 
possible for assessing the accuracy of the solutions. Com- 
parisons between the two software packages enabled us 
to evaluate the implementations of the theoretical mod- 
els within a particular software. Likewise, the results of 
independent radio tracking observations obtained for the 
different spacecraft and analysis programs have enabled 
us to compare our results to infer realistic error levels 
from differences in data sets and analysis methods. Our 
analysis of the Galileo and Ulysses missions (reported in 
Sections IV CI and IV D|) was done partially for this pur- 
pose. 
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B. Accuracy of consistency /model tests 

a. Consistency of solutions: A code that models 
the motion of solar system bodies and spacecraft in- 
cludes numerous lengthy calculations. Therefore, the 
software used to obtain solutions from the Doppler data 
is, of necessity, very complex. To guard against poten- 
tial errors in the implementation of these models, we 
used two software packages; JPL's ODP/ Sigma model- 
ing software 41, 54| and The Aerospace Corporation's 
POEAS/CHASMP software package Izllzl. The differ- 
ences between the JPL and Aerospace orbit determina- 
tion program results are now examined. 

As discussed in Section IIV Fl in estimating parame- 
ters the CHASMP code uses a standard variation of pa- 
rameters method whereas ODP uses the Cowell method 
to integrate the equations of motion and the variational 
equations. In other words, CHASMP integrates six first- 
order differential equation, using the Adams-Moulton 
predictor-corrector method in the orbital elements. Con- 
trariwise, ODP integrates three second-order differential 
equations for the accelerations using the Gauss-Jackson 
meth od. (For more details on these methods see Ref. 

[ml.) 

As seen in our results of Sections IVbnd IVTl agreement 
was good; especially considering that each program uses 
independent methods, models, and constants. Internal 
consistency tests indicate that a solution is consistent at 
the level of one part in 10^^. This implies an acceleration 
error on the order of no more then one part in 10^ in ap. 

b. Earth orientation parameters: In order to check 
for possible problems with Earth orientation, CHASMP 
was modified to accept Earth orientation information 
from three different sources. (1) JPL's STOIC pro- 
gram that outputs UTIR-UTC, (2) JPL's Earth Orienta- 
tion Parameter files (UTl-UTC), and (3) The International 
Earth Rotation Service's Earth Orientation Parameter 
file (UTl-UTC). We found that all three sources gave vir- 
tually identical results and changed the value of ap only 
in the 4th digit p^ . 

c. Planetary ephemeris: Another possible source of 
problems is the planetary ephemeris. To explore this a 
fit was first done with CHASMP that used DE200. The 
solution of that fit was then used in a fit where DE405 
was substituted for DE200. The result produced a small 
annual signature before the fit. After the fit, the maneu- 
ver solutions changed a small amount (less then 10%) 
but the value of the anomalous acceleration remained 
the same to seven digits. The post-fit residuals to DE405 
were virtually unchanged from those using DE200. This 
showed that the anomalous acceleration was unaffected 
by changes in the planetary ephemeris. 

This is pertinent to note for the following subsection. 
To reemphasize the above, a small "annual term" can be 
introduced by changing the planetary ephemerides. This 
annual term can then be totally taken up by changing 
the maneuver estimations. Therefore, in principle, any 
possible mismodeling in the planetary ephemeris could be 



at least partially masked by the maneuver estimations. 

d. Differences in the codes' model implementations: 
The impact of an analyst's choices is difficult to address, 
largely because of the time and expense required to pro- 
cess a large data set using complex models. This is espe- 
cially important when it comes to data editing. It should 
be understood that small differences are to be expected 
as models differ in levels of detail and accuracy. The an- 
alysts' methods, experience, and judgment differ. The 
independence of the analysis of JPL and Aerospace has 
been consistently and strictly maintained in order to pro- 
vide confidence on the validity of the analyses. Acknowl- 
edging such difficulties, we still feel that using the very 
limited tests given above is preferable to an implicit as- 
sumption that all analysts' choices were optimally made. 

Another source for differences in the results presented 
in Table is the two codes' modeling of spacecraft re- 
orientation maneuvers. ODP uses a model that solves 
for the resulted change in the Doppler observable Aw 
(instantaneous burn model). This is a more convenient 
model for Doppler velocity measurements. CHASMP 
models the change in acceleration, solves for Aa (finite 
burn model), and only then produces a solution for Av. 
Historically, this was done in order to incorporate range 
observations (for Galileo and Ulysses) into the analysis. 

Our best handle on this is the no-corona results, espe- 
cially given that the two critical Pioneer 10 Interval III 
results differed by very httle, 0.02 x lO^'* cm/s^. This 
data is least affected by maneuver modeling, data edit- 
ing, corona modeling, and spin calibration. Contrari- 
wise, for the other data, the differences were larger. The 
Pioneer Interval I and II results and the Pioneer 11 re- 
sults differed, respectively, by (0.21, 0.23, 0.25) in units 
of 10~^ cm/s^. In these intervals models of maneuvers 
and data editing were crucial. Assuming that these er- 
rors are uncorrelated, we compute their combined effect 
on anomalous acceleration ap as 

^consist /model = ±0.13 X 10~^ Cm/s^. (48) 

e. Mismodeling of maneuvers: A small contribution 
to the error comes from a possible mismodeling of the 
propulsion maneuvers. In Section flV El we found that for 
a typical maneuver the standard error in the residuals is 
ctq ~ 0.095 mm/s. 

Then we would expect that in the period between two 
maneuvers, which on average is r = 11.5/28 year, the ef- 
fect of the mismodeling would produce a contribution to 
the acceleration solution with a magnitude on the order 
of ^fljian = (Jo/t = 0.07 X 10~® cm/s^. Now let us assume 
that the errors in the Pioneer Doppler residuals are nor- 
mally distributed around zero mean with the standard 
deviation of Sa^^an that constitute a single measurement 
accuracy. Then, since there are N — 28 maneuvers in the 
data set, the total error due to maneuver mismodeling is 

a^^ = ^ = 0.01 X 10-« cm/sl (49) 
V N 
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/. Mismodeling of the solar corona: Finally, recall 
that our number for mismodeling of the solar corona, 
±0.02 X 10~^ cm/s^, was already explained in Section 

ivirn 



C. Apparent annual/diurnal periodicities in the 
solution 

In Ref. Ilj we reported, in addition to the constant 
anomalous acceleration term, a possible annual sinusoid. 
If approximated by a simple sine wave, the amplitude 
of this oscillatory term is about 1.6 x 10~^ cm/s^. The 
integral of a sine wave in the acceleration, ap, with an- 
gular velocity lo and amplitude yields the following 
first-order Doppler amplitude in two-way fractional fre- 
quency: 

At/ 2Aq 



c OJ 



(50) 



The resulting Doppler amplitude for the annual angular 
velocity 2 x 10~'^ rad/s is Aj//j/ = 5.3 x 10"^^ At 
the Pioneer downlink S-band carrier frequency of ^ 2.29 
GHz, the corresponding Doppler amplitude is 0.012 Hz 
(i.e. 0.795 mm/s). 

This term was first seen in ODP using the BSF method. 
As we discussed in Section lTV Gl treating ap as a stochas- 
tic parameter in JPL's batch-sequential analysis allows 
one to search for a possible temporal variation in this 
parameter. Moreover, when many short interval times 
were used with least-squares CHASMP, the effect was 
also observed. (See Figure [HI in Section IVll ') 

The residuals obtained from both programs are of the 
same magnitude. In particular, the Doppler residuals 
are distributed about zero Doppler velocity with a sys- 
tematic variation ~ 3.0 mm/s on a time scale of '--^ 3 
months. More precisely, the least-squares estimation 
residuals from both ODP / Sigma and CHASMP are dis- 
tributed well within a half-width taken to be 0.012 Hz. 
(See, for example, Figure|51) Even the general structures 
of the two sets of residuals are similar. The fact that 
both programs independently were able to produce simi- 
lar post-fit residuals gives us confidence in the solutions. 

With this confidence, we next looked in greater de- 
tail at the acceleration residuals from solutions for ap. 
Consider Figure ITTI which shows the ap residuals from 
a value for ap of (7.77 ± 0.16) x 10"^ cm/s^. The data 
was processed using ODP / Sigma with a batch-sequential 
filter and smoothing algorithm. The solution for ap 
was obtained using 1-day batch sizes. Also shown are 
the maneuver times. At early times the annual term is 
largest. During Interval II, the interval of the large spin- 
rate change anomaly, coherent oscillation is lost. During 
Interval III the oscillation is smaller and begins to die 
out. 

In attempts to understand the nature of this annual 
term, we first examined a number of possible sources, in- 
cluding effects introduced by imprecise modeling of ma- 
neuvers, the solar corona, and the Earth's troposphere. 




FIG. 17: ODP 1-day batch-sequential acceleration residuals 
using the entire Pioneer 10 data set. Maneuver times are 
indicated by the vertical dashed lines. 



We also looked at the influence of the data editing strate- 
gies that were used. We concluded that these effects 
could not account for the annual term. 

Then, given that the effect is particularly large in the 
out-of-the-ecliptic voyage of Pioneer 11 J^, we focused 
on the possibility that inaccuracies in solar system mod- 
eling are the cause of the annual term in the Pioneer so- 
lutions. In particular, we looked at the modeling of the 
Earth orbital orientation and the accuracy of the plane- 
tary ephemeris. 

g. Earth's orientation: We specifically modeled the 
Earth orbital elements Ap and Aq as stochastic parame- 
ters. {Ap and Aq are two o f the Set HI elements defined 
by Brouwer and Clemence |l23l |.') Sigma was applied to 
the entire Pioneer 10 data set with ap, Ap, and Aq deter- 
mined as stochastic parameters sampled at an interval of 
five days and exponentially correlated with a correlation 
time of 200 days. Each interval was fit independently, 
but with information on the spacecraft state (position 
and velocity) carried forward from one interval to the 
next. Various correlation times, 0-day, 30-day, 200-day, 
and 400-day, were investigated. The a priori error and 
process noise on Ap and Aq were set equal to 0, 5, and 
10 /Ltrad in separate runs, but only the 10 /ivad case re- 
moved the annual term. This value is at least three orders 
of magnitude too large a deviation when compared to 
the present accuracy of the Earth orbital elements. It is 
most unlikely that such a deviation is causing the annual 
term. Furthermore, changing to the latest set of EOP 
has very little effect on the residuals. [We also looked 
at variations of the other four Set HI orbital elements, 
essentially defining the Earth's orbital shape, size, and 
longitudinal phase angle. They had little or no effect on 
the annual term.] 

h. Solar system modeling: We concentrated on In- 
terval III, where the spin anomaly is at a minimum and 
where ap is presumably best determined. Further, this 
data was partially taken after the DSN's Block 5 hard- 
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ware implementation from September 1994 to August 

1995. As a result of this implementation the data is less 
noisy than before. Over Interval III the annual term is 
roughly in the form of a sine wave. (In fact, the mod- 
eling error is not strictly a sine wave. But it is close 
enough to a sine wave for purposes of our error analy- 
sis.) The peaks of the sinusoid are centered on conjunc- 
tion, where the Doppler noise is at a maximum. Look- 
ing at a CHASMP set of residuals for Interval III, we 
found a 4-parameter, nonlinear, weighted, least-squares 
fit to an annual sine wave with the parameters amplitude 
v^,t. = (0.1053 ± 0.0107) mm/s, phase (-5.3° ± 7.2°), 
angular velocity Wa.t = (0.0177 ± 0.0001) rad/day, and 
bias (0.0720 ± 0.0082) mm/s. The weights eliminate 
data taken inside of solar quadrature, and also account 
for different Doppler integration times Tc according to 
a = (0.765 mm/s) [(60 s)/Tc]^/^. This rule yields post-fit 
weighted RMS residuals of 0.1 mm/s. 

The amplitude, Ua.t., and angular velocity, Wa.t., of the 
annual term results in a small acceleration amplitude of 
fla.t. = Wa.t.Wa.t. = (0.215 ± 0.022) X 10"^ cm/s^. We will 
argue below that the cause is most likely due to errors in 
the navigation programs' determinations of the direction 
of the spacecraft's orbital inclination to the ecliptic. 

A similar troubling modeling error exists on a much 
shorter time scale that is most likely an error in the space- 
craft's orbital inclination to the Earth's equator. We 
looked at CHASMP acceleration residuals over a limited 
data interval, from 23 November 1996 to 23 December 

1996, centered on opposition where the data is least af- 
fected by solar plasma. As seen in Figure [TBI there is a 
significant diurnal term in the Doppler residuals, with pe- 
riod approximately equal to the Earth's sidereal rotation 
period (23'^56'"04".0989 mean solar time). 
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FIG. 18; CHASMP acceleration residuals from 23 November 
1996 to 23 December 1996. A clear modeling error is repre- 
sented by the solid diurnal curve. (An annual term maximum 
is also seen as a background.) 

After the removal of this diurnal term, the RMS 



Doppler residuals are reduced to amplitude 0.054 mm/s 
for Tc = 660 s {a^/iy = 2.9 x lO^^^ = 1000 s). The 

amplitude of the diurnal oscillation in the fundamental 
Doppler observable, v^.t. , is comparable to that in the 
annual oscillation, Wa.t., but the angular velocity, Wd.t., 
is much larger than Wa.t.- This means the magnitude 
of the apparent angular acceleration, a^x. = v^x.^ix. = 
(100.1 ± 7.9) X 10^^ cm/s^, is large compared to ap. Be- 
cause of the short integration times, Tc = 660 s, and long 
observing intervals, T 1 yr, the high frequency, diur- 
nal, oscillation signal averages out to less than 0.03 x 10~* 
cm/s^ over a year. This intuitively helps to explain why 
the apparently noisy acceleration residuals still yield a 
precise value of ap. 

Further, all the residuals from CHASMP and 
ODP/ Sigma are essentially the same. Since ODP and 
CHASMP both use the same Earth ephemeris and the 
same Earth orientation models, this is not surprising. 
This is another check that neither program introduces 
serious modeling errors of its own making. 

Due to the long distances from the Sun, the spin- 
stabilized attitude control, the long continuous Doppler 
data history, and the fact that the spacecraft communica- 
tion systems utilize coherent radio-tracking, the Pioneers 
allow for a very sensitive and precise positioning on the 
sky. For some cases, the Pioneer 10 coherent Doppler 
data provides accuracy which is even better than that 
achieved with VLBI observing natural sources. In sum- 
mary, the Pioneers are simply much more sensitive de- 
tectors of a number of solar system modeling errors than 
other spacecraft. 

The annual and diurnal terms are very likely different 
manifestations of the same modeling problem. The mag- 
nitude of the Pioneer 10 post-fit weighted RMS residuals 
of « 0.1 mm/s, implies that the spacecraft angular posi- 
tion on the sky is known to < 1.0 milliarcseconds (mas). 
(Pioneer 11, with 0.18 mm/s, yields the result « 1.75 
mas.) At their great distances, the trajectories of the 
Pioneers are not gravitationally affected by the Earth. 
(The round-trip light time is now 24 hours for Pioneer 
10.) This suggests that the sources of the annual and 
diurnal terms are both Earth related. 

Such a modeling problem arises when there are er- 
rors in any of the parameters of the spacecraft orien- 
tation with respect to the chosen reference frame. Be- 
cause of these errors, the system of equations that de- 
scribes the spacecraft's motion in this reference frame 
is under-determined and its solution requires non-linear 
estimation techniques. In addition, the whole estima- 
tion process is subject to Kalman filtering and smooth- 
ing methods. Therefore, if there are modeling errors in 
the Earth's ephemeris, the orientation of the Earth's spin 
axis (precession and nutation), or in the station coordi- 
nates (polar motion and length of day variations), the 
least-squares process (which determines best-fit values of 
the three direction cosines) will leave small diurnal and 
annual components in the Doppler residuals, like those 
seen in Figures [171181 
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Orbit determination programs are particularly sensi- 
tive to an error in a poorly observed direction |124j . If 
not corrected for, such an error could in principle signif- 
icantly affect the overall navigational accuracy. In the 
case of the Pioneer spacecraft, navigation was performed 
using only Doppler tracking, or line-of-sight observations. 
The other directions, perpendicular to the line-of-sight or 
in the plane of the sky, are poorly constrained by the data 
available. At present, it is infeasible to precisely parame- 
terize the systematic errors with a physical model. That 
would have allowed one to reduce the errors to a level 
below those from the best available ephemeris and Earth 
orientation models. A local empirical parameterization is 
possible, but not a parameterization over many months. 

We conclude that for both Pioneer 10 and 11, there are 
small periodic errors in solar system modeling that are 
largely masked by maneuvers and by the overall plasma 
noise. But because these sinusoids are essentially uncor- 
related with the constant ap, they do not present im- 
portant sources of systematic error. The characteristic 
signature of ap is a line ar d rift in the Doppler, not an- 
nual/diurnal signatures |125| |. 

i. Annual/diurnal mismodeling uncertainty: We 
now estimate the annual term contribution to the error 
budget for ap. First observe that the standard errors for 
radial velocity, Vr, and acceleration, a^, are essentially 
what one would expect for a linear regression. The 
caveat is that they are scaled by the root sum of squares 
(RSS) of the Doppler error and unmodeled sinusoidal 
errors, rather than just the Doppler error. Further, 
because the error is systematic, it is unrealistic to 
assume that the errors for Vr and can be reduced by 
a factor l/-\/iV, where N is the number of data points. 
Instead, averaging their correlation matrix over the data 
interval, T, results in the estimated systematic error of 

c^l = ^ <^l. = ^ (4 + <... + <.,. ) • (51) 

(Tt = 0.1 mm/s is the Doppler error averaged over T (not 
the standard error on a single Doppler measurement). 

4 and a^^ ^ are equal to the amplitudes of correspond- 
ing unmodeled annual and diurnal sine waves divided by 
V2. The resulting RSS error in radial velocity determi- 
nation is about cr„^ = (cf^+cr^^j +^^dt )^^^ ~ ^-^"^ mm/s 
for both Pioneer 10 and 11. Our four interval values of 
ap were determined over time intervals of longer than a 
year. At the same time, to detect an annual signature in 
the residuals, one needs at least half of the Earth's orbit 
complete. Therefore, with T = 1/2 yr, Eq. H51|l results 
in an acceleration error of 

We use this number for the systematic error from the 
annual/diurnal term. 



X. ERROR BUDGET AND FINAL RESULT 

It is important to realize that our experimental observ- 
able is a Doppler frequency shift, i.e., Ai/(t). [See Figure 
Inland Eq. ifT^ .] In actual fact it is a cycle count. We 
interpret this as an apparent acceleration experienced by 
the spacecraft. However, it is possible that the Pioneer 
effect is not due to a real acceleration. (See Section IxJ ^ 
Therefore, the question arises "In what units should we 
report our errors?" The best choice is not clear at this 
point. For reasons of clarity we chose units of accelera- 
tion. 

The tests documented in the preceding sections have 
considered various potential sources of systematic error. 
The results of these tests are summarized in Table El 
which serves as a systematic "error budget." This bud- 
get is useful both for evaluating the accuracy of our so- 
lution for ap and also for guiding possible future efforts 
with other spacecraft. In our case it actually is hard to 
totally distinguish "experimental" error from "system- 
atic error." (What should a drift in the atomic clocks 
be called?) Further, there is the intractable mathemati- 
cal problem of how to handle combined experimental and 
systematic errors. In the end we have decided to treat 
them all in a least squares uncorrelated manner. 

The results of our analyses are summarized in Table ITU 
There are two columns of results. The first gives a bias, 
hp, and the second gives an uncertainty, zbtxp. The con- 
stituents of the error budget are listed separately in three 
different categories: 1) systematics generated external to 
the spacecraft; 2) on-board generated systematics, and 
3) computational systematics. Our final result then will 
become some average 

ap = ap(exper) ±crp, (53) 

where, from Eq. ap(exper) = (7.84 ± 0.01) x 10"^ 

cm/s^. 

The least significant factors of our error budget are in 
the first group of effects, those external to the spacecraft. 
From the table one sees that some are near the limit of 
contributing. But in totality, they are insignificant. 

As was expected, the on-board generated systematics 
are the largest contributors to our total error budget. 
All the important constituents are listed in the second 
group of effects in Table ITU Among these effects, the ra- 
dio beam reaction force produces the largest bias to our 
result, 1.10 X 10~* cm/s^. It makes the Pioneer effect 
larger. The largest bias/uncertainty is from RTG heat 
reflecting off the spacecraft. We argued for an effect as 
large as (—0.55 ±0.55) x 10"^ cm/s^. Large uncertainties 
also come from differential emissivity of the RTGs, ra- 
diative cooling, and gas leaks, ±0.85, ±0.48, and ±0.56, 
respectively, xlO~^ cm/s^. The computational system- 
atics are listed in the third group of Table E] 

Therefore, our final value for ap is 

ap = (8.74 ± 1.33) x 10"^ cm/s^ 

~ (8.7 ± 1.3) X 10"^ cm/s^ (54) 
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TABLE II; Error Budget: A Summary of Biases and Uncertainties. 



Item Description of error budget constituents 



Bias 

10"** 



Uncertainty 
i/s^ IQ-** cm/s^ 



1 Systematics generated external to the spacecraft: 

a) Solar radiation pressure and mass 

b) Solar wind 

c) Solar corona 

d) Electro-magnetic Lorentz forces 

e) Influence of the Kuiper belt's gravity 

f) Influence of the Earth orientation 

g) Mechanical and phase stability of DSN antennae 

h) Phase stability and clocks 

i) DSN station location 

j) Troposphere and ionosphere 

2 On-board generated systematics: 

a) Radio beam reaction force 

b) RTG heat reflected off the craft 

c) Differential emissivity of the RTGs 

d) Non-isotropic radiative cooling of the spacecraft 

e) Expelled Helium produced within the RTGs 

f ) Gas leakage 

g) Variation between spacecraft determinations 

3 Computational systematics: 

a) Numerical stability of least-squares estimation 

b) Accuracy of consistency/model tests 

c) Mismodeling of maneuvers 

d) Mismodeling of the solar corona 

e) Annual/diurnal terms 



+0.03 



+ 1.10 
-0.55 



+0.15 



+0.17 



±0.01 
± < 10"^ 
±0.02 
± < 10-* 
±0.03 
±0.001 
± < 0.001 
± < 0.001 
± < 10-^ 
± < 0.001 



±0.11 
±0.55 
±0.85 
±0.48 
±0.16 
±0.56 
±0.17 



±0.02 
±0.13 
±0.01 
±0.02 
±0.32 



Estimate of total bias/error 



+0.90 



±1.33 



The effect is clearly significant and remains to be ex- 
plained. 



XI. POSSIBLE PHYSICAL ORIGINS OF THE 
SIGNAL 

A. A new manifestation of known physics? 

With the anomaly still not accounted for, possible ef- 
fects from applications of known phy sics have been ad- 
vanced. In particular, Crawford |l26l | suggested a novel 
new effect: a gravitational frequency shift of the radio 
signals that is proportional to the distance to the space- 
craft and the density of dust in the intermediate medium. 
In particular, he has argued that the gravitational inter- 
action of the S-band radio signals with the interplanetary 
dust may be responsible for producing an anomalous ac- 
celeration similar to that seen by the Pioneer spacecraft. 
The effect of this interaction is a frequency shift that 



is proportional to the distance and the square root of 
the density of the medium in wh ich i t travels. Similarly, 
Didon, Perchoux, and Courtens |l27j | proposed that the 
effect comes from resistance of the spacecraft antennae as 
they transverse the interplanetary dust. This is related 
to more general ideas that an asteroid or comet belt, 
with its associated dust, might cause the effect by gravi- 
tational interactions (see Section FVII El) or resistance to 
dust particles. 

However, these ideas have problems with known prop- 
erties of the interplanetary medium that were outlined in 
Section IVII El In particular, infrared observations rule 
out more than 0.3 Earth mass from Kuiper Belt dust in 
the trans-Neptunian region [osL lo^ . Ulysses and Galileo 
measurements in the inner solar syste m find very few dust 
grains in the 10"^* — 10"^^ kg range |128| . The density 
varies greatly, up and down, within the belt (which pre- 
cludes a constant force) and, in any event, the density is 
not large enough to produce a gravitational acceleration 
on the order of ap l95 l -|97 | . 
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One can also speculate that there is some unknown in- 
teraction of the radio signals with the solar wind. An 
experimental answer could be given with two different 
transmission frequencies. Although the main communi- 
cation link on the Ulysses mission is S-up/X-down mode, 
a small fraction of the data is S-up/S-down. We had 
hoped to utilize this option in further analysis. However, 
using them in our attempt to study a possible frequency 
dependent nature of the anomaly, did not provide any 
useful results. This was in part due to the fact that X- 
band data (about 1.5 % of the whole data available) were 
taken only in the close proximity to the Sun, thus pro- 
hibiting the study of a possible frequency dependence of 
the anomalous acceleration. 



B. Dark matter or modified gravity? 

It is interesting to speculate on the unlikely possibility 
that the origin of the anomalous signal is new physics 
jl29| |. This is true even though the probability is that 
some "standard physics" or some as-yet-unknown sys- 
tematic will be found to explain this "acceleration." The 
first paradigm is obvious. "Is it dark matter or a modifi- 
cation of gravity?" Unfortunately, neither easily works. 

If the cause is dark matter, it is hard to understand. A 
spherically-symmetric distribution of matter which goes 
as p ^ r^^ produces a constant acceleration inside the 
distribution. To produce our anomalous acceleration 
even only out to 50 AU would require the total dark 
matter to be greater than 3 x lO^*Af0. But this is in 
conflict with the accuracy of the ephemeris, which allows 
only of order a few times 10~^Mq of dark matter even 
within the orbit of Uranus . ( A 3- cloud neutrino model 
also did not solve the problem |13(]| |.') 

Contrariwise, the most commonly studied possible 
modification of gr avity (at various scales) is an added 
Yukawa force |13.1| . Then the gravitational potential is 



V{r) = - 



GMm 
(1 + a)r 



1 



~r/A 



(55) 



where a is the new coupling strength relative to Newto- 
nian gravity, and A is the new force's range. Since the 
radial force is — —drV{r) = ma, the power series for 
the acceleration yields an inverse-square term, no inverse- 
r term, then a constant term. Identifying this last term 
as the Pioneer acceleration yields 



ap = 



aai 



2{l + a) X^' 



(56) 



where ai is the Newtonian acceleration at distance ri = 1 
AU. (Out to 65 AU there is no observational evidence of 
an r term in the acceleration.) Eq. H56|) is the solution 
curve; for example, a = — 1 x 10"'^ for A = 200 AU. 

It is also of interest to consider some of the recent pro- 
posals to modify gravity, as alternatives to dark matter 
|l32j| - p'35j |. Consider Milgrom's proposed modification 



of gravity '13^, where the gravitational acceleration of 
a massive body is a cx for some constant Oq ^ a 
and a cx 1/r for ag 3> a. Depending on the value of H, 
the Hubble constant, oq ~ ap! Indeed, as a number of 
people have noted. 



UH = cH ^ 8 X 10^^ cm/s 



(57) 



if H = 82 km/s/Mpc. 

Of course, there are (fundamental and deep) theoreti- 
cal problems if one has a new force of the phenomenolog- 
ical types of those above. Even so, the deep s pace data 
piques our curiosity. In fact, Capozziello et al. |136| | note 
the Pioneer anomaly in their discussion of astrophysical 
structures as manifestations of Yukawa coupling scales. 
This ties into the above discussion. 

However, any universal gravitational explanation for 
the Pioneer effect comes up against a hard experimen- 
tal wall. The anomalous acceleration is too large to 
have gone undetected in planetary orbits, particularly 
for Earth and Mars. NASA's Viking mission provided 
ra dio-rangin g measurements to an accuracy of about 12 
m |l37lll38j| . If a planet experiences a small, anomalous, 
radial acceleration, a a, its orbital radius r is perturbed 

by 



Ar = 



r aA 

UN 



(58) 



where I is the orbital angular momentum per unit mass 
and a^r is the Newtonian acceleration at r. (The right 
value in Eq. H58|l holds in the circular orbit limit.) 

For Earth and Mars, Ar is about —21 km and —76 
km. However, the Viking data determines the difference 
between the Mars and Earth orbital radii to about a 100 
m accuracy, and their sum to an accuracy of about 150 
m. The Pioneer effect is not seen. 

Further, a perturbation in r produces a perturbation 
to the orbital angular velocity of 



Auj = 



2laA 29 aA 



GMp 



(59) 



The determination of the synodic angular velocity (loe — 
lum) is accurate to 7 parts in 10"^^, or to about 5 ms accu- 
racy in synodic period. The only parameter that could 
possibly mask the spacecraft-determined an is {GMq). 
But a large error here would cause inconsistencies with 
the overall planetary ephemeris |^ [Also, there 

would be a pr oblem with the advance of the perihelion 
of Icarus |l39j .] 

We conclude that the Viking ranging data limit any 
unmodeled radial acceleration acting on Earth and Mars 
to no more than 0.1 x 10^* cm/s^. Consequently, if the 
anomalous radial acceleration acting on spinning space- 
craft is gravitational in origin, it is not universal. That 
is, it must affect bodies in the 1000 kg range more than 
bodies of planetary size by a factor of 100 or more. This 
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would be a strange violation of the Principle of Equiva- 
lence |14(]| |. (Similarly, the Aw results rule out the uni- 
versality of the at time-acceleration model. In the age of 
the universe, T, one would have ~ 0.7 T.) 

A new dark matter model was recently proposed by 
Munyaneza and Viollier |l4lj to explain the Pioneer 
anomaly. The dark matter is assumed to be gravitation- 
ally clustered around the Sun in the form of a spherical 
halo of a degenerate gas of heavy neutrinos. However, al- 
though the resulting mass distribution is consistent with 
constraints on the mass excess within the orbits of the 
outer planets previously mentioned, it turns out that the 
model fails to produce a viable mechanism for the de- 
tected anomalous acceleration. 



C. New suggestions stimulated by the Pioneer 
effect 

Due to the fact that the size of the anomalous acceler- 
ation is of order cH , where H is the Hubble constant (see 
Eq. (|57|l '). the Pioneer results have stimulated a number 
of new physics sugge stions. For example, Rosales and 
Sanchez-Gomez 142] propose that ap is due to a local 
curvature in light geodesies in the expanding spacetime 
universe. They argue that the Pioneer effect represents 
a new cosmological Foucault experiment, since the solar 
system coordinates are not true inertial coordinates with 
respect to the expansion of the universe. Therefore, the 
Pioneers are mimicking the role that the rotating Earth 
plays in Foucault's experiment. Therefore, in this picture 
the effect is not a "true physical effect" and a coordinate 
transformation to the co-moving cosmological coordinate 
frame would entirely remove the Pioneer e ffect . 

From a similar viewpoint, Guruprasad 143] finds ac- 
commodation for the constant term while trying to ex- 
plain the annual term as a tidal effect on the physical 
structure of the spacecraft itself. In particular, he sug- 
gests that the deformations of the physical structure of 
the spacecraft (due to external factors such as the effec- 
tive solar and galactic tidal forces) combined with the 
spin of the spacecraft are directly responsible for the de- 
tected annual anomaly. Moreover, he proposes a hypoth- 
esis of the planetary Hubble's flow and suggests that Pio- 
neer's anomaly does not contradict the existing planetary 
data, but supports his new theory of relativistically elas- 
tic space-ti me. 

0stvang |l44j further exploits the fact that the grav- 
itational field of the solar system is not static with re- 
spect to the cosmic expansion. He does note, however, 
that in order to be acceptable, any non-standard expla- 
nation of the effect should follow from a general theoret- 
ical framework. Even so, 0stvang still presents quite a 
radical model. This model advocates the use of an ex- 
panded PPN-framework that includes a direct effect on 
local scales due to the cosmic space-time expansion. 

Belayev |l45j considers a Kaluza-Klein model in 5 di- 
mensions with a time-varying scale factor for the com- 



pactified fifth dimension. His comprehensive analysis led 
to the conclusion that a variation of the physical con- 
stants on a cosmic time scale is responsible for the ap- 
pearance of the anomalous acceleration observed in the 
Pioneer 10/11 trac king data. 

Modanese |l46l | considers the effect of a scale- 
dependent cosmological term in the gravitational ac- 
tion. It turns out that, even in the case of a static 
spherically-symmetric source, the external solution of his 
modified gravitational field equations contains a non- 
Schwartzschild-like component that depends on the size 
of the test particles. He argues that this additional term 
may be relevant to the observed anomaly. 

Mansouri, Nasseri and Khorrami |147| argue that there 
is an effective time variation in the Newtonian gravi- 
tational constant that in turn may be related to the 
anomaly. In particular, they consider the time evolution 
of G in a model universe with variable space dimensions. 
When analyzed in the low energy limit, this theory pro- 
duces a result that may be relevant to the long-range 
acceleration discussed here. A similar analysis was per- 
formed by Sidharth 148], who also discussed cosmologi- 
cal models with a time- varying Newtonian gravitational 
constant. 

Inavov |l49l | suggests that the Pioneer anomaly is pos- 
sibly the manifestation of a superstrong interaction of 
photons with single gravitons that form a dynamical 
background in the solar system. Every gravitating body 
would experience a deceleration effect from such a back- 
ground with a magnitude proportional to Hubble's con- 
stant. Such a deceleration would produce an observable 
effect on a solar system scale. 

All these ideas produce predictions that are close to 
Eq. (|S7|l, but they certainly must be judged against dis- 
cussions in the following two subsections. 

In a different framework. Foot and Volkas jl5d|. sug- 
gest the anomaly can be explained if there is mirror mat- 
ter of mirror dust in the solar system, this could pro- 
duce a drag force and not violate solar-system mass con- 
straints. 

Several scalar-field ideas have also appeared. Mbelek 
and Lachieze-Rey 451] have a model based on a long- 
range scalar field, which also predicts an oscillatory de- 
cHne in ap beyond about 100 AU. This model does ex- 
plain the fact that ap stays approximately constant for 
a long period (recall that Pioneer 10 is now past 70 AU). 
From a similar standpoint Calchi Novati et al. (15Z] dis- 
cuss a weak-limit, scalar-tensor extension to the standard 
gravitational model. However, before any of these pro- 
posals can be seriously considered they must explain the 
precise timing data for millisecond binary pulsars, i.e., 
the gravitational radiation indirectly observed in PSR 
1913-1-16 by Hulse and Taylor |l53l |. Furthermore, there 
should be evidence of a distance-dependent scalar field if 
it is uniformly coupled to ordinary matter. 

Consoli and Siringo |l54j and Consoli |l55l | consider the 
Newtonian regime of gravity to be the long wavelength 
excitation of a scalar condensate from electroweak sym- 
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metry breaking. They speculate that the self-interactions 
of the condensate co uld b e the origins of both Milgrom's 
inertia modification |l82l llSRl ] and also of the Pioneer 
effect. 

Capozziello and Lanibiase |l5fij argue that flavor os- 
cillations of neutrinos in the Brans-Dicke theory of grav- 
ity may produce a quantum mechanical phase shift of 
neutrinos. Such a shift would produce observable effects 
on astrophysical/cosmo logical length and time scales. In 
particular, it results in a variation of the Newtonian grav- 
itational constant and, in the low energy limit, might be 
relevant to our study. 

Motivated by th e work of Mannheim |l33lll57j . Wood 
and Moreau |158| | investigated the theory of conformal 
gravity with dynamical mass generation. They argue 
that the Higgs scalar is a feature of the theory that cannot 
be ignored. In particular, within this framework they find 
one can reproduce the standard gravitational dynamics 
and tests within the solar system, and yet the Higgs fields 
may leave room for the Pioneer effect on small bodies. 

In summary, as highly speculative as all these ideas 
are, it can be seen that at the least the Pioneer anomaly 
is influencing the phenomenological discussion of modern 
gravitational physics and quantum cosmology |159l | . 



D. Phenomenological time models 

Having noted the relationships ap ^ c at oi Eq. H16(l 
and that of Eq. (|57|l . we were motivated to try to 
think of any (purely phenomenological) "time" distor- 
tions that might fortuitously fit the CHASMP Pioneer 
results shown in Figure |H1 In other words, are Eqs. 157|l 
and/or (|16|1 indicating something? Is there any evidence 
that some kind of "time acceleration" is being seen? 

The Galileo and Ulysses spacecraft radio tracking data 
was especially useful. We examined numerous "time" 
models searching for any (possibly radical) solution. It 
was thought that these models would contribute to the 
definition of the different time scales constructed on the 
basis of Eq. ^ and discussed in the Section FlV Bl The 
nomenclature of the standard time scales ^54j-,5_5,] was 
phenomenologically extended in our hope to find a desir- 
able quality of the trajectory solution for the Pioneers. 

In particular we considered: 

i) Drifting Clocks. This model adds a constant accel- 
eration term to the Station Time (ST) clocks, i.e., in the 
ST-UTC (Universal Time Coordinates) time transforma- 
tion. The model may be given as follows: 

AST = STreceived ^ ST^gnt > AST + — ftdocks ' AST^ 

(60) 

where STreceived and STsent are the atomic proper times 
of sending and receiving the signal by a DSN antenna. 
The model fit Doppler well for Pioneer 10, Galileo, and 
Ulysses but failed to model range data for Galileo and 
Ulysses. 



ii) Quadratic Time Augmentation. This model adds a 
quadratic-in-time augmentation to the TAI-ET (Interna- 
tional Atomic Time — Ephemeris Time) time transforma- 
tion, as follows 

1 9 

ET ET+-aET-ET^. (61) 

The model fits Doppler fairly well but range very badly. 

iii) Frequency Drift. This model adds a constant fre- 
quency drift to the reference S-band carrier frequency: 

/i^ fi , "fr.drift • TAI\ 

i^s-band(i) = t-of^H j- (62) 

The model also fits Doppler well but again fits range 
poorly. 

iv) Speed of Gravity. This model adds a "light time" 
delay to the actions of the Sun and planets upon the 
spacecraft: 

Vav = c(l + "^P-s--- ■ ~ . (63) 

The model fits Pioneer 10 and Ulysses well. But the 
Earth flyby of Galileo fit was terrible, with Doppler resid- 
uals as high as 20 Hz. 

All these models were rejected due either to poor fits 
or to inconsistent solutions among spacecraft. 

E. Quadratic in time model 

There was one model of the above type that was espe- 
cially fascinating. This model adds a quadratic in time 
term to the light time as seen by the DSN station. Take 
any labeled time Ta to be 

Ta = - to -> - to + ^Ot (^a ~ ^o) ■ (64) 

Then the light time is 

ATAI TAIj-eceived 1'^lsent ^ 

^ ATAI + ioq^ad • (tAi2,,,i,3, - TAlLnt) • (65) 

It mimics a line of sight acceleration of the spacecraft, 
and could be thought of as an expanding space model. 
Note that Oquad affects only the data. This is in contrast 
to the at of Eq. H16|) that affects both the data and the 
trajectory. 

This model fit both Doppler and range very well. Pi- 
oneers 10 and 11, and Galileo have similar solutions al- 
though Galileo solution is highly correlated with solar 
pressure; however, the range coefficient of the quadratic 
is negative for the Pioneers and Galileo while positive for 
Ulysses. Therefore we originally rejected the model be- 
cause of the opposite signs of the coefficients. But when 
we later appreciated that the Ulysses anomalous accelera- 
tion is dominated by gas leaks (see Section IV D 2|) . which 
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makes the different-sign coefficient of Ulysses meaning- 
less, we reconsidered it. 

The fact that the Pioneer 10 and 11, Galileo, and 
Ulysses are spinning spacecraft whose spin axis are pe- 
riodically adjusted so as to point towards Earth turns 
out to make the quadratic in time model and the con- 
stant spacecraft acceleration model highly correlated and 
therefore very difficult to separate. The quadratic in time 
model produces residuals only slightly 20%) larger 
than the constant spacecraft acceleration model. How- 
ever, when estimated together with no a priori input i.e., 
based only the tracking data, even though the correlation 
between the two models is 0.97, the value Cquad deter- 
mined for the quadratic in time model is zero while the 
value for the constant acceleration model ap remains the 
same as before. 

The orbit determination process clearly prefers the 
constant acceleration model, ap, over that the quadratic 
in time model, Cquad of Eq. H65() . This implies that a 
real acceleration is being observed and not a pseudo ac- 
celeration. We have not rejected this model as it may 
be too simple in that the motions of the spacecraft and 
the Earth may need to be included to produce a true 
expanding space model. Even so, the numerical relation- 
ship between the Hubble constant and ap, which many 
people have observed (cf. Section FXI C|l . remains an in- 
teresting conjecture. 



XII. CONCLUSIONS 

In this paper we have discussed the equipment, theo- 
retical models, and data analysis techniques involved in 
obtaining the anomalous Pioneer acceleration ap. We 
have also reviewed the possible systematic errors that 
could explain this effect. These included computational 
errors as well as experimental systematics, from systems 
both external to and internal to the spacecraft. Thus, 
based on further data for the Pioneer 10 orbit determi- 
nation (the extended data spans 3 January 1987 to 22 
July 1998) and more detailed studies of all the systemat- 
ics, we can now give a total error budget for our analysis 
and a latest result of ap = (8.74 ± 1.33) x 10"^ cm/s^ 

This investigation was possible because modern radio 
tracking techniques have provided us with the means 
to investigate gravitational interactions to an accuracy 
never before possible. With these techniques, relativis- 
tic solar-system celestial mechanical experiments using 
the planets and interplanetary spacecraft provide critical 
new information. 

Our investigation has emphasized that effects that pre- 
viously thought to be insignificant, such as rejected ther- 
mal radiation or mass expulsion, are now within (or near) 
one order of magnitude of possible mission requirements. 
This has unexpectedly emphasized the need to carefully 
understand all systematics to this level. 

In projects proposed for the near future, such as a 
Doppler measurement of the solar gravitational deflection 



using the Cassin i spa cecraft |l60j and the Space Interfer- 
ometry Mission |l6l| . navigation requirements are more 
stringent than those for current spacecraft. Therefore, all 
the effects we have discussed will have to be well-modeled 
in order to obtain sufficiently good trajectory solutions. 
That is, a better understanding of the nature of these 
extra small forces will be needed to achieve the stringent 
navigation requirements for these missions. 

Currently, we find no mechanism or theory that ex- 
plains the anomalous acceleration. What we can say with 
some confidence is that the anomalous acceleration is a 
line of sight constant acceleration of the spacecraft to- 
ward the Sun [t^- Even though fits to the Pioneers ap- 
pear to match the noise level of the data, in reality the 
fit levels are as much as 50 times above the fundamental 
noise limit of the data. Until more is known, we must ad- 
mit that the most likely cause of this effect is an unknown 
systematic. (We ourselves are divided as to whether "gas 
leaks" or "heat" is this "most likely cause.") 

The arguments for "gas leaks" are: i) All spacecraft 
experience a gas leakage at some level, ii) There is enough 
gas available to cause the effect, iii) Gas leaks require no 
new physics. However, iv) it is unlikely that the two 
Pioneer spacecraft would have gas leaks at similar rates, 
over the entire data interval, especially when the valves 
have been used for so many maneuvers. [Recall also that 
one of the Pioneer 11 thrusters became inoperative soon 
after launch. (See Section rilBD ] v) Most importantly, 
it would require that these gas leaks be precisely pointed 
towards the front JJj of the spacecraft so as not to cause 
a large spin-rate changes. But vi) it could still be true 
anyway. 

The main arguments for "heat" are: i) There is so 
much heat available that a small amount of the total 
could cause the effect, ii) In deep space the spacecraft 
will be in approximate thermal equilibrium. The heat 
should then be emitted at an approximately constant 
rate, deviating from a constant only because of the slow 
exponential decay of the Plutonium heat source. It is 
hard to resist the notion that this heat somehow must 
be the origin of the effect. However, iii) there is no solid 
explanation in hand as to how a specific heat mechanism 
could work. Further, iv) the decrease in the heat supply 
over time should have been seen by now. 

Further experiment and analysis is obviously needed 
to resolve this problem. 

On the Pioneer 10 experimental front, there now exists 
data up to July 2000. Further, there exists archived high- 
rate data from 1978 to the beginning of our data arc in 
Jan 1987 that was not used in this analysis. Because this 
early data originated when the Pioneers were much closer 
in to the Sun, greater effort would be needed to perform 
the data analyses and to model the systematics. 

As Pioneer 10 continues to recede into interstellar 
space, its signal is becoming dimmer. Even now, the 
return signal is hard to detect with the largest DSN an- 
tenna. However, with appropriate instrumentation, the 
305-meter antenna of the Arecibo Observatory in Puerto 
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Rico will be able to detect Pioneer's signal for a longer 
time. If contact with Pioneer 10 can be maintained with 
conscan maneuvers, such further extended data would be 
very useful, since the spacecraft is now so far from the 
Sun. 

Other spacecraft can also be used in the study oi ap. 
The radio Doppler and range data from the Cassini mis- 
sion could offer a potential contribution. This mission 
was launched on 15 October 1997. The potential data 
arc will be the cruise phase from after the Jupiter flyby 
(30 December 2000) to the vicinity of Saturn (just before 
the Huygens probe release) in July 2004. Even though 
the Cassini spacecraft is in three-axis-stabilization mode, 
using on-board active thrusters, it was built with very so- 
phisticated radio-tracking capabilities, with X-band be- 
ing the main navigation frequency. (There will also be S- 
and K-band links.) Further, during much of the cruise 
phase, reaction wheels will be used for stabilization in- 
stead of thrusters. Their use will aid relativity exper- 
iments at solar conjunction and gravitational wave ex- 
periments at solar opposition. (Observe, however, that 
the relatively large systematic from the close in Cassini 
RTGs will have to be accounted for.) 

Therefore, Cassini could yield important orbit data, in- 
dependent of the Pioneer hyperbolic-orbit data. A simi- 
lar opportunity may exist, out of the plane of the ecliptic, 
from the proposed Solar Probe mission. Under considera- 
tion is a low-mass module to be ejected during solar flyby. 
On a longer time scale, the reconsidered Pluto/Kuiper 
mission (with arrival at Pluto around 2029) could even- 
tually provide high-quality data from very deep space. 

All these missions might help test our current mod- 
els of precision navigation and also provide a new test 
for the anomalous ap. In particular, we anticipate that, 
given our analysis of the Pioneers, in the future preci- 
sion orbital analysis may concentrate more on systemat- 
ics. That is, data/systematic modeling analysis may be 
assigned more importance relative to the astronomical 
modeling techniques people have concentrated on for the 
past 40 years )l6^ - |l64| . 

Finally, we observe that if no convincing explanation 
is to be obtained, the possibility remains that the effect 
is real. It could even be related to cosmological quanti- 
ties, as ha]s been intimated. [See Eq. and Sec. IXII 
especially the text around Eq. (|57|l .] This possibility ne- 
cessitates a cautionary note on phenomenology: At this 
point in time, with the limited results available, there is 
a phcnomcnological equivalence between the ap and at 
points of view. But somehow, the choice one makes af- 
fects one's outlook and direction of attack. If one has to 
consider new physics one should be open to both points 
of view. In the unlikely event that there is new physics, 
one does not want to miss it because one had the wrong 
mind set. 
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APPENDIX A: APPENDIX 

In Tabic IIIII we give the hyperbolic orbital parame- 
ters for Pioneer 10 and Pioneer 11 at epoch 1 January 
1987, 01:00:00 UTC. The semi-major axis is a, e is the 
eccentricity, / is the inclination, 11 is the longitude of the 
ascending node, lo is the argument of the perihelion, Mg 
is the mean anomaly, /o is the true anomaly at epoch, 
and To is the heliocentric radius at the epoch. The direc- 
tion cosines of the spacecraft position for the axes used 
are (a, (3, 7). These direction cosines and angles are re- 
ferred to the mean equator and equinox of J2000. The 
ecliptic longitude £0 f^nd latitude 69 a-re also listed for an 
obliquity of 23° 26 21. 4119. The numbers in parentheses 
denote realistic standard errors in the last digits. 



Parameter 


Pioneer 10 


Pioneer 11 


a [km] 


-1033394633(4) 


-1218489295(133) 


e 


1.733593601(88) 


2.147933251(282) 


I [Deg] 


26.2488696(24) 


9.4685573(140) 


n [Deg] 


-3.3757430(256) 


35.5703012(799) 


u [Deg] 


-38.1163776(231) 


-221.2840619(773) 


Mo [Deg] 


259.2519477(12) 


109.8717438(231) 


/o [Deg] 


112.1548376(3) 


81.5877236(50) 


ro [km] 


5985144906(22) 


3350363070(598) 


Q 


0.3252905546(4) 


-0.2491819783(41) 


P 


0.8446147582(66) 


-0.9625930916(22) 


7 


0.4252199023(133) 


-0.1064090300(223) 


4 [Deg] 


70.98784378(2) 


-105.06917250(31) 


60 [Deg] 


3.10485024(85) 


16.57492890(127) 



TABLE III: Orbital parameters for Pioneer 10 and Pioneer 
11 at epoch 1 January 1987, 01:00:00 UTC. 
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We give the errata for this paper in both the arXiv |gr-qc/0104064| 
[Phys. Rev. D 65, 082004 (2002)] formats. 



version 4] and published 



in 
o 
o 



PACS numbers: 04.80.-y, 95.10.Eg, 95.55.Pe 



Sec. I-C, p. (4), para. 4, hue 6 (published version 
only): 

Delete "load sharing." 



Eq. (29): Change 1.24 x 10"" to 4.4 x 10 
3 lines down: Change 10~^ to 10~^. 
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Sec. IV-F, 3rd para, from end, p. (16), lines 3-4 
, (published version only) 

■ Delete "the" in front of "JPL's" and one of the two 
If. ' successive "from"s. 



Table II, p. 43 (40), item 1. (b): 
10"^ should be 10"^ 



Sec. XI-C, p. 44 (41), 4 lines below Eq. (57): 
10-3 should be lO^^. 




"Inavov" should be "Ivanov. 



Sec. XI-C, p. 45 (42), para. 7, line 1: 




Line 2: Change /© to /©/c and to w^. 



Eq. (29): Change w^/c to w^. 



Sec. Vn-B, p. 29 (27): 



■Mashoon" should be "Mashhoon. 



Ref. [20], p. 49 (45): 




